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ABSTRACT 


In  part  I,  Schrodinger ’ s  spin  one  exchange  interaction, 

2 

P j_ j  =  (~ i*i§ j )  ’*■  -1,  is  introduced  to  form  a  Hamiltonian  for  a 

nearest  neighbour  lattice  model  of  ferromagnetism.  The  physical  rele¬ 
vance  of  the  Hamiltonian  is  discussed.  Group  theoretic  methods  to¬ 
gether  with  the  cluster  expansion  technique  for  obtaining  the  high 
temperature  partition  function  and  initial  susceptibility  series  are 
introduced.  The  first  eight  terms  of  the  susceptibility  series  are 
obtained  for  the  face-centred  cubic  lattice.  Analysis  by  the  ratio  and 
Pade  approximant  methods  yields  the  critical  temperature  and  the  criti¬ 
cal  index  y  which  is  found  to  be  very  nearly  5/^  and  thus  certainly 
different  from  the  spin  one  Heisenberg  model  index. 

In  part  II,  the  Kubo  linear  response  theory  is  applied  to  the 
quantum  mechanical  Ising  model  to  obtain  a  general  expression  which  is 
valid  for  all  regular  lattices  for  the  perpendicular  response  function 
in  terms  of  spin  correlations.  For  the  honeycomb  and  plane  square  lat¬ 
tices,  exact  expressions  are  found  for  the  frequency  dependent  initial 
perpendicular  susceptibility  Xj(w,T)  valid  for  all  temperatures.  A 
discrepancy  arises  between  the  zero  frequency  result  for  the  plane 
square  lattice  and  Fisher’s  static  calculation.  This  is  resolved  by 
using  an  extended  calculation  which  explicitly  illustrates  the  differ¬ 
ence  between  isothermal  and  adiabatic  susceptibilities. 
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Finally  the  magnetic  inelastic  scattering  of  neutrons  from  a  spin 
z  Ising  system  is  considered.  For  the  plane  square  lattice,  the  exact 
energy  distribution  of  the  scattering  for  all  temperatures  is  obtained 
and  found  to  vary  most  markedly  in  the  critical  region.  Some  qualitative 
projections  are  also  made  concerning  three  dimensional  systems. 
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CHAPTER  I 
INTRODUCTION 


A.  Survey  of  Physical  Systems 

The  study  of  phase  transitions  and  critical  points  cuts  across 
many  of  the  familiar  boundaries  of  the  sciences.  It  attracts  the  inter¬ 
est  of  the  solid-state  physicist,  the  low-temperature  physicist,  the 
physical  chemist,  the  chemical  engineer  and  the  metallurgist.  One  rea¬ 
son  lies  in  the  large  variety  of  physical  systems  that  exhibit  critical 
phenomena.  They  include  simple  fluids,  ferromagnets,  binary  solutions, 
antiferromagnets  and  binary  alloys.  It  is  one  of  the  purposes  of  these 
opening  paragraphs  to  point  out  that  despite  the  variety  of  systems, 
their  quantitative  behavior  displays  remarkable  similarities.  It  is 
one  of  the  challenges  of  this  subject  to  reconcile  the  variety  with  the 
similarity  and,  as  we  proceed, we  shall  see  to  what  extent  this  has  been 
achieved. 

It  should  be  pointed  out  that  such  systems  as  ferroelectrics , 
superconductors  and  superfluids  have  been  deliberately  left  out.  They 
stand  apart  in  the  first  case  because  of  the  dominant  nature  of  the 
long  range  coulomb  forces  and  in  the  latter  because  of  the  essentially 
quantum-mechanical  origin  of  the  phenomena. 

In  figure  1.1  a  number  of  (p,v)  isotherms  for  a  simple  fluid  are 


Fig.  1.1 


Typical  isotherms  for  a  simple  fluid. 
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sketched.  It  is  an  experimental  fact  that  for  temperatures  greater  than 

the  critical  temperature,  T  ,  one  can  pass  by  compression  from  the  gas- 

eous  state  to  the  liquid  state  without  any  anomalies  in  the  derivatives 

of  the  free  energy.  However  for  temperatures  close  to  Ttl  there  exists 

c 

a  well-defined  maximum  in  the  isothermal  compressibility, 

=  “  v  ^ap^T  »  (1.1) 

which  occurs  at  the  critical  density  p,  =  l/v  .  As  T  is  approached, 
this  maximum  increases  in  intensity  (corresponding  to  the  flattening  of 
the  isotherms  in  figure  1.1)  and  we  characterize  the  critical  point  by 
the  statement, 

Kp(p  =  pc,  T—*Tc+)—?oo  #  (1.2) 

In  addition,  the  critical  point  is  characterized  by  the  'closure* 
of  the  coexistence  curve  (see  figure  l.l)  i.e.  by 

(pL  -  PG)-^°  •  (T->TC.)  (1.3) 

In  a  very  careful  and  accurate  experiment,  Weinberger  and  Schneider 
(1952)  measured  this  density  discontinuity  for  Xenon  and  obtained  the 
result, 

(pL  -  pG)  =  A(TC  -  TW  .  <T-»T8-) 


(1.4) 
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with 


0  =  0.345  ±  0.015  • 


(1.5) 


There  is  an  understandable  tendency  in  this  field  to  anticipate  that  the 
critical  indices,  of  which  1  0'  is  our  first  example,  should  equal  ra¬ 
tional  numbers.  In  that  spirit  (1.5)  strongly  suggests  that  a  1 l/3 ’  rd 
power  law  describes  the  behavior  of  the  coexistence  curve. 

Finallytwe  add  that  the  critical  point  is  also  associated  with  a 
striking  thermal  anomaly.  In  recent  experiments  on  the  specific  heat 
of  argon,  Bagatskii, Voronel  and  Gusak  (1963)  reach  the  con¬ 

clusion  that 


(1.6) 


The  above  experimental  facts  indicate  the  critical  region  is  of 


primary  importance  and  interest.  This  is  more  than  substantiated  when  we 
consider  ferromagnets  which  would  appear  to  bear  very  little  similarity 
in  their  physical  properties  to  simple  fluids.  A  ferromagnet  is  char¬ 
acterized  by  a  spontaneous  or  residual  magnetization  which  persists  at 
low  temperatures  even  when  the  magnetizing  field  is  reduced  to  zero,i.e. 


(1.7) 


where  M(T,H)  is  the  equilibrium  magnetization  in  the  field  H. 

The  critical  temperature  (known  as  the  Curie  point  in  this  context) 


1.4 


is  associated  with  the  following  behavior,  (see  figure  1.2) 


M0(T )  — >  0  (T — >  T  -) 

\J  c 

Mn(T)  =0  T  >  T 
0  c 


In  addition,  the  initial  susceptibility  ,  given  by 


Z0(T)  =  (§f) 


T  ,H  =  0 


(1.8) 


(1.9) 


is  positive  and  as  T — jT  +,  it  diverges  to  infinity  (see  figure  1.2). 

c 

Reflection  of  the  magnetization  curve  in  figure  1.2  about  the  Mq = 0  axis 
and  anticlockwise  rotation  through  an  angle  of  90°  suggests  an  analogy 
with  the  gas-liquid  coexistence  curve  in  figure  1.1.  In  figure  1.3  an 
ideal  magnetization  curve  for  a  ferromagnet  is  drawn  and  again  a  compari¬ 
son  with  figure  1.1  suggests  that  the  susceptibility  of  the  ferromagnet 
and  the  compressibility  of  the  gas-liquid  system  are  analogous  variables. 
One  difference,  however,  is  that  the  discontinuity  in  the  magnetization, 
i.e.  a  two  phase  region,  occurs  only  for  zero  field  while  p  -  pG  is 
non-zero  for  all  p  <  p,.  (In  both  instances,  of  course,  T  must  be  less 

than  T  ). 
c 

That  the  above  comparison  is  somewhat  more  than  an  analogy  is  de¬ 
monstrated  by  the  experiments  of  Heller  and  Benedek  (1965).  Using  mag¬ 
netic  resonance  techniques,  they  measured  the  spontaneous  magnetization 


Spontaneous  magnetization  and  initial  susceptibility 
of  a  f erromagnet  near  its  Curie  point. 
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Fig  1.3 


Ideal  magnetization  curve  of  a  ferromagnet 


1.4b 
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of  the  insulating  ferromagnet  EuS  and  found 

Mq(T)  =  A(TC  -  T)P  (T  — >  T^-)  (1.10) 

where  3  =  0.33  i  0.01.  A  one  third  law  seems  also  to  apply  to  ferro- 
magnets. 

A  binary  fluid  mixture  of  two  components  A  and  B,  in  which  A  mol¬ 
ecules  attract  A  and  B  attract  B,  will  undergo  a  phase  separation  at 
low  temperatures  into  an  A- rich  phase  and  a  B-rich  phase.  According  to 
J.S,  Rowlinson  (1959) »  the  coexistence  curve,  representing  in  this  case 
the  difference  in  molar  fractions  of  the  two  phases,  again  follows  a 
one  third  power  law  to  a  good  approximation. 

The  three  systems  we  have  just  considered  are  examples  in  which 
•like*  attracts  ’like'.  The  molecules  of  the  gas  attract  one  another, 
the  'up'  spins  of  the  ferromagnet  attract  up  spins  and  the  A  molecules 
of  the  binary  fluid  attract  other  A  molecules.  If  we  consider  cases  in 
which  ’like1  attracts  "unlike"  e.g.  the  antiferromagnet,  we  again  find 
significant  similarities  in  behavior.  Measurements  by  Heller  and  Bene- 
dek  (1962)  of  the  'sublattice  magnetization'  of  the  antiferromagnet  MnF^ 
indicate  once  more  that  a  one  third  power  law  describes  its  temperature 
dependence. 

This  discussion  could  continue  with  other  examples  to  further 
emphasize  the  remarkable  similarity  in  behavior  of  these  seemingly  dis¬ 
tinct  physical  systems.  However,  we  shall  conclude  this  section  by  con- 
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sidering  what  all  this  evidence  suggests. 

With  the  possible  exception  of  Bose-Einstein  condensation,  there 
is  no  theory  of  a  phase  transition  that  does  not  relate  the  occurence  of 
the  critical  point  to  the  existence  of  some  intermolecular  interaction. 
It  seems  clear  the  interaction  is  a  necessary  prerequisite  for  a  phase 
transition  to  take  place.  One  of  the  most  important  questions  in  this 
field  is  how  important  are  the  details  of  the  interaction  in  determining 
the  precise  behavior  of  coexistence  curves,  susceptibilities  etc.  in  the 
critical  region.  The  experimental  evidence  we  have  just  reviewed  sug¬ 
gests  that  many  of  the  details  are  not  very  important.  On  the  other 
hand,  there  must  exist  some  common  feature  in  the  interactions  present 
in  all  these  systems  that  is  principally  responsible  for  the  observed 
behavior. 

In  many  respects,  the  models  that  theoreticians  have  examined, 
particularly  the  Ising  model,  are  attempts  to  remove  all  irrelevant  de¬ 
tail  from  the  interactions  and  calculate  the  consequences  of  the  one  im¬ 
portant  feature  that  remains. 

We  shall  return  to  this  question  later.  But  before  we  become 
embroiled  in  the  consideration  of  these  details  we  shall  first  give  a 
qualitative  description  of  the  critical  region  and  why  it  occurs.  Re¬ 
search  in  this  field  is  dominated  by  mathematical  and  computational  dif¬ 
ficulties  as  later  chapters  in  this  thesis  will  show  and  it  is  important 
to  have  a  reasonably  clear  physical  picture  of  critical  phenomena  before 
embarking  on  detailed  studies. 


. 
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B.  Description  of  the  Critical  Region 

Because  of  the  similarities  indicated  in  the  last  section  between 
the  various  physical  systems  exhibiting  phase  transitions,  it  is  possible 
to  give  a  qualitative  description  in  a  number  of  different  contexts,  all 
of  which  would  be  equivalent.  One  could  use  the  ’language*  of  simple 
fluids,  of  homogeneous  binary  mixtures  or,  as  we  shall  choose  in  this 
instance,  of  ferromagnetism  and  the  essential  ideas  would  be  the  same. 
Personal  preference  is  dictated  largely  by  one’s  own  background  and  ex¬ 
perience. 

We  shall  consider  a  simple  picture  of  a  ferromagnet  in  which  a 
large  number  of  ’spins’  are  located  at  the  sites  of  a  regular  lattice. 

At  zero  temperature,  all  the  spins  are  aligned  in  the  same  direction 
(throughout  this  thesis  we  are  referring  to  a  single  domain)  giving  rise 
to  the  spontaneous  magnetization.  As  the  temperature  rises,  excited 
states  of  the  system  are  populated  in  which  a  gradually  increasing  num¬ 
ber  of  the  spins  are  reversed  in  direction  (we  are  assuming  in  this  sim¬ 
plified  description  that  the  spins  can  occupy  only  two  antiparallel  ori¬ 
entations).  Physical  intuition  perhaps  suggests  that  the  trend  from  all 
aligned  spins  to  a  random  mixture  of  ’up’  spins  and  ’down’  spins  should 
have  a  smooth  temperature  dependence  in  which  the  final  random  state  is 
reached  in  the  limit  of  infinite  temperature.  Certainly  from  a  theoreti¬ 
cal  standpoint  this  is  what  happens  in  any  finite  system.  However,  this 
is  not  what  takes  place  in  real  ferromagnets  or  in  theoretically  infinite 
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systems  (i.e.  N,V— *«>  such  that  N/V  =  const.). 

A  convenient  description  of  what  we  believe  takes  place  as  the 
temperature  increases  can  be  given  in  terms  of  the  long  range  order 
parameter.  This  is  a  measure  of  the  tendency  of  any  spin  to  align  it- 
self  in  the  same  direction  as  another  spin  in  the  limit  that  the  dis¬ 
tance  separating  the  spins  becomes  infinite.  At  zero  temperature,  the 
long  range  order  is  unity.  As  the  temperature  rises,  it  begins  to  de¬ 
crease,  reflecting  the  increasing  significance  of  the  disrupting  effect 
of  the  thermal  energy  of  the  spins.  However,  rather  than  gradually  de¬ 
crease  to  zero  in  the  infinite  temperature  limit,  the  correlation  be¬ 
tween  the  two  infinitely  separated  spins  becomes  abruptly  zero  at  a 
finite  critical  temperature.  The  spontaneous  magnetisation  per  spin  is 
often  defined  as  proportional  to  the  square  root  of  the  long  range  order 
and.  so  consultation  of  figure  1,2  indicates  how  the  latter  varies  from 
unity  at  zero  temperature  to  zero  at  the  critical  temperature. 

We  can  regard  this  whole  process  as  a  competition  between  the 
ferromagnetic  interaction  and  the  thermal  energy  of  the  spins.  As  the 
temperature  rises  the  thermal  energy  increases  and  at  the  critical 
point  it  finally  and  abruptly  overcomes  the  ferromagnetic  interaction 
tending  to  align  the  infinitely  separated  spins. 

Let  us  now  approach  the  critical  point  from  the  infinite  tempera¬ 
ture  limit.  There,  the  spins  are  randomly  orientated  and  behave  inde¬ 
pendently.  Thermal  energy  completely  dominates  their  motion  and  we  say 
the  ferromagnet  is  behaving  as  a  paramagnet.  However,  as  the  temperature 
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is  lowered,  the  spin-spin  interaction  begins  to  make  itself  felt  and 
short  range  order  gradually  becomes  significant.  By  this  is  meant  there 
exists  some  degree  of  correlation  between  spins  which  are  nearest  neigh¬ 
bours  on  the  lattice  or  not  too  far  separated.  If  one  spin  points  in  a 
given  direction  there  is  a  probability  greater  than  one  half  that  its 
near  neighbours  will  be  aligned  in  the  same  direction. 

As  the  temperature  continues  to  drop,  the  correlation  length,  the 
maximum  separation  of  appreciably  correlated  spins,  increases  until  in 
the  critical  region  it  reaches  macroscopic  dimensions  becoming  infinite 
at  the  critical  temperature.  The  long  range  order  has  set  in.  This  in¬ 
creasing  correlation  between  the  spins  is  reflected  macroscopically  in  a 
rapidly  increasing  magnetic  susceptibility  which  becomes  infinite  at  the 
critical  point.  That  the  spins  should  tend  to  cooperate  so  rapidly  that 
even  at  a  non-zero  temperature  two  spins  infinitely  far  apart  are  ’aware* 
of  each  other  is  an  indication  of  how  significant  is  the  interaction  be¬ 
tween  the  spins.  Clearly  no  non-zero  critical  temperature  could  occur 
in  the  absence  of  this  interaction.  Such  phenomena  are  often  referred 
to  as  ’cooperative'  phenomena  and  we  can  see  how  descriptive  the  term  is. 

We  hope  by  calculating  the  critical  behavior  of  the  susceptibility 
with  various  types  of  interaction  to  determine  how  much  this  'spreading* 
of  the  correlation  throughout  the  spin  system  depends  on  the  details  of 
the  interaction.  The  magnetic  susceptibility  is  therefore  an  important 
'mirror'  of  the  critical  behavior  and  much  of  this  thesis  is  concerned 


with  its  calculation  for  certain  models. 
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Finally,  another  interesting  feature  of  the  critical  region  is 
the  existence  of  large  scale  fluctuations.  Near  the  critical  temperature, 
there  occur  correlated  deviations  of  well-separated  spins  from  their 
equilibrium  orientations.  These  fluctuations  are  not  only  spread  over 
large  regions  of  the  ferromagnet  but  persist  over  long  periods  of  time. 
They  arise  because  of  the  delicate  balance  between  the  thermal  energy  and 
the  spin-spin  interaction  that  exists  in  the  critical  region.  In  fluid 
systems  they  give  rise  to  critical  opalescence  and  the  analogous  phenome¬ 
non  in  magnetic  systems  is  the  critical  scattering  of  neutrons.  There 
has  been  much  theoretical  and  experimental  activity  in  this  area  in  recent 

4 

years  (Green  and  Sengers  1965). 

C.  Brief  Historical  Survey 

One  of  the  standard  treatments  of  many- body  systems  and  phase  tran¬ 
sitions  which  has  existed  for  at  least  sixty  years  and  is  still  used  at 
the  present  time  is  the  ’mean  field’  or  'internal  field'  approach.  It 
is  an  attempt  to  replace  the  perhaps  complicated  pairwise  interactions 
betxgeen  the  particles  (or  spins)  by  a  uniform  internal  field  with  which 
each  particle  interacts  directly.  This  field  is  calculated  self-consist¬ 
ent  ly  so  that  the  average  interaction  of  each  particle  with  its  neigh¬ 
bours  is  preserved.  It  represents  the  first  attempt  to  theoretically 
predict  a  phase  transition.  The  Weiss  theory  of  ferromagnetism  (1907), 
the  Van  der  Waals'  equation  of  state  (1873)  and  the  Bragg-Williams  approx¬ 
imation  for  binary  alloys  (193^)  are  all  examples  of  this  approach  applied 
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to  different  physical  systems. 

Such  theories  do  succeed  in  predicting  the  existence  of  a  criti¬ 
cal  point  but  fail  in  giving  the  observed  behavior  of  the  thermodynamic 
variables  in  the  critical  region.  For  instance,  they  predict  one  half 
power  laws  for  spontaneous  magnetization  and  coexistence  curves  which 
disagree  with  experimental  results.  For  compressibility  and  suscepti¬ 
bility,  they  give  the  famous  Curie-Weiss  law,  namely, 


Kp  = 


B- 


T-T 


=  X0(T)  , 


(T 


’Tc±) 


(1.11) 


which  is  again  not  in  accord  with  experiment  for  temperatures  close  to 
critical.  The  essential  reason  for  their  failure  is  that  the  very  con¬ 
cept  of  'mean  field*  neglects  consideration  of  short-range  and  long- 
range  order  which  play  such  an  important  role  in  determining  the  char¬ 
acteristics  of  the  critical  region.  It  should  be  mentioned  that, in  the 
context  of  binary  alloys,  Bethe  (1935)  attempted  to  improve  on  the  Bragg- 
Williams  approximation  by  taking  nearest  neighbour  correlations  into  ac¬ 
count  but  the  improvement  in  the  results  was  very  limited.  As  one  might 

expect,  he  did  obtain  a  ’tail*  to  the  specific  heat  curve  above  T  which 

c 

is  not  predicted  by  'mean-field*  arguments.  However,  the  finite  discon¬ 
tinuity  in  the  specific  heat  at  T  ,  characteristic  of  the  B.-W.  approx- 

V 

imation,  remains. 

Practically  all  modern  theories  of  critical  phenomena  that  have 


. 
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had  any  success  in  comparison  with  experiment  have  been  based  on  a  very 
small  number  of  models.  The  oldest  and  still  most  popular  is  the  Ising 
model  introduced  and  solved  exactly  in  one  dimension  by  Ising  (1925). 

It  can  be  characterized  by  the  Hamiltonian, 


=  -J  , 

(ij)  J 


(1.12) 


where  i  and  j  label  the  sites  of  the  lattice,  and  a  is  a  two-valued 
variable  (usually  chosen  to  be  +1  or  -l).  We  shall  assume  in  all  our 
considerations  that  the  summation  in  (1.12)  is  over  nearest  neighbour 
pairs  of  sites.  This  is  in  direct  contrast  to  the  mean  field  calcula¬ 
tions  which  are  equivalent  to  a  system  with  infinitely  long  range  but 
infinitely  weak  interactions.  The  Ising  model  will  be  discussed  in 
more  detail  in  the  next  section. 

Ising' s  one  dimensional  solution  did  not  exhibit  a  non-zero 
critical  temperature,  but  some  years  later  Peierls  (1936)  predicted  a 
phase  transition  could  take  place  in  Ising  systems  of  dimension  greater 
than  one.  Then  Kramers  and  Wannier  (1941)  deduced  the  location  of  the 
critical  temperature  for  the  case  of  the  simple  quadratic  lattice  in 
zero  magnetic  field.  Finally,  Onsager  (1944)  solved  the  problem  ex¬ 
actly  on  this  lattice  and  obtained  a  logarithmic  singularity  of  the 
specific  heat.  It  is  difficult  to  exaggerate  the  importance  of  his  sol¬ 
ution.  Apart  from  later  solutions  on  other  two-dimensional  lattices,  it 
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stands  alone  as  an  exact  solution  of  a  many- body  system  which  exhibits 
a  critical  phase  transition.  It  proved  conclusively  that  such  critical 
behavior  could  be  predicted  by  a  model  which  included  only  nearest  neigh¬ 
bour  interactions. 

Since  then,  a  great  amount  of  effort  has  gone  into  the  search  for 
a  similar  exact  solution  of  the  three-dimensional  Ising  model  and  has 
met  with  singularly  little  success.  Adopting  a  different  approach  and 
developing  what  has  now  become  known  as  the  'exact*  series  expansion 
technique,  the  King's  College  group  in  London  have,  however,  achieved 
great  success  in  analysing  the  behavior  of  the  Ising  model  and  others 
near  their  critical  temperatures.  This  approach  has  been  very  compre¬ 
hensively  reviewed  by  Domb  (I960).  Part  I  of  this  thesis  is  an  example 
of  an  exact  series  expansion  calculation.  These  series,  be  they  high 
temperature  or  low  temperature,  are  perturbation  expansions  about  in¬ 
finite  or  zero  temperature  in  appropriate  variables.  A  finite  number 
of  the  coefficients  (as  large  as  can  be  calculated  within  reasonable 
time  limits)  are  found  exactly  and  from  the  truncated  series,  it  turns 
out  to  be  possible  to  make  very  accurate  predictions  of  the  behavior  of 
these  models  near  T  .  Further  details  will  be  given  in  later  chapters. 

Another  important  landmark, particularly  with  regard  to  ferromag¬ 
netism, was  the  introduction  by  Heisenberg  (1928)  of  what  is  now  known 
as  the  Heisenberg  model  of  ferromagnetism.  It  is  described  by  the 
Hamiltonian  , 


W  =  -J 


£ 

(ij) 


O.  ,C7. 


(1.13) 
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which  differs  from  (1.12)  in  that  the  a' s  are  not  scalars  but  quantum, 
mechanical  spin  operators.  It  is  considered  a  more  realistic  model  of  a 
ferromagnet  than  the  Ising  model  but  it  turns  out  to  be  much  harder  to 
solve.  Indeed,  no  exact  solution  in  any  dimension  has  been  obtained  and 
only  approximate  techniques  such  as  exact  series  expansions  for  high 
temperatures  and  spin-wave  theory  for  low  temperatures  have  been  success¬ 
ful  in  predicting  the  behavior  of  this  model. 

D.  Discussion  of  Models 

In  setting  up  a  theory  of  a  many- body  system,  some  may  feel  that 
the  ultimate  goal  is  to  solve  the  many-body  Schrodinger  equation  where 
all  contributions  to  the  Hamiltonian  have  been  included.  Such  a  task  is 
not  only  impossible  but  quite  useless.  Once  the  basic  principles  of 
wave-mechanics  have  been  impressively  verified  by  the  results  of  an  exact 
calculation  on  some  very  simple  system  such  as  the  hydrogen  atom,  it 
serves  little  purpose  to  seek  the  same  degree  of  accuracy  on  more  and 
more  complicated  systems.  If  that  were  the  only  task  of  theoretical  phys¬ 
ics,  it  would  indeed  be  a  dull  and  pedestrian  discipline.  Rather,  when 
confronted,  by  a  complicated  many  body  system,  the  physicist  must  use  his 
intuition  plus  admittedly  certain  compulsions  arising  from  inevitable 
mathematical  limitations  to  simplify  the  problem.  He  ignores  the  details 
he  feels  are  of  little  relevance  to  the  phenomena  he  wishes  to  predict 
and  exaggerates  those  of  primary  importance.  Comparison  of  the  results 
of  his  theory  with  experiment  will  ultimately  tell  him  how  successful 
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his  1  pruning 1  has  been.  In  modern  physics,  our  aim  is  not  to  produce  an 
exact  copy  on  paper  of  nature  and  the  above  procedure  is  the  only  means 
by  which  we  can  improve  our  understanding. 

The  commonly  used  models  of  critical  phenomena,  namely  the  Ising 
model  and  the  Heisenberg  model, can  be  regarded  as  rather  successful  exam¬ 
ples  of  this  procedure.  The  Ising  model,  in  particular,  is  a  case  where 
almost  the  maximum  amount  of  simplification  has  taken  place.  It  is  dif¬ 
ficult  to  think  of  any  further  modifications  one  could  make  to  (1.12)  to 
simplify  it  further  that  does  not  render  the  model  quite  trivial  and  un¬ 
interesting.  One  frequently  reads  of  how  artificial  the  Ising  model  is 
but  it  is  not  easy  to  reconcile  this  point  of  view  with  the  great  ad¬ 
vances  in  our  understanding  of  phase  transitions  that  have  come  about 
from  extensive  study  of  its  properties.  It  can  be  used  as  a  model  of  any 
system  involving  lattice  sites  with  two  states  per  site.  The  lattice  gas 
model  in  which  the  volume  of  the  gas  is  divided  into  cells  which  can 
either  be  empty  or  occupied  by  one  gas  molecule,  the  ferromagnetic  and 
antiferromagnetic  in  which  each  site  has  either  an  'up'  spin  or  a  'down* 
spin  and  the  binary  alloy  with  either  A  or  B  type  molecules  at  the  sites 
are  all  examples  to  which  the  Ising  model  can  readily  be  applied.  In 
terms  of  this  model,  then,  we  can  understand  the  similarity  between  the 
different  physical  systems  outlined  in  section  A. 

As  an  example  of  the  predictions  of  the  model,  the  index  *3'  de¬ 
scribing  the  behavior  of  the  coexistence  and  magnetization  curves  in  the 
region  of  the  critical  point  is  predicted  by  three  dimensional  Ising 


1 
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models  to  be  given  by  (Fisher  and  Essam  1963), 

(3  =  0.312  +  0.002,  (1.14) 

which  is  remarkably  close  to  the  experimental  value  of  almost  one  third 
quoted  in  section  A.  This  result  strongly  indicates  that  the  model  in¬ 
corporates  most  of  the  essential  features  determining  the  nature  of  the 
critical  region,  but  the  small  discrepancy  with  experiment  is  significant. 
One  of  the  outstanding  theoretical  questions  still  remaining  is  exactly 
what  has  been  omitted  in  setting  up  the  Ising  model  Hamiltonian. 

One  could  go  on  and  describe  many  of  the  other  interesting  results 
given  by  the  model.  However,  we  shall  summarize  by  noting  that  in  the 
past  the  vast  amount  of  work  done  was  concerned  with  equilibrium  proper¬ 
ties  such  as  specific  heat,  susceptibility  or  compressibility  and  nature 
of  coexistence  curves.  Recently  there  has  been  a  shift  in  emphasis  to  a 
study  of  the  non-equilibrium  or  dynamic  properties  of  the  model.  Dis¬ 
cussions  of  critical  point  scattering  by  Fisher  and  Burford  (1967)  and 
isothermal  frequency  dependent  susceptibility  by  Glauber  (1963)  are  exam¬ 
ples. 

The  second  part  of  this  thesis  is  an  example  of  such  a  study.  We 
determine  exactly  for  two-dimensional  Ising  models  the  frequency  de¬ 
pendent  perpendicular  susceptibility  (in  the  context  of  magnetism  a  meas¬ 
ure  of  the  response  of  the  system  to  a  magnetic  field  applied  perpendicu¬ 
lar  to  the  direction  along  which  the  spins  can  align)  and  also  the  in- 
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elastic  neutron  scattering  cross-section.  Just  as  the  equilibrium  prop¬ 
erties  of  the  model  have  so  closely  mirrored  the  properties  of  real  sys¬ 
tems,  we  hope  that  such  studies  of  the  dynamical  features  of  the  model 
will  lead  to  further  understanding  of  critical  phenomena. 

Thus  far,  in  our  discussion  of  models,  we  have  concentrated  on 
the  Ising  model  since  it  is  the  ’prototype1  of  all  such  models.  However, 
particularly  with  regard  to  ferromagnetism,  the  Heisenberg  model  is  at 
least  equally  important.  Its  physical  origin  lies  in  the  quantum  me¬ 
chanical  exchange  forces  that  exist  between  systems  of  electrons, and  its 
discovery  by  Heisenberg  (1928)  provided  physicists  with  a  much  sought- 
after  source  of  ferromagnetism.  A  simple  derivation  of  the  model  fol¬ 
lowing  closely  Heisenberg’s  original  arguments  will  be  given  in  the  next 
chapter.  This  will  lead  naturally  to  the  introduction  of  the  spin  one 
exchange  model  of  ferromagnetism  the  study  of  which  forms  the  first 
part  to  this  thesis.  The  introduction  of  this  model  is  an  attempt  to 
extend  the  very  limited  number  of  models  which  have  been  investigated. 

One  of  the  principal  aims  is  to  shed  some  more  light  on  the  question 
that  has  been  mentioned  several  times  throughout  this  introduction  i.e. 
how  important  is  the  nature  of  the  interaction  in  determining  the  char¬ 
acteristics  of  the  critical  region. 

Chapter  three  contains  a  number  of  preliminaries  to  the  main  cal¬ 
culation,  namely  an  introduction  to  high  temperature  series  expansions 
and  some  discussion  of  graph  theory.  Detailed  consideration  is  given  in 
chapter  four  to  the  two  main  techniques  of  evaluating  exact  series  expan- 
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sions,  the  direct  method  and  the  linked  cluster  method, and  the  applica¬ 
tion  of  the  latter  to  the  spin  one  exchange  model  is  described.  Chapter 
five  discusses  the  group  theoretic  techniques  employed  to  simplify  the 
calculation  and  part  one  of  the  thesis  concludes  with  detailed  analysis 
and  discussion  of  the  results  in  chapter  six. 

Part  two  begins  with  an  introduction  to  the  concept  of  the  per¬ 
pendicular  susceptibility  of  the  Ising  model  and  how  it  can  yield  dynamic 
features  of  the  model.  Also  in  chapter  seven,  the  Kubo  formalism  on  which 
the  whole  calculation  is  based  is  discussed  in  general  terms  and  then 
applied  to  the  problem  of  the  perpendicular  susceptibility.  The  next 
chapter  contains  some  discussion  of  the  interesting  exact  results  obtain¬ 
ed  for  the  plane  square  and  honeycomb  lattices.  The  resolution  of  some 
apparent  discrepancies  is  then  described  in  chapter  nine.  In  chapter 
ten,  the  same  techniques  used  in  this  calculation  are  applied  to  deter¬ 
mine  the  inelastic  neutron  scattering  cross-section  of  the  Ising  model. 
Finally  in  chapter  eleven  a  general  survey  of  the  results  obtained  in 
both  parts  of  the  thesis  is  given. 
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PART  I 


THE  SPIN  ONE  EXCHANGE 


MODEL  OF  FERROMAGNETISM 
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CHAPTER  II 

THE  EXCHANGE  INTERACTION 


We  shall  consider  a  quantum  mechanical  system  of  two  identical 
particles  (labelled  1  and  2)  of  spin  S  with  an  interaction  potential 
Vi £  which  is  symmetric  under  exchange  of  the  two  particles.  Our  aim  is 
to  develop  an  expression  for  what  is  known  as  the  exchange  interaction 
particularly  for  the  cases  S  =  \  and  1, though  all  our  arguments  can 
easily  be  extended  to  higher  spin. 

We  introduce  the  operator  P-^  which  represents  the  exchange  of 
the  two  particles.  It  has  two  eigenvalues,  +1  for  symmetric  states 
and  -1  for  antisymmetric  states.  Since  both  total  spin  Jl2  =  ~1  +  -2 
and  total  z  component  of  spin  S-^>  =  S^  +  s|  commute  with  P^£»  we  can 
conclude  that  any  simultaneous  eigenstate  of  and  S-^  has  a  definite 

symmetry  (i.e.  is  an  eigenstate  of  P^).  Indeed  any  non-degenerate 
eigenfunction  of  a  commuting  set  of  symmetric  operators  will  likewise 
be  either  symmetric  or  antisymmetric.  In  addition,  since  P12  commutes 
with  *  (S*  +  S*)  ±  i(S^  +  Sj?),  the  symmetry  character  of  the  si¬ 
multaneous  eigenstate  is  independent  of  S and  is  determined  solely  by 

the  value  of  ^3.2*  This  means  that  any  eigenstate  of  even  if  it  is 

2 

not  also  an  eigenstate  of  S12  will  have  a  definite  symmetry  since  it 
can  be  considered  as  a  linear  combination  of  the  simultaneous  eigen¬ 
states  of  jSjj,  and.  . 


Such  a  linear  combination  will  have  the  same 


\ 


2.2 


symmetry  character  as  each  of  its  components. 

One  can  also  show  that, as  the  total  spin  quantum  number  S-^ 
decreases  from  2S  to  0,  the  corresponding  sets  of  eigenstates  alternate 
in  symmetry  character  with  the  * 2S  *  states  always  being  symmetric,  A 
simple  proof  due  to  J,  Hamilton  is  contained  in  a  paper  on  the  general 
connection  between  exchange  and  spin  by  Schrodinger  (1941),  We  shall 
illustrate  it  for  the  case  of  spin  1  particles. 

We  try  by  inspection  to  write  down  systematically  all  the  eigen¬ 
states  of  of  definite  symmetry.  Then  by  a  process  of  elimination, 
we  determine  which  eigenvalue  of  S-^>  belongs  to  that  state  and  hence 
the  symmetry  character  associated  with  that  eigenvalue.  Let  a(i),  0(i), 

y(i)  represent  the  eigenstates  of  belonging  to  the  eigenvalues  1, 

2 

0,  -1  respectively,  For  the  case  of  =  2,  there  is  only  one  possible 
eigenstate  of  definite  symmetry,  namely  a(l)a(2)»  This  state  must  also 
have  S-^2  =  2  and  therefore  all  eigenstates  of  with  =  2  are  sym¬ 

metric,  For  s=  1,  we  can  have  either  ot(l)0(2)  +  a(2)0(l)  or 
a(l)0(2)  -  a(2)0(l).  The  former  must  be  associated  with  =  2  since 
as  we  have  said,  all  such  states  are  symmetric  and  hence  the  latter  is 

a  state  whose  total  spin  quantum  number  is  1.  We  conclude  that  all 

2 

states  with  =  1  are  antisymmetric.  When  Sj2~  there  exist  three 
possibilities  viz.  a(l)y(2)  +  a(2)y(l),  a(l)y(2)  -  a(2)y(l)  and  0(l)0(2) 


i.e.  two  symmetric  and  one  antisymmetric  states.  Now  one  symmetric 
state  must  be  associated  with  S-^  =  2,  the  antisymmetric  state  with 
S-j 2  =  1  and  hence  the  state  with  =  0  must  be  symmetric.  Thus  as  we 
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stated  earlier,  the  symmetry  character  alternates  from  S-^>  =  2  to 
S12  =  0,  For  the  spin  \  case,  we  would  similarly  find  that  S^>  =  1 
states  are  symmetric  (the  triplet  states)  and  for  =  0,  the  state  is 
antisymmetric  (the  singlet  state).  This  conclusion  that  the  symmetry 
character  and  total  spin  quantum  number  of  an  eigenstate  are  closely 
linked  forms  an  important  step  in  the  argument  that  Heisenberg  (1928) 
used  to  derive  an  expression  for  the  exchange  operator  P12*  Our  deri¬ 
vation  will  closely  follow  an  account  given  by  Van  Vleck  (1932). 

Let  us  first  suppose  that  the  two  particles  do  not  interact  and 
for  the  moment  we  do  not  consider  spin.  The  Schrodinger  wave  equation 
for  the  system  then  is 


+  YoY(2) 


V(l)  -  V(2)j  tp  =  0 


(2.1) 


where  m  is  the  mass  of  each  particle  and  V(i)  is  the  potential  to 
which  the  ith  particle  is  subjected  independently  of  the  presence  of  the 
other  particle.  The  solutions  of  this  equation  are  of  course  of  the 
form  f 

^(1,2)  =  \(DV2)  i  (2.2) 


associated  with  the  energy  eigenvalue  W  =  +  Wm  where  ^  and  <J.'m  are 

solutions  of  the  single  particle  Schrodinger  equation  belonging  to  the 
eigenvalues  and  Wm  respectively.  An  alternative  solution,  however, 
exists  associated  with  the  same  eigenvalue,  namely, 


Vl2)  =  V2)VD  , 


(2.3) 


where  (2.3)  has  been  obtained  from  (2.2)  by  exchanging  the  two  particle 
This  exchange  degeneracy  is  removed  when  a  symmetric  interaction  be¬ 
tween  two  particles  V-j^  is  introduced.  It  is  a  simple  exercise  to  show 
that  the  two  eigenstates  now  become 


^  sym  “  ^1  + 


/2 


(2.4) 


and 


^ant  =  ~  (^1  ~  ^ij) 


(2.5) 


associated  respectively  with  the  eigenvalues. 


Wo  +  K12  +  J12  1 


W  = 


(2.6) 


WQ  +  Ki2  ~ 


(b) 


where  WQ  is  the  energy  in  absence  of  the  -  interaction  term  V^> » 


k32  = 


hV12Vvldv2  =  'hlV12IJ'lIdvldv2 


(2.7) 


and 


J12  = 


t;IV12t,'lIdvldv2  . 


(2.8) 
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We  now  consider  the  inclusion  of  spin.  For  spin  j  particles,  the 
Pauli  exclusion  principle  demands  that  the  total  wave  function  be  anti¬ 
symmetric  under  exchange  of  the  two  particles.  In  the  absence  of  any 
coupling  between  space  variables  and  spin  variables,  the  total  wave 
function  will  be  the  product  of  a  space  dependent  part  and  a  spin  de¬ 
pendent  part.  Then  the  spin  wave  function  associated  with  (2.4)  must  be 
antisymmetric  and  the  one  associated  with  (2.5)  must  be  symmetric  i.e. 
the  eigenvalues  of  P-^  are  -1  and  +1  respectively  where  now  P^  repre¬ 
sents  the  exchange  of  only  the  spins  of  the  particles.  In  other  words, 
has  eigenvalues  +1  and  -1  for  states  where,  from  (2.6),  V-^,  has 
eigenvalues  anc*  ^12  +  ^12®  can  therefore  set  up  the  im¬ 

portant  operator  identity. 


V12  =  K12  “  J12P12  • 


(2.9) 


Heisenberg's  principal  achievement  was  essentially  to  demonstrate  (2.9) 
for  the  case  of  electrons  where  is  the  familiar  coulomb  interaction. 
We  should  emphasize  that  this  exchange  operator  is  simply  a  manifesta¬ 
tion  in  spin  space  of  the  inter-electron  coulomb  repulsion  which  comes 
about  principally  because  of  the  Pauli  exclusion  principle.  The  often 
used  phrase  that  the  exchange  interaction  has  no  classical  analogue  is 
therefore  a  little  misleading. 

The  spin  exchange  operator  P-jj?  cun  be  put  in  a  much  more  familiar 
form.  For  the  antisymmetric  spin  state  of  two  spin  -§  particles,  we 


2.6 


found  that  S-^  must  equal  0.  As  a  consequence  the  operator  S-^.Sg  given 


= 


1  (s2 
2  ~12 


A  -  &> 


(2.10) 


1  (S  2 

2  '-12 


2) 
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(2.11) 


must  have  an  eigenvalue  of  -  f-  in  that  state.  Similarly  for  the  sym- 

Hr 

metric  state,  S-^.Sg  has  an  eigenvalue  of  2.  It  is  therefore  possible  to 
set  up  another  important  operator  identity  of  the  form  , 


P!2  ”  a  +  b§]_»S2  , 


(2.12) 


where  a  and  b  are  constants  determined  by  the  simultaneous  equations  , 


The  solution  yields  the  familiar  Dirac  relation  , 


(2.13) 


^  +  2Si'§2 


(2.14) 


It  is  in  this  form  that  the  Heisenberg  exchange  interaction  is  usually 
written  i.e.  neglecting  a  relatively  unimportant  constant, 
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V12  =  -%*S2  a  (2.15) 

The  purpose  of  this  detailed  discussion  has  been  not  only  to  emphasize 
the  physical  origins  of  exchange  but  to  introduce  the  spin  one  exchange 
operator  which  is  obtained  by  a  natural  extension  of  equations  (2,9)  to 
(2,15)  to  the  spin  one  case.  For  the  symmetric  spin  state,  S-^>  must  be 
2  or  0  and  1  for  the  antisymmetric  state.  Proposing  for  P^2  a  form, 

P12  =  a(S1#S2)2  +  b(Sl9S2)  +  c  ,  (2.16) 

we  obtain  the  operator  identity  , 

P12  “  +  (~la~2^  “  1  *  (2.17) 

and  hence  a  spin  one  exchange  interaction  which  can  be  written  in  the 
form  , 

V-^  =  ~JL(S19S2)2  +  (Sl0S2)-l]  .  (2.18) 

A  model  where  the  basic  interaction  is  given  by  (2.18)  will  be  investi¬ 
gated  in  the  following  chapters. 

The  Heisenberg  exchange  given  by  (2.15)  has  also  been  investi¬ 
gated  for  higher  spin  systems.  In  these  cases,  it  represents  the  sum 
of  all  the  single  electron  exchanges  between  the  electrons  in  the  spin 
system  with  the  assumption  that  the  interaction  constant,  J,  is  the 
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same  for  all  such  exchanges.  For  a  spin  1  system  (2.18)  represents  the 
simultaneous  exchange  of  the  system  of  electrons  in  a  relative  spin  1 
state.  In  other  words,  if  the  electrons  are  so  coupled  that  they  can 
be  considered  as  a  single  spin  1  system  rather  than  as  a  system  of  in¬ 
dependent  electrons,  we  might  expect  (2.18)  to  be  significant. 

There  are  other  important  reasons  for  considering  a  spin  1  ex¬ 
change  model.  For  spin  1  systems,  the  most  general  form  of  isotropic 
spin-spin  interaction  is  given  by  (2.16),  Since  S-^.Sg  has  only  three 
distinct  eigenvalues,  it  is  easy  to  show  that  (S^.Sg)11  for  n  >  2  can  be 
expressed  in  the  form  (2.16).  While  such  interactions  have  not  been 
considered  previously  by  theoreticians  in  this  field,  it  should  be  in¬ 
teresting  to  see  what  the  most  general  form  of  interaction  for  a  spin  1 
system  predicts  particularly  in  comparison  to  the  spin  1  Heisenberg  model. 
The  reason  that  the  particular  form  (2.17)  is  chosen  is  partly  for  argu¬ 
ments  given  in  the  previous  paragraph  and  partly  since  as  we  shall  find 
later,  the  spin  1  exchange  operator  allows  for  considerable  mathematical 
simplification.  As  we  have  said  earlier,  we  can  also  hope  to  shed  some 
light  on  the  question  of  how  important  is  the  nature  of  the  interaction 
in  determining  the  details  of  a  phase  transition. 

Finally,  let  us  mention  an  important  case  where  an  interaction  of 
the  form  (2.16)  arises.  In  spin  notation,  the  components  Qi^(i, j = x,y,z) 
of  the  quadrupole  moment  tensor  are  given  by, 


+S‘^Si 
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(2.19) 
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Segre  (1953)  has  shown  that  a  quadrupole-quadrupole  interaction  of  the 
form  9 

=  J  2  ,  (2.20) 

ij 


can  be  written  as, 


JL9(S1.S2)  4-  ^  (Sp*^) 


2  2 


] 


(2.21) 


Now  Nakamura  (1955)  has  argued  that  the  principal  interaction  in  solid 
hydrogen  with  a  high  concentration  of  ortho-molecules  (>  60$)  is  an 
electric  quadrupole-quadrupole  interaction  where  the  interaction  con¬ 
stant,  J,  is  given  explicitly  by, 


6e2o£ 

j  = - ±  .  (2.22) 

25  R3 

»e’  is  the  electronic  charge,  Qjj  the  quadrupole  moment  of  the  ortho- 
hydrogen  molecule  and  R  is  the  distance  between  interacting  quadrupoles. 
Using  data  given  by  Nakamura,  we  find 

J  =  0.654k°K  (2.23) 

where  k  is  Boltzmann’s  constant.  That  the  interaction  constant  can  be 
evaluated  is  in  direct  contrast  with  most  magnetic  systems.  Another 
feature  worth  noting  is  that  the  dependence  of  J  would  make  the 

r3 
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normal  nearest  neighbour  interaction  approximation  a  good  one  in  this 
instance. 

While,  of  course,  (2.21)  does  not  have  its  physical  origins  in 
exchange,  it  is  nevertheless  formally  an  example  of  (2.16).  Unfortu¬ 
nately  for  statistical-mechanical  calculations,  the  interaction  constant, 
as  we  have  seen  from  (2.22),  is  positive  and  leads  to  an  ’antiferromag¬ 
netic1  type  of  interaction.  The  case  we  shall  investigate  in  the  fol¬ 
lowing  chapters  is  the  one  where  the  ferromagnetic  state  is  preferred. 
Experience  has  shown  that  series  expansion  techniques  give  more  accurate 
results  for  ferromagnetic  systems. 
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CHAPTER  III 


AN  INTRODUCTION  TO  HIGH  TEMPERATURE  SERIES  EXPANSIONS 


A,  Review  of  Statistical  Mechanics 

As  we  have  just  indicated  in  the  previous  chapter,  the  spin  one 
exchange  model  of  ferromagnetism  is  characterized  by  the  Hamiltonian, 


(3.1) 


X  =  X0  + 


(3.2) 


where  m  is  the  magnetic  moment  associated  with  each  spin  and  H  represents 
an  external  magnetic  field  aligned  along  the  z  direction  of  our  coordi¬ 
nate  system.  In  71q,  i  and  j  label  the  sites  of  a  regular  lattice  (in 


our  case  we  shall  consider  the  face-centred  cubic  lattice)  and  the  sum¬ 


mation  is  taken  over  all  nearest  neighbour  pairs  of  sites.  The  inter¬ 
action  constant  J  is  considered  positive  so  that  the  system  is  ferromag¬ 
netic.  The  interaction  differs  by  a  constant  from  that  introduced  in 
(2.18).  This  alteration  has  no  effect  on  the  final  results  but  it  does 
introduce  significant  simplification  into  the  calculation  by  making  the 
interaction  matrix  traceless.  74'-^  represents  the  interaction  of  the  sys¬ 


tem  with  an  external  magnetic  field.  We  shall  always  consider  the 


/ 
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limiting  case  of  this  field  tending  to  zero,, 

Our  aim  is  to  proceed  from  the  microscopic  description  of  our  sys¬ 
tem  as  given  by  (3*2)  to  calculate  the  familiar  macroscopic  thermodynamic 
variables  such  as  specific  heat,  entropy  and  magnetic  susceptibility. 

The  whole  formalism  of  statistical  mechanics  is  set  up  primarily  for 
this  purpose.  For  systems  with  a  fixed  number  of  particles  N,  in  thermal 
contact  with  a  heat  bath,  the  function  which  provides  this  important 
bridge  between  the  microscopic  and  the  macroscopic  is  the  canonical  par¬ 
tition  function  given  for  quantum  mechanical  systems  by, 


Z(T,H)  =  trace (exp (-3?$  2  * 


(3.3) 


where  0  =  l/kT  (k  is  Boltzmann’s  constant)  and  U  is  the  Hamiltonian  ma¬ 
trix  operator  given  for  our  model  by  (3.2).  The  fundamental  equation  in 
statistical  mechanics  links  the  partition  function  to  the  Helmholtz  free 
energy  F,  namely, 


F  =  -kT  log  Z  . 


(3.4) 


Simple  thermodynamics  then  gives  for  the  entropy  S, 


S  =  JL  (kT  log  Z)  , 


(3.5) 


for  the  magnetic  specific  heat  C, 
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C  ^  M  (log  Z)A 

and  for  the  magnetic  susceptibility  X  » 

.2 

%  =  (kT  log  Z)  . 
dH 


(3.6) 


(3.7) 


It  should  be  noted  that  in  (3.5-7)  no  indication  is  given  of  whether  the 
derivatives  are  taken  at  constant  magnetic  field  or  constant  magnetiza¬ 
tion.  We  shall  always  be  considering  the  H  =  0  limit  which  for  a  ferro¬ 
magnetic  in  the  high  temperature  region  above  the  critical  point  is 
identical  with  the  M=0  condition. 

As  mentioned  in  chapter  I,  in  the  high  temperature  series  expan¬ 
sion  method,  (3.3)  is  expanded  as  a  power  series  in  inverse  temperature. 
If  we  introduce  the  convenient  notations, 


P=2  [(Si-S  /  +  (S^S.)  -  i]. 


nn 


3 


Q  =  XI  Sf  , 
i 


(3.8) 


a  = 


mH 

J  ’ 


substitute  (3.2)  in  (3*3)  and  expand  the  exponential,  we  obtain 


■ 
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Z(T,H) 


4r-i  trace(P  i-aQ)r  T,„ 
-2_i  TT  “  K 

r  =  0  rI 


(3.9) 


oo 


^ — j  trace(P+aQ) 
trace  I  +  2_j  - rj 

r  =  1 


Kr 


(3.10) 


I  is  the  unit  matrix  which  for  a  system  of  N  spin  1  particles  has  dimen< 
sions  3W  x  3N.  (3.10)  can  then  be  written  in  the  form, 


o°  ,, 

Z(T,H)  =  3w(l  +  2]  yy  Kr) 

r  =  1 


(3.11) 


where 


^r 


trace  (P  +  gQ)r 


(3.12) 


To  obtain  the  entropy  etc,,  we  are  primarily  interested  in  log  Z(T,H) 
which  from  (3*11)  is  given  by, 

00  n 

log  Z(T,H)  =  N  log  3  +  log[l  +  22  -f  Kr]  (3.13) 

r  =  1  1 


°°  K 

=  m  log  3  +  rr  Kr  .  (3.14) 

r  =  1  1 


In  the  language  of  statistics,  the  coefficients  A.r  are  the  cumulants 
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corresponding  to  the  moments  p  ,  Explicit  relations  giving  A.^  as  com¬ 
binations  of  all  the  u  with  s  ^  r  are  to  be  found  in  appendix  A  for 

r  s 

r  ^  9  and  for  the  particular  case  of  p-^  =  0„  Alternatively  the  A.^  can 
be  regarded  as  those  parts  of  the  p  that  are  directly  proportional  to 
N.  Since  log  Z  is  an  extensive  quantity,  Kr  must  be  proportional  to  N 
and,  in  fact, 


Hr 


•<*  terms  proportional  to 


etc, , 


(3.15) 


these  latter  terms  cancelling  in  the  transformation  from  (3*13)  to  (3.H) • 
For  the  initial  magnetic  susceptibility,  Xq,  using  (3® 7) »  (3.8) 
and  (3.1H),  we  obtain 


S 

r  =  1 


1 l 


rl 


(3.15) 


where 


d2Ar(ot) 

A.1' (a)  =  lim  . -  u 

r  a— ^  0  8a2 


(3.16) 


From  (3.8),  (3.12)  and  appendix  A,  we  find 


^1  ~  Hi  -  0  ~  All 


(3.17) 


and,  as  a  result,  (3.15)  is  frequently  written  in  terms  of  the  reduced 


initial  sus 


(3.18) 


ll  ^  I  I —  t  \ 

Again  A.r  can  either  be  regarded  as  that  part  of  ]jpr  is  defined  in  an 
analogous  manner  to  A.""]  which  is  proportional  to  N  or  as  a  particular 


combination  of  all  the  p.g  and  pg  for  s  <:  r.  These  combinations  are  also 
given  in  appendix  A  for  r  q. 

The  object  of  this  technique  is  to  calculate  the  exact  numerical 
value  of  as  many  of  the  Ar  and  A.1^  as  possible.  The  former  will  give  the 
series  for  the  entropy,  specific  heat  and  internal  energy  in  zero  external 
field  (if  we  put  ct=0  in  (3.12))  and  the  latter  the  initial  magnetic  sus¬ 
ceptibility.  The  first  term  in  each  case  reproduces  the  behavior  in  the 
infinite  temperature  limit.  As  more  terms  are  evaluated,  the  series  will 
represent  a  good  approximation  to  the  exact  behavior  at  lower  and  lower 
temperatures.  Ideally  to  obtain  information  close  to  the  critical  tem¬ 
perature,  we  should  simply  calculate  a  sufficient  number  of  terms  of  the 
series.  However,  the  great  labour  involved  in  finding  these  coefficients 
strictly  limits  the  number  which  it  is  possible  to  calculate.  Fortu¬ 
nately  there  exist  techniques  which  enable  us  to  predict  the  behavior  of 
the  series  near  the  critical  temperature  when  the  truncated  series  by  it¬ 
self  would  represent  a  bad  approximation. 


First,  however,  we  shall  describe  the  two  important  and  distinct 
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methods  that  are  used  to  calculate  the  coefficients.  They  are  known  as 
the  direct  method  and  the  cluster  expansion  method.  In  both  cases, 
an  integral  part  of  the  procedure  is  to  set  up  a  correspondence  between 
the  many  traces  of  products  of  spin  variables  that  contribute  to  any 
given  A.r,  p.r  etc.  and  various  types  of  linear  graphs.  As  a  necessary 
preliminary,  we  shall  give  precise  definitions  of  many  of  the  terms 
that  now  form  part  of  the  established  'vernacular'  of  graph  theory.  We 
shall  follow  very  closely  those  given  by  Sykes  et  al,  (1966)  who  in  turn 
used  Ore  (1962)  as  their  principal  reference. 

B,  Concepts  in  Graph  Theory 

A  linear  graph  consists  of  a  set  of  vertices  (points)  V,  and  asso¬ 
ciated  pairings  of  these  vertices.  These  pairings  are  represented  by 
the  lines  (edges  or  bonds )  of  the  graph.  An  undirected  graph  is  one  in 
which  the  order  of  the  vertices  in  a  pair  is  not  considered.  We  shall  be 
concerned  throughout  solely  with  undirected  graphs. 

A  graph  is  connected  if  for  every  pair  of  vertices  there  is  at 
least  one  sequence  of  lines  called  a  path  such  that  the  first  line  is 
incident  upon  one  of  the  vertices  and  the  last  line  incident  on  the 
other,  consecutive  lines  always  having  a  common  vertex  and  no  line  being 
repeated  in  the  sequence.  A  disconnected  graph  is  one  in  which  there 
exists  at  least  one  pair  of  vertices  between  which  there  is  no  such  path. 
The  graph  drawn  below  is  an  example  of  a  disconnected  undirected  linear 
graph  of  three  connected  components  for  which  six  vertices  constitute 
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the  vertex  set  and  four  edges  the  edge  setc 


The  degree  (or  order)  of  a  vertex  is  equal  to  the  number  of  lines 
which  are  incident  on  that  vertex.  The  above  graph  has  three  vertices 
of  order  two,  two  of  order  one  and  one  of  zero  order.  Vertices  of  order 
greater  than  two  are  known  as  nodes  while  those  of  degree  two  or  less 
are  called  antinodes.  A  bridge  is  a  simple  path  connecting  two  nodes,  a 
simple  path  being  one  which  traverses  a  given  vertex  no  more  than  once. 
The  suppression  of  a  vertex  of  degree  two  is  defined  as  the  removal  of 
the  vertex  and  the  identification  of  the  two  lines  incident  upon  it. 
Conversely  the  insertion  of  a  vertex  of  degree  two  is  the  operation  of 
replacing  a  line  incident  upon  two  such  vertices  by  two  lines  and  one 
vertex,  each  line  incident  upon  one  of  the  original  vertices  and  the  two 
lines  together  incident  upon  the  new  vertex.  The  two  operations  are  il¬ 
lustrated  below  with  the  vertex  in  question  marked  by  an  asterisk. 


z  INSERTION 
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Two  graphs  are  isomorphic  if  there  is  a  one  to  one  correspondence 
between  their  vertices  such  that  two  vertices  are  joined  by  a  line  in 
one  of  them  if  and  only  if  they  are  joined  by  a  line  in  the  other.  The 
two  graphs  drawn  below  are  isomorphic  and  the  vertices  are  labelled  to 
show  the  one  to  one  correspondence. 


Two  graphs  are  said  to  be  homeomorphic  if  one  of  them  may  be 
derived  from  the  other  by  the  insertion  or  suppression  of  any  number  of 
vertices  of  degree  two.  Such  graphs  are  said  to  be  of  the  same  topo¬ 
logical  type  or  have  the  same  topology.  Each  graph  is  a  realisation  of 
the  topology.  For  example,  (a)  is  the  schematic  representation  of  a 
certain  topology  and  (b)  and  (c)  are  two  of  the  infinite  number  of  dis¬ 
tinct  realisations  of  that  topology. 


(a) 


(b) 


(c) 
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If  we  denote  by  n,  the  number  of  connected  components  of  a  graph, 
e  the  number  of  lines  and  v»  the  number  of  vertices,  then  we  define  the 
cyclomatic  index  c,  of  the  graph  to  be 

c  =  e  -  v  -h  n  .  (3.19) 

This  is  a  very  important  concept  in  the  classification  of  graphs.  It 
should  be  noted  that  all  homeomorphic  graphs  have  the  same  cyclomatic 
index.  For  example  all  connected  linear  chains  have  cyclomatic  index  zero. 

If  g  and  g  are  two  graphs  with  disjoint  vertex  sets  and  disjoint 
edge  sets,  we  denote  by  g^  U  gg  the  graph  whose  vertex  and  edge  sets  are 
the  unions  of  the  respective  sets  of  g  and  g  .  For  example, 

x  0 
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The  graph  H  is  a  subgraph  of  a  graph  G  when  the  vertex  set  V(H) 
is  contained  in  the  vertex  set  V(G)  and  all  the  edges  of  H  are  edges  of 
G.  A  graph  H  is  a  section  graph  of  G  if  the  vertex  set  V(H)  is  contained 
in  the  vertex  set  V(G)  and  if  the  lines  of  H  are  all  the  lines  of  G 
which  connect  the  vertices  of  H.  For  example  (b)  below  is  both  a  sub¬ 
graph  and  a  section  graph  of  (a)  but  (c)  is  only  a  subgraph. 


, 
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Any  subgraph  of  G  that  is  isomorphic  with  a  graph  g  is  said  to 

represent  a  weak  embedding  of  g  in  G.  Likewise  any  section  graph  of  G 
which  is  isomorphic  with  g  is  said  to  represent  a  strong  embedding  of  g 
in  G.  Clearly  a  strong  embedding  is  also  a  weak  embedding  but  not  vice 
versa.  The  lattice  constant  of  a  graph  g  is  defined  in  two  senses.  The 
weak  lattice  constant  of  g  on  a  graph  G,  denoted  by  (g;G),  is  equal  to 
the  number  of  subgraphs  of  G  which  are  isomorphic  with  g.  Correspond¬ 
ingly,  the  strong  lattice  constant  of  g  on  G,  denoted  by  |_g ;  G J  is  the 
number  of  section  graphs  of  G  isomorphic  with  g.  For  example 


=  0 


=  2 


It  turns  out  that  the  strong  lattice  constants  arise  in  low  tern- 
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perature  calculations.  In  our  case,  we  shall  be  concerned  only  with  the 
weak  lattice  constants  and  in  succeeding  chapters  this  will  always  be 
implied  even  if  specific  reference  to  'weak'  is  not  made. 

It  is  necessary  to  generalise  the  concept  of  lattice  constant  to 
the  case  when  G  is  an  infinite  lattice.  Clearly  as  defined  the  lattice 
constant  for  this  case  would  be  either  zero  or  infinity.  If  we  consider 
first  a  lattice  of  N  sites  and  ignore  boundary  effects,  the  lattice  con¬ 
stant  can  be  expressed  as  a  function  of  N.  Connected  graphs  have  lat¬ 
tice  constants  proportional  to  N  and  since  for  the  most  part,  we  shall  be 
concerned  with  connected  graphs  we  redefine  lattice  constant  to  mean  the 
number  of  embeddings  per  site  when  the  term  is  applied  to  an  infinite 
lattice. 

To  conclude  these  introductory  remarks  on  some  of  the  relevant 
concepts  of  graph  theory,  we  emphasize  that  thus  far  it  has  been  implic¬ 
itly  assumed  that  a  pair  of  vertices  could  be  connected  by  no  more  than 
one  edge.  The  graph  concept  can  be  extended  to  allow  a  pair  of  vertices 
to  be  joined  by  several  distinct  edges.  This  is  particularly  important 
to  our  discussion  of  the  ‘direct’  method  of  series  analysis.  This  to¬ 
gether  with  our  calculation  using  the  'cluster  expansion'  method  will  be 
discussed  in  detail  in  the  next  chapter. 
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CHAPTER  IV 

THE  EVALUATION  OF  SERIES  COEFFICIENTS 


A,  The  Direct  Method 

It  is  clear  from  (3.8)  and  (3.10)  that  in  order  to  calculate  the 
high  temperature  series  for  the  zero  field  partition  function,  we  have 
to  evaluate 


trace(Pn)  =  trace (P^  .  -  i)]  n  •  (4.1) 

n.n.  J 

It  is  very  convenient  to  set  up  a  correspondence  between  the  products  of 
spin  operators  that  contribute  to  (4.1)  for  a  given  'n1  and  sets  of 
linear  graphs  of  n  bonds.  The  graphs  are  unlabelled  and  the  particular 
sites  to  which  the  spin  operators  are  attached  is  of  no  relevance.  Only 
the  topology  of  the  spin  products  is  important.  Thus  for  n=2,  both 
(Pf2  -  1/3)  and  (P23  -  l/3)^  are  represented  by  the  graph 
and  (Pp2  "  l/ 3 ) ( ^23  “  l/3)  will  be  represented  by  the  graph 
whatever  the  angle  between  the  bonds  (12)  and  (23)  on  the  lattice  under 
consideration.  There  are  then  three  main  steps  in  the  calculation,  namely 

a)  for  a  given  1  n',  to  enumerate  all  possible  linear  graphs  con¬ 
sisting  of  n  bonds, 

b)  to  calculate  the  lattice  constants  of  these  graphs  and 
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c)  to  evaluate  the  traces  of  all  the  products  of  spin  variables 
corresponding  to  each  graph. 

The  enumeration  of  graphs  has  been  considerably  simplified  by  our 
modification  of  the  nearest  neighbour  interaction  from  P .  .  to  P.  .  -  l/3. 
In  the  latter  form,  the  following  rearrangement  is  possible,  namely, 


m=  1 


m  m 
a .  a  . 
1  J 


* 


(4.2) 


where  the  ara(m  =  1,2... 12)  are  all  traceless  matrices.  (They  are  listed 
in  appendix  B.)  The  Heisenberg  interaction  is  of  this  form  where  'a1  is 
the  spin  operator  S  and  m  =  x,  y  and  z. 

The  spin  space  of  a  collection  of  spins  is  the  direct  product  of 
the  individual  spin  spaces.  Any  product  of  spin  operators  belonging  to 
different  lattice  sites  is  actually  a  direct  product  of  the  correspond¬ 
ing  spin  matrices.  For  example,  any  term  from  (4.2)  when  written  in  a 
matrix  representation  takes  the  form, 

a™  a”  =  (a”  x  I). (I  x  a”)  (4.3) 

=  (a™  .  I)  x  (I.a“)  (4.4) 

m  m 

=  a-^  x  a^ 


(4.5) 


where  on  the  right  hand  side  of  (4.3-5)  the  1 .'  refers  to  ordinary 
matrix  multiplication  and  the  1  x’  to  direct  product  multiplication. 

From  a  well  known  property  of  the  trace  we  can  immediately  conclude  that 


m 


trace  (at'1  x  a™)  =  (trace  a™)(traee  a“) . 


rrn 


(4.6) 


(4.3-6)  can  easily  be  extended  to  the  general  case  of  a  product  of  any 
number  of  spin  variables  operating  on  different  lattice  sites.  Since 
the  a*s  have  all  be  constructed  to  be  traceless,  these  results  have  the 
important  consequence  that  any  product  involving  only  one  1  am’  on  at 
least  one  lattice  site  will  automatically  be  traceless  *  Such  a  situa¬ 
tion  arises  in  all  spin  products  corresponding  to  graphs  with  at  least 
one  vertex  of  order  one.  As  an  example,  in  P ^  a  typical  product  cor- 
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responding  to  the  graph 


...  2 


is  given  by, 


trace ( aqa2a2 a2  a4a4a3a3 ) 


(4.7) 


=  (trace  a^) (trace  a^a^a^) (trace  a^a^)( trace  a^a^) 


=  0 


(4.8) 


since  trace  a^  =  0  . 


(4.9) 


'  ' 
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We  can  therefore  ignore  all  such  graphs  (often  referred  to  as  open  graphs 
or  configurations )  in  our  enumeration  and  this  represents  an  important 
simplification.  Nevertheless  the  number  of  graphs  becomes  quite  large  as 
the  number  of  bonds  increases.  Rushbrooke  and  Wood  (1958)  using  this 
method  on  the  Heisenberg  model  report  that,  even  with  further  simplifica¬ 
tions  they  were  able  to  introduce,  they  enumerated  a  total  of  120  graphs 
with  a  maximum  of  six  bonds  of  which  91  were  disconnected. 

In  enumerating  these  graphs  no  account  is  taken  of  the  specific 
lattice  structure.  The  features  of  a  given  lattice  arise  in  the  calcula¬ 
tion  of  the  lattice  constants.  For  instance,  for  loose-packed  lattices 
such  as  the  body-centred  cubic,  all  graphs  which  contain  polygons  with  an 
odd  number  of  sides  e.g.  triangles,  pentagons,  tetrahedra  etc.  cannot 
occur  and  their  lattice  constants  are  zero.  The  fundamental  parameters 
that  characterize  any  given  lattice  are  its  coordination  number  1 q '  and 
the  lattice  constants  of  the  closed  connected  graphs  i.e.  connected 
graphs  whose  vertices  are  all  of  degree  greater  than  or  equal  to  two. 

For  the  well  known  lattices,  many  of  these  constants  are  listed  by  Domb 
(i960)  in  his  major  review  article  on  critical  phenomena.  All  other 
lattice  constants  namely  of  disconnected  and/or  open  graphs  can  be  ex¬ 
pressed  in  terms  of  these  parameters.  The  techniques  used  are  very 
clearly  described  by  Domb  (i960)  and  he  also  lists  many  of  these  addi¬ 
tional  constants.  It  is  important  to  note  that  the  lattice  constant 
for  a  connected  graph  is  directly  proportional  to  N,  the  number  of  sites 
in  a  given  regular  lattice.  For  a  disconnected  graph,  however,  con- 
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sisting  of  n  connected  components,  the  lattice  constant  has  terms  propor¬ 
tional  to  Nn,  N  ^...N.  As  we  mentioned  earlier,  in  transforming  from  Z 
to  log  Z,  all  terms  involving  powers  of  11  higher  than  the  first  cancel 
so  that  for  most  purposes  we  are  primarily  interested  only  in  that  part 
of  the  lattice  constant  linear  in  JM. 

The  most  difficult  step  by  far  in  this  method  is  to  enumerate  and 
calculate  the  traces  of  all  the  different  products  of  spin  variables 
associated  with  the  same  graph.  The  fundamental  problem  is  simply  that 
the  permutation  operators  that  constitute  our  interaction  do  not  commute 
e.g. 

|>12,  p23J  4  o  (4.io) 

or  equivalently, 

[a™,  ai-J  ^  0  for  a 11  m  ^  n  ,  (4.11) 

Because  of  this,  for  example,  the  terms  (P-^  ~  l/3)(?23  ”  -  l/3) 

and  (P-j-2  -  l/3)(P^^  -  I/3XP23  -  l/ 3)  both  associated  with  the  graph 

have  different  traces.  The  number  of  distinct  arrangements 

of  the  interactions  all  corresponding  to  the  same  graph  that  have  dif¬ 
ferent  traces  increases  enormously  as  the  number  of  bonds  in  the  graph 
increases.  For  instance,  Rushbrooke  and  V/ood  (1958)  point  out  that 


graphs  such  as 


give  rise  to  as  many  as  3^0 
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separate  terras.  The  work  then  is  very  tedious  and  lengthy  and  it  is  a 
triumph  of  ingenuity  and  perseverance  that  Rushbrooke  and  Wood  (1958) 
obtained  by  this  method  six  terms  of  the  specific  heat  series  for  the 
general  spin  Heisenberg  model. 

We  should  note  that  this  great  difficulty  is  of  course  absent  in 
such  classical  models  as  the  Ising  model.  As  a  result,  Hunter  (1967) 
has  calculated  twelve  terms  of  the  specific  heat  series  for  the  Ising 
model  on  the  face-centred  cubic  lattice  by  the  direct  method. 

Thus  far  our  remarks  have  centred  round  the  calculation  of  the 
partition  function.  For  the  initial  susceptibility  given  by  (3.9)  and 
(3.10),  we  have  to  calculate  k” .  Now  from  (3.12),  is  given  by, 

|i"  =  -pjp  tracetP^Q2)  ,  (4.12) 

since  P  and  Q  commute.  As  with  the  calculation  of  trace  (Pn)  so  also 
with  (4.12)  we  have  a  correspondence  between  terms  in  the  expansion  of 
the  product  and  linear  graphs.  Only  now  our  graphs  are  decorated  by  the 
addition  of  two  crosses.  These  crosses  occur  on  the  sites  occupied  by 
the  terms  from  Q  .  Rushbrooke  and  Wood  (1958)  discuss  three  theorems 
which  enable  them  to  calculate  the  susceptibility  series  without  evalu¬ 
ating  any  further  traces  than  were  required  for  the  partition  function 
series.  They  were  thus  able  to  determine  seven  terms  of  the  suscepti¬ 
bility  series. 


'  .  •  -- 


B.  The  Cluster  Expansion  Method 


It  seems  that  the  work  of  Rushbrooke  and  Wood  represents  the 
maximum  amount  of  information  that  can  be  obtained  practically  by  this 
method  about  quantum  mechanical  models  where  the  interactions  do  not  com¬ 
mute.  It  was  with  this  background  that  Domb  (i960)  first  tentatively 
suggested  the  cluster  expansion  approach  (which  is  an  expansion  in  terms 
of  connected  graphs  or  clusters  only).  Together  with  the  application  of 
group  theory  it  represents  a  major  advance  in  technique  and.  enabled  Domb 
and  Wood  (1965)  to  obtain  for  the  spin  l/2  Heisenberg  model  only,  an 
additional  three  terms  in  the  susceptibility  series  for  loose  packed  lat¬ 
tices  and  two  terms  for  the  specific  heat  series.  Subsequently,  Rush¬ 
brooke  et  al.  (1966)  obtained  a  total  of  eleven  terms  in  the  suscepti¬ 
bility  series  and  nine  terms  in  the  specific  heat  series  for  loose  packed 
lattices  and  ten  terms  and  eight  terms  respectively  for  the  face-centred 
cubic  lattice. 

It  is  also  the  approach  that  we  have  used  in  our  calculation  on 
the  spin  one  exchange  model.  It  has  enabled  us  to  calculate  more  terms 
and  consequently  make  more  accurate  predictions  regarding  critical  be¬ 
havior  than  has  previously  been  possible  for  any  other  spin  one  model, 
particularly  the  spin  one  Heisenberg  model. 

At  least  two  detailed  and  rigorous  derivations  of  the  cluster  ex¬ 
pansion  method  exist  in  the  literature  (Rushbrooke  (1964)  and  Sykes  et 
al.  (I966)).  Rather  than  reproduce  one  or  both  of  these  derivations,  we 


■ 
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shall  try  to  illustrate  as  clearly  as  possible  how  the  cluster  expansion 
can  be  considered  merely  as  a  rearrangement  of  the  direct  method  expansion. 

In  the  latter,  we  have  just  found  that  the  coefficients  of  the 
zero-field,  ’log  Z’  expansion  for  a  regular  lattice  L,  are  expressed  in 
the  form, 

\  =  H  2  [ooeff.  of  N  in  (gr;L)]W(gr)  ,  (4.13) 

gr 

where  the  summation  is  over  all  graphs  of  r  edges  including  disconnected 
graphs  and  multiline  graphs  i.e.  graphs  with  more  than  one  distinct  edge 
between  at  least  one  pair  of  vertices.  For  instance,  for  r=2,  three 

graphs  may  contribute  namely, 

The  W(gr)  are  the  weights  associated  with  each  graph  which  in  this  case 
are  the  sums  of  all  the  traces  of  products  of  spin  operators  which  cor¬ 
respond  to  each  graph. 

The  first  step  in  the  rearrangement  is  to  eliminate  effectively 
any  consideration  of  the  disconnected  graphs.  To  do  this,  we  realise 
that  the  coefficient  of  N  in  the  lattice  constant  of  a  disconnected 
graph  is  the  sum  of  parts  each  of  which  can  be  associated  with  one  of 
the  connected  graphs  formed  from  an  overlap  of  the  connected  components 
of  the  disconnected  graph.  An  overlap  of  a  number  of  graphs  in  this 
context  means  the  identification  of  at  least  one  of  the  vertices  from 
each  graph.  To  illustrate,  the  only  two  overlaps  of  the  disconnected 


, 
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graph 


tice  constant  of 


number  q  is  given  by. 


are 


and 


The  lat- 


on  a  regular  lattice  L  with  coordination 


;L)  =  (B.  .  2q  +  1)  . 


(4.14) 


The  coefficient  of  N  in  (4.14)  can  be  split  in  two  parts, 


f  (q  -  1)  =  -  ( 


;L)  , 


(4.15) 


and 


b) 


Nq 

w 


=  -I( 

2 


;L)  . 


In  (4,13)  therefore,  we  need  not  explicitly  consider  /  / 

instead,  alter  W(  )  by  the  contribution  -W(  ^  ) 


(4.16) 


but 


and 


W( 


)  by  the  amount  -  l/2  W(  J  /  ).  In  other  words,  (4.13) 


can  be  rewritten  in  the  form, 


\  =  2?  (gr;L)w»  (gr)  , 

gr 


(4.17) 
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where  now  the  primed  summation  is  over  only  connected  graphs  and  the 
weights  of  the  disconnected  graphs  have  been  appropriately  distributed 
among  the  connected  graphs.  For  instance  in  (4.17),  we  would  find 


a)  W‘( 


)  =  W( 


)  -  W( 


(4.18) 


b)  W'(  )  =  W(  .  .  )-|w(  J  J  )  (4.19) 

independently  of  the  lattice  L. 

Thus  (3.14)  can  be  written  as, 


log  Z(T)  = 


oo 


Z] 

r  =  0 


v  4-  (sr;L)w  (g_)  , 

e  r! 

fer 


(4.20) 


where,  to  repeat,  only  connected  graphs  are  considered  explicitly.  The 
final  step  in  the  rearrangement  is  to  group  together  all  terms  in  (4.20) 
associated  with  graphs  of  the  same  basic  type.  Such  graphs  have  the  same 
number  of  vertices  and  the  same  topology  and  differ  only  in  the  multi¬ 
plicity  of  their  bonds.  The  basic  graph  of  such  a  collection  is  the 

graph  with  no  multiple  bonds.  For  instance, 

etc.  are  graphs  of  the  same  basic  t 


basic  graph.  To  see  that  this  regrouping  is  so  vital,  we  introduce  the 


'  * 


. 
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concept  of  labelled  lattice  constant  (g;G)*  which  represents  the  number 
of  weak  embeddings  of  a  labelled  graph  g  on  a  graph  G  whose  sites  are 
themselves  distinguishable .  Here  embeddings  are  counted  separately 
which  differ  not  in  the  particular  sites  of  G  which  are  identified  with 
the  vertices  of  g  but  only  in  the  particular  order  in  which  the  identi¬ 
fication  is  made.  An  important  relation  is  that 

(g;G)  =  ■h.ff-h  (4.21) 

(g;g)* 

where  (g;g)*  is  often  known  as  the  symmetry  number  of  g.  For  example, 


( 


(4.22) 


and. 


)*  =  2  . 


(4.23) 


Let  C^,  i=0,l,2,3. . . ,  denote  the  connected  basic  graphs  or  clus¬ 
ters  as  we  shall  now  call  them.  For  our  purposes  it  will  prove  conven¬ 
ient  to  label  the  clusters  in  order  of  increasing  number  of  bonds,  then 
increasing  number  of  vertices  for  a  given  number  of  bonds  and  finally  by 
some  arbitrary  convention  for  equal  numbers  of  bonds  and  vertices.  If 
g^  is  a  graph  belonging  to  the  same  basic  type  as  CL  ,  then  we  have  the 
important  identity, 


(gr;L)»  =  (c^L)* 


(4.24) 


*  bctA 

. 

,  • 

s:  r »  .j  ex 
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From  (4.21)  we  can  therefore  write, 


(gr;L) 


(C±;L)* 

(gr;sr)* 


(4.25) 


(c^c.)*  (c±;L)* 

(gr;gr)*  (ci;ci)* 


(4.26) 


Again  using  (4.21)  we  obtain  from  (4.25) 


(C. ;C.)* 

(g  -L)  =  -■1-  ■  1  (C.;L)  . 

(gr;gr)* 


Substitution  of  (4.27)  in  (4,20)  then  yields, 


log  Z(T )  =  L  (C±;L)  <f>  ±(K) 
i 


f^OO  =  (C^cp* 


Kr  W'(gr) 

rl  (gr;gr)» 


(4.2?) 


(4.28) 


(4.29) 


In  (4.28)  the  summation  is  over  connected  clusters  only  and  in  (4.29)  it 
is  over  all  graphs  of  the  same  basic  type  as  Ch  , 

(4.28)  constitutes  the  basic  relatior  in  the  cluster  expansion 
method .  The  outstanding  feature  of  this  method  is  that  the  ?4(K)  are 
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not  found  by  considering  (4.29)  which  would  be  just  as  difficult  and  la¬ 
borious  as  the  direct  method  itself.  Rather  we  realise  that  (4.28)  can 
be  applied  to  any  finite  system  as  well  as  the  regular  lattice  L.  This 
means  we  can  write 


(4.30) 


The  ^(K)  are  obtained  by  successively  solving  (4.30)  for  increasing  i 
or  in  other  words,  they  are  found  recursively  from  the  relation, 


?<i+1(K)  =  log  Z(C±+1;K)  - 


T.  (c,;C.  ,)  ?S.(K)  . 
j  -  o  J  h  J 


(4.31) 


It  is  important  to  note  that  the  ^(K)  depend  only  on  and  not  on  the 
system  to  which  the  cluster  expansion  (4.28)  is  being  applied.  Thus  the 
<f>  .  (K)  obtained  from  (4.30)  and  (4.31)  are  the  same  (K)  we  substitute 
in  (4.28)  to  obtain  log  Z  for  any  regular  lattice.  The  only  feature  in 
(4.28)  that  characterizes  the  particular  lattice  under  consideration  is 
the  lattice  constants  (Ch;L). 

To  illustrate  the  procedure,  let  us  consider  i  =  0->3.  (4.30) 


gives 
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log  Z  ( 


,  K)  =jf0(K)  , 


log  Z  ( 


log  Z  ( 


it)  =  2^0(K)  +^X(K)  , 


K)  =  3^o (K)  +  2^_(K)  +  ?S2(K)  , 


(4.32) 


log  Z  (  ^  \4  ,  K)  =  3^q(K)  +  3(^00  +  3^2(K)  +  ?L(K)  . 


The  left-hand  side  of  (4,32)  is  evaluated  in  each  case  as  a  series  in  K 
of  the  form  (3.14).  Since  these  are  the  partition  functions  of  finite 
clusters,  their  calculation  presents  no  problem  in  principle  and  can  be 
performed  on  a  computer.  For  instance  for  the  Hamiltonian 

is  given  by, 


TV-  =  - J(P-|  p  -i-  -H  P^  “  l)  » 


(4.33) 


whose  matrix  representation  has  dimensions  27  x  27.  To  obtain  p  ,  we 
merely  have  to  evaluate  trace  (4r)  and  then  use  the  cumulant-moment 
identities  in  appendix  A  to  obtain  log  Z(  J *\^  ,  K)  as  a  series  in  K. 

We  then  solve  (4.32)  for  ^(K),  i  =  0-?3,  and  to  obtain  their  contribu¬ 
tion  to  log  Z(K)  for  any  given  regular  lattice  we  substitute  in  (4.28) 
using  the  appropriate  lattice  constants. 

A  completely  analogous  procedure  can  be  followed  to  obtain  the 
initial  magnetic  susceptibility.  For  instance,  corresponding  to  (4.28) 


4.15 


and  (4.29),  we  would  have, 


X0(T)  =  2  (Ci:L)  ^(K) 


i 


(4.34) 


^i(K)  =  (ciici)»  r. 

&'r 


(4.35 


(r  +  2)l  (gr;gr)* 


From  (3.18)  and  (4.12)  we  see  that  the  graphs  g  in  (4.35)  have  r  edges 
with  two  of  the  vertices  decorated,  with  crosses. 

Finally,  an  important  feature  of  both  the  <f>  .(K)  and  ^(K)  is  that 
the  lowest  power  of  K  in  their  series  that  can  have  a  non-zero  coefficient 


is  K°  with  b  equal  to  the  number  of  bonds  in  C^.  This  is  easily  under¬ 


stood  since  in  both  (4.29)  and  (4.35)  the  summation  is  over  all  graphs  of 
the  same  basic  type  as  C^.  The  graph  with  the  smallest  number  of  bonds 
in  this  collection  is,  of  course,  Ch  itself  which  contributes  to  that 
power  of  K  which  equals  its  number  of  bonds.  This  means  then  that  to  ob¬ 
tain  the  exact  coefficients  in  (4.28)  and  (4.34)  for  any  infinite  lattice 
up  to  K^,  we  have  to  determine  at  most  the  series  for  ^(K)  and  ^(K)  up 
to  K  corresponding  to  all  graphs  with  a  maximum  of  b  bonds. 


In  summary,  while  the  direct  method  involves  consideration  of  both 


connected.,  disconnected  and  multiline  graphs,  the  cluster  expansion  is 
concerned  only  with  connected  basic  graphs.  Whereas  the  calculation  of 
weights  in  the  direct  method  involves  an  enumeration  of  many  different 
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traces  corresponding  to  any  given  graph,  in  the  cluster  expansion  method 
the  weights  ?H(K)  and  Vh(K)  present  no  such  difficulties  and  can  be  cal¬ 
culated  by  computer.  The  advantages  then  are  obvious. 

The  only  practical  difficulty  arises  in  the  calculation  of  the 

traces  for  the  larger  finite  clusters.  In  our  calculation,  for  instance, 

7 

to  obtain  the  coefficient  of  K  ,  it  is  necessary  to  consider  clusters 
with  up  to  seven  vertices.  The  corresponding  matrices  are  2187  x  2187! 

It  is  quite  impossible  for  the  computing  facilities  available  to  us  to 
manipulate  such  large  matrices.  However,  since  our  interaction  Hamil¬ 
tonian  consists  of  a  sum  of  permutation  operators,  we  were  able  to  make 
use  of  the  powerful  formalism  of  group  theory  applied  to  the  particular 
case  of  the  symmetric  group.  This  essentially  enabled  us  to  perform  a 
similarity  transformation  on  our  matrix  representation  so  that  our  ma¬ 
trices  are  'block  diagonalised '  or  'reduced1.  As  an  indication  of  the 
success  of  this  approach,  for  the  2187  x  2187  matrices  the  largest  ma¬ 
trix  we  in  fact  considered  was  35  x  35. 

In  the  next  chapter,  we  shall  describe  how  group  theory  has  been 
applied  to  this  problem  both  to  determine  log  Z  and  X  . 
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CHAPTER  V 

APPLICATION  OF  GROUP  THEORY 


A.  Concepts  in  Group  Theory 

As  a  preliminary  we  shall  define  a  number  of  the  terms  in  group 
theory  that  we  shall  use  frequently  in  this  chapter.  A  much  more  com¬ 
plete  account  of  these  basic  concepts  can  be  found  in  many  books  on  the 
subject  of  which  Wigner's  (1959)  is  perhaps  the  best  known. 

A  group  is  a  set  of  objects  (the  elements  of  the  group)  among 
which  a  law  of  combination  is  defined  having  the  following  character¬ 
istics. 

a)  It  satisfies  the  closure  property;  namely,  if  A  and  B  are 
elements  of  the  group  and  under  the  law  of  combination  we  obtain 

AB  =  C  ,  (5.1) 

then  C  must  also  be  an  element  of  the  group. 

b)  The  associative  law  must  hold  i.e. 

( AB)C  =  A(BC)  .  (5.2) 

c)  There  exists  among  the  elements  of  the  group  an  identity 
element  E  such  that 


XE  =  EX  =  X 


(5.3) 
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where  X  is  any  element  of  the  group, 

d)  To  every  element  X  of  the  group  there  corresponds  a  reciprocal 
element  X"1  of  the  group  such  that 

XX"1  =  X"3^  =  E  .  (5.4) 

The  above  four  conditions  are  often  known  as  the  group  axioms  or  postu¬ 
lates. 

The  number  of  elements  in  a  group  is  known  as  the  order  of  the 

group. 

If  A  and  B  are  two  elements  of  a  group  such  that 

B  =  X-1AX  ^  (5.5) 

where  X  is  some  element  of  the  group,  then  B  is  said  to  be  conjugate  to 
A.  We  define  the  class  of  A  to  be  the  collection  of  elements  X“1AX 
where  X  runs  over  all  the  elements  of  the  group  i.e.  the  class  of  A  is 
all  those  elements  that  are  conjugate  to  A. 

A  set  of  elements  of  a  group  which  itself  forms  a  group  is  called 
a  subgroup. 

Two  groups  G  and  G*  are  said  to  be  isomorphic  if  the  elements  of 
G  are  in  one-to-one  correspondence  with  the  elements  of  G'  such  that,  if 
X  and  Y  of  G  correspond  to  X*  and  Y1  of  G* ,  then  XY  of  G  corresponds  to 
X'Y'  of  G',  On  the  other  hand  G  is  said  to  be  homomorphic  to  G1  if  each 
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element  of  G'  corresponds  to  at  least  one  element  of  G  such  that  if  X 
and  Y  of  G  correspond  to  X’  and  Y’  of  G*  then  XY  of  G  corresponds  to  X*Y' 
of  G’ .  Evidently  while  isomorphism  is  a  reciprocal  property,  homomorphism 
is  not. 

A  representation  of  a  group  G  is  a  matrix  group  onto  which  the 
group  to  be  represented  is  homomorphic  i.e.  to  every  element  X  of  G  there 
corresponds  a  matrix  D(X)  such  that 

D(X)D( Y)  =  D(XY)  (5.6) 


where  the  law  of  combination  of  the  matrices  is  ordinary  matrix  multi¬ 
plication.  It  is  clear  that  these  matrices  must  be  square  and  the  num¬ 
ber  of  rows  or  columns  is  known  as  the  dimension  of  the  representation. 
If  the  matrix  group  is  isomorphic  to  the  group  it  represents,  then  the 
representation  is  said  to  be  faithful.  Two  representations  which  are 
related  by  a  similarity  transformation  are  said  to  be  equivalent  repre¬ 
sentations.  Such  representations  have  identical  traces  for  correspond¬ 
ing  matrices  and  in  the  formalism  of  group  theory  this  means  there  is  no 
essential  difference  between  them. 

From  two  representations  D(X)  and  D' (X)  which  may  or  may  not  be 
equivalent  it  is  possible  to  form  a  third  representation  D"(X)  given  by, 


D"(X)  = 


D(X)  0  \ 

0  D'(X)/ 


(5.7) 


' 


where  as  before  X  represents  any  element  of  the  group.  D"(X)  and  any 
equivalent  representation  are  said  to  form  reducible  representations  of 
the  group.  In  the  form  (5.7)  we  say  D"(X)  is  a  direct  sum  of  D(X)  and 
D’(X).  Representations  for  which  a  transformation  to  the  form  (5.7)  is 
not  possible  are  said  to  be  irreducible.  If  D(X)  and  D'(X)  are,  however, 
irreducible  representations  then  D"(X)  alone  is  said  to  be  in  completely 
reduced  form.  This  concept  can  be  naturally  extended  to  representations 
of  the  form  (5.7)  but  containing  more  than  two  irreducible  representa¬ 
tions  along  the  principal  diagonal. 

If  X  and.  I  are  two  elements  of  a  group  G  belonging  to  the  same 
class  then  it  is  an  easy  matter  to  show  that 

trace  D(X)  =  trace  D(Y)  (5.8) 

where  D(X)  is  a  representation  of  G.  The  traces  of  the  matrices  cor¬ 
responding  to  a  set  of  elements  of  G  chosen  one  from  each  class  is  known 
as  the  character  table  of  the  representation.  Evidently  equivalent  rep¬ 
resentations  have  the  same  character  table.  Indeed  for  most  purposes  in 
group  theory  all  one  needs  to  know  about  a  representation  is  its  charac¬ 
ter  table. 

One  of  the  important  features  of  the  formalism  is  that  associated 
with  any  group  there  exists  a  finite  number  of  inequivalent  irreducible 
representations  usually  specified  by  their  character  tables.  The  main 
task  for  which  the  formalism  is  used  in  its  applications  is  to  reduce 
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to  its  completely  reduced  form  any  given  reducible  representation  for 
which  perhaps  only  the  character  table  is  known.  How  many  times  each 
irreducible  representation  occurs  in  the  completely  reduced  form  is  de¬ 
termined  by  the  following  important  result  from  the  representation  theory 
of  groups  (e.g.  Wigner  1959 ): 

n(r)  =  i  Z  trX(R)trP(E)  (5.9) 

B  R 


where  y  labels  the  irreducible  representations  and  1  the  reducible  rep¬ 
resentation,  n(y)  is  the  number  of  times  the  y  irreducible  representa¬ 
tion  occurs,  g  is  the  order  of  the  group,  trV(R)  is  the  trace  of  the  ma¬ 
trix  in  the  representation  y  corresponding  to  the  element  R  and  the 
summation  is  over  all  the  elements  of  the  group.  It  should  be  noted 
that  in  (5*9)  the  n(y)  are  obtained  only  from  a  knowledge  of  the  charac¬ 
ter  tables  of  the  representations. 

B.  Application  to  the  Symmetric  Group 

Our  task  in  this  section  is  two-fold.  Firstly  we  wish  to  deter¬ 
mine  the  completely  reduced  form  of  a  particular  reducible  representation 
of  the  symmetric  group  which  means  essentially  using  (5*9).  However, 
for  our  purposes  a  knowledge  of  only  the  number  of  times  a  given  irre¬ 
ducible  representation  occurs  is  not  sufficient.  We  require  in  addition 
to  determine  the  explicit  form  of  the  matrices  in  each  irreducible  rep- 
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resentation  that  occurs.  Then  the  determination  of  the  coefficients  in 
the  finite  cluster  expansions  for  log  Z  and  X  presents  few  difficulties. 

Before  discussing  both  these  steps  in  detail,  we  shall  introduce 
some  of  the  relevant  features  of  the  symmetric  group.  The  elements  of 
the  symmetric  group  TT^  of  the  Nth  degree  are  the  permutations  or  the 
interchanges  of  N  objects.  The  order  of  the  group  is  then  NJ  A  Permu¬ 
tation  which  replaces  1  by  2  by  (Xg  ....  and  finally  N  by  is  de¬ 


noted  by 


(1  2  3  ...  N] 

\  cc-^  a2  #  •  •  cc-j^  j 


An  interesting  property  of  all  Nl  permu¬ 


tations  is  that  they  can  be  resolved  into  a  product  of  cycles.  A  cycle 
(rl  r2  •••  r^)  is  a  permutation  which  replaces  every  element  r^  by  the 
element  following  it  r^  ^ ,  except  the  last  element  r^  which  is  replaced 
by  the  first  element  r^.  The  cycle  (r^  r^  . ..  r^)  is  then  identical  with 


the  permutation 


1  rl  r2  •••  rA. 
r2  r3  *  *  *  rl 


.  As  an  example  the  permutation 

112  3  4  S  ^  o  \ 

3456172/  °an  ^*e  rearran£e6  interchanging  columns  (which  of 
course  does  not  alter  the  permutation)  to  give 


'13  5  2  4  6  7 
,3514672 


from 


which  it  is  obvious  that  the  decomposition  into  cycles  is 


/ 1  2  3  4  5  6  7  ] 
\  3  4  5  6  1  7  2  / 


(1  3  5)(2  467). 


(5.10) 


The  cycle  structure  of  this  permutation  is  denoted  by  (3,4)  the  numbers 
denoting  the  length  of  each  cycle  in  the  permutation.  To  generalise,  the 
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cycle  structure  of  a  permutation  from  TT^  is  given  by  a  partition  of  the 

positive  integer  N  into  an  ordered  set  of  positive  integer  whose  sum  is 

, 

N.  Typically  the  partition  (1  ,  2  ,  3  \ ),  where  we  have 

W 

£  r«  =  K  ,  (5.11) 

r  =  l 

represents  a  permutation  containing  cycles  involving  one  object, 
cycles  involving  two  objects,  . ..,  and  finally  cycles  involving  N 
objects. 

One  can  see  almost  by  inspection  that  permutations  with  the  same 
cycle  structure  belong  to  the  same  class  of  the  symmetric  group.  Indeed, 
the  number  of  classes  in T1j,j  is  equal  to  the  number  of  partitions  of  N 
(the  order  of  the  integers  in  the  partition  is  of  no  consequence)  and  we 
label  each  class  by  its  corresponding  partition.  In  particular  the  in¬ 
teraction  (2,18)  with  which  we  are  primarily  concerned  is  the  transposi¬ 
tion  or  simple  permutation  involving  only  two  objects  and  all  such  simple 
permutations  belong  to  the  class  ofTTN  denoted  by  ( 1^”^  2). 

One  of  the  well  known  results  of  the  general  representation  the¬ 
ory  of  groups  is  that  the  number  of  irreducible  representations  associ¬ 
ated  with  a  group  is  equal  to  the  number  of  its  classes.  In  the  case  of 
the  symmetric  group TT^ ,  it  is  also  possible  to  set  up  a  one-to-one  cor¬ 
respondence  between  the  irreducible  representations  and  the  partitions 
of  N  (e.g.  Wigner  1959).  We  then  also  label  the  irreducible  representa- 
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tions  by  the  partitions,  only  now  we  write  the  numbers  in  descending  order 
to  avoid  ambiguity.  The  character  tables  of  all  the  irreducible  repre¬ 
sentations  ofTTj,,  for  N  =  2,  3  ...  7  are  reproduced  from  Littlewood  (1940) 
in  appendix  C, 

In  the  case  of  a  cluster  of  N  vertices,  the  reducible  representa¬ 
tion  (of  dimension  3^)  that  we  consider  is  the  direct  product  representa¬ 
tion  in  which  the  space  of  the  N  spin  one  systems  is  constructed  as  the 
direct  product  of  N  single  spin  spaces.  It  is  then  possible  to  deter¬ 
mine  almost  by  inspection  the  character  table  of  the  representation.  As 
an  example  of  the  simple  procedure  used  we  shall  consider  the  case  of 

N=3  in  detail.  We  represent  the  three  linearly  independent  single  spin 

3 

one  states  by  a,  0,  y.  The  3  =  2 7  basic  states  of  the  direct  product 

representation  are  given  by  6e.^  where  6e'$  =  a,0  or  y .  Since  there  are 
nine  states  of  the  form  66e.,  the  trace  of  the  simple  transposition  oper¬ 
ator  which  interchanges  the  first  two  symbols  and  hence  all  members  of 
the  class  (12)  is  nine.  Likewise  there  are  three  states  of  the  form  666, 
so  that  the  trace  of  the  operator  that  permutes  all  three  symbols  i.e. 
the  class  (3)  is  three.  The  character  table  of  our  reducible  represen¬ 
tation  ofTT^  is  given  by  , 

(l3)  (12)  (3) 

27  9  3  (5.12) 


In  similar  fashion  the  character  table  of  the  direct  product  reducible 
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representation  oi  Tljj  for  K  <  7  can  easily  be  obtained  and  the  results 
are  given  in  appendix  C. 

We  are  now  in  a  position  to  apply  (5.9)  to  determine  n(y)  Tor 
each  irreducible  representation  ofTT^  for  N  ^  7.  The  calculation  is 
straightforward  and  the  results  obtained  are  reproduced  in  appendix  C. 

To  obtain  the  series  for  log  Z  for  any  given  finite  cluster  of  N  verti¬ 
ces  we  can  write 

=  trace  (Pn)  (5.13) 

=  Y  trace  (P*0  (5.1*0 

V  Y 

where  P^  is  the  matrix  representing  P  in  the  irreducible  representation 
y  of  TT^.  Then  using  the  cumulant-moment  identities  given  in  appendix  A, 
we  immediately  obtain  log  Z  in  the  form  (3.1*0.  The  great  simplification 
introduced  in  the  transition  from  (5*13)  1°  (5.1*0  is  worth  emphasizing. 
For  N  =  7  for  instance,  the  dimension  of  the  matrices  in  (5.13)  Is 
2187  =  3?  whereas  the  largest  matrix  in  (5.1*0  has  a  dimension  of  35. 

There  remains  however,  the  task  of  determining  the  explicit  form 
of  the  permutation  matrices  (i.e.  of  those  matrices  representing  the 
class  (l^-^  2)  of  TT^)  in  any  irreducible  representation.  The  technique 
we  use  is  based  largely  on  a  method  first  introduced  by  Yamanouchi 
(1937).  11  is  essentially  an  inductive  procedure  in  which,  in  order  to 

calculate  the  irreducible  representations  of  TT^  we  assume  a  complete 


} 
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knowledge  of  the  irreducible  representations  of  both  IT  ^  ^  an<^  TT  p_2* 
Let  us  suppose  that  we  wish  to  determine  the  details  of  the  ir¬ 
reducible  representation  D  of  TT  j,.  associated  with  the  partition 
(ff,  f 2 j  ...»  fk)  where 


a)  f-]_  £  *2  ^  ^  ^  fk  ^ 0  » 


b)  i;  f 

r  =  1  r 


=  N  . 


(5.15) 


The  matrices  of  D  which  do  not  contain  the  index  N  i.e,  belong  to  the 

,  -n-  o  -rr  are  formed  by  the  direct  sum  of  the  corresponding 

subgroup  TT  N-]_  of  TT  N  ** 

matrices  in  the  irreducible  representations  of  TT  k  ^  associated  with  the 


partitions , 


( f i~l *  f2*  •••*  * 


(fr  f2-l. 


•  ♦  •  » 


> 


f  ^““1  f  •  •  •  9 


(5.16) 


1  *  ^2  *  0  •  •  t 


fk-D. 
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excepting  these  partitions  in  which  the  f's  are  not  in  decreasing  order. 

The  irreducible  representation  associated  with  (f-^,  .  .  ♦ 

f  .  .^f^.)  will  be  denoted  by  Dr*  We  shall  use  the  suffix  1  r‘  to 

label  the  rows  and  columns  of  D  in  which  the  matrices  occur. 

If  we  consider  only  the  elements  belonging  to  the  subgroup  TT * 
then  Dr  itself  will  take  the  form  of  a  direct  sum  of  the  irreducible 
representations  Drs  of  TT^  2  where  Drg  corresponds  to  the  partition 
(fl»  f 2 ^  •  .  . ,  fr-l,  .  .  .,  fg-1,  .  .  .^f^).  Tlle  Slim  is  for  fixed  r  and 

s  varying  from  1  .  .  .  k  where,  as  before,  cases  in  which  the  f’s  are 

not  in  decreasing  order  after  subtraction  are  omitted. 

These  irreducible  representations  fall  into  three  categories, 

a)  Drr  (only  if  f  fr+1+2)  occurs  once, 

b)  Drs  =  Dsr  arise  twice  if  both  do  not  violate  the  decreasing 
order  of  the  integers  in  the  partitions, 

c)  D  is  possible  but  D  is  not  e.g.  if  s  =  r  -  1  and  f  =  f  . 

'  rs  sr  0  r-i  r 

The  rows  and  columns  within  the  group  previously  labelled  r  are 
now  subdivided  and  denoted  by  the  suffix  rs.  The  rs  rows  and  columns 
indicate  at  which  point  on  the  principal  diagonal  the  matrix  associated 
with  Drs  will  be  placed. 

The  problem  remains  of  calculating  the  permutation  matrices  in¬ 
volving  the  symbol  N.  However,  it  is  sufficient  to  determine  only 

PN-1,M  since 

Pi,N  =  PN-1,N  Pi , N-l  PN-1,N 
(for  all  1  ^  N) 


(5.1?) 
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and  we  have  just  described  how  the  ^  ^  are  obtained. 

Since  ^  ^  commutes  with  all  the  elements  of  TT^  £»  we  make  use 
of  another  of  the  important  theorems  of  the  representation  theory  of 
groups  known  as  Schur's  lemma  which  states: 

If  we  consider  two  irreducible  representations  D'L(X)  and 
D 2(X)  of  the  same  group,  of  dimension^  and  ^  respec¬ 
tively  where  X  is  a  typical  element  in  that  group  and  there 
exists  a  matrix  M  with  rows  and  columns  such  that 


MD*-(X)  =  D2(X)H 


(5.18) 


for  all  X,  then  we  conclude  either 

1)  M  =  0,  provided  ^  4^2* 

2)  M  =  0  if  IT^X)  and  D^(X)  are  not  equivalent  represen¬ 
tations  w±th-£-^  =^/2’ 

3)  if  D^CX)  and  D^(X)  are  equivalent  representations  and 
e1=e2=e.  M  is  a  non-singular  %  dimensional  square 
matrix,  or 

4)  if  D^CX)  and  D^(X)  are  the  same  representation,  M  is 
given  by, 

M  =  XI  . 


(5.19) 


where  X  is  a  constant  and  I  is  the  €>  dimensional  unit 


matrix. 


81*6  j- 
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If  now  we  subdivide  the  rows  and  columns  of  ^  j,j  and  label  them 
1  rs'  just  as  we  did  for  the  elements  of  TT^  £»  we  obtain,  since  ^ 

commutes  with  TT^  £» 


V  P  =  P  V 

rs » rs  rs,pq  rs,pq  pq.pq 


(5.19) 


where  V  „  =  V  _  6  is  the  representation  matrix  for  any  ele- 

rs,pq  rs,rs  rs.pq  ^  v 

ment  of  Tf\T  o  in  the  irreducible  representation  D„_  and  P„_  is  the 

x  s  rs  i 

block  in  the  subdivided  form  of  P^  ^  p  occurring  at  the  intersection  of 
rows  rs  and  columns  pq. 

According  to  Schur's  lemma,  we  must  have  that 


P 


rs  ,pq 


=  0 


(5.20) 


unless  pq  =  rs  or  sr  in  which  case 


and 


P  =  A. 

rs,rs  rs,rs 


I 

rs  ,rs 


(5.21) 


^rs,sr  ~  ^rs.sr  ^rs,sr 


(5.22) 


where  Ia  ^  is  the  unit  matrix  occurring  at  the  intersection  of  rows  a 

2 

and  columns  b.  Since  N  =  I*  we  have  f°r  indices  corresponding  to 

the  types  a)  and  c)  discussed  above  that 


■  >1 
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and  for  type  b) , 


/V, 


rs  ,rs 


2 


\k 

2 

-h 

k 

rs,rs 

rs ,  sr 

(5.23) 


(5.24) 


Thus  we  require  only  the  diagonal  elements  k  to  construct 

I*S  f  rs 

%-l.N-  T°  do  this,  we  first  label  by  P  the  matrix  representing  the 

transposition  of  the  ith  and  jth  symbols  in  the  irreducible  representa- 

N 

tion  D,  It  can  easily  be  shown  that  I,  P.  .  i.e.  the  sum  of  all  the 

i  <  i 
l\l-2 

elements  of  the  class  (1  ~ ^  2)  commutes  with  all  members  of  TTm  and 
therefore  by  Schur’s  lemma  we  can  write 


N 

Z  =  cl  (5.25) 

i<  j  J 


where  c  is  some  constant.  Taking  traces  of  both  sides  we  immediately 
obtain  the  identity, 


MMI  t  =  od 
2 


(5.26) 


where  d  is  the  dimension  of  D  and  T  is  the  trace  of  any  of  the  matrices 

N— 2 

belonging  to  the  class  (1  2)  in  the  representation  D  (given  in  appen 

dix  C).  c  is  then  given  by, 

N(N-1)T 


c 


2d 


(5.27) 
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Similarly  in  the  irreducible  representation  D  ,  we  can  write 


N  -  1 

E  ] 

i<  j  =  1 


,(r) 

ij 


=  c  1 
r  r 


(5.28) 


and  in  the  irreducible  representation  D  , 

±  5 


N  -  2  /  n 

s  =  “i--1 

i  <  j  =  1  J 


rs  rs 


(5.29) 


where  I  and  I  „  are  unit  matrices  of  dimension  d  and  d  respectively, 
r  r  s  r  i  o 

Corresponding  to  (5.2 7)  we  have 


(N-l)(N-2)T 


2d. 


_  (N-2)(N-3) 

°rs  "  2drs  rs  » 


where  and  Tps  are  the  traces  of  the  permutation  operators  in  the  rep 

resentations  D  and  D  . 

i  ±  s 

It  is  clear  from  (5.25)  and  (5.28)  that  the  matrix  given  by, 


= 


N  -  1 

Ei  Pi-n 


N  N  -  1 

E  pij  -  E  Pij  . 

i<j  =  l  J  i <  j  =  1  J 


(5.3D 


is  diagonal  with  sub-matrices  (c-cr)lr  r  along  the  principal  diagonal, 


/ 
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Similarly  from  (5.28)  and  (5.29)  we  can  see  that  ^  defined  by, 


N  -  2 

AN-1  =  Pi,N-l 


(5.32) 


is  a  diagonal  matrix  with  sub-matrices  (c  -  c  )l  _  occurring  on  the 

X  X  S  X  S  j  X  S 

principal  diagonal. 

From  (5.31)  we  can  write  immediately  that 


N  -  2 


ANPW-1,N 


=  u 


i  =  1 


Pi,NPN-l,H  +  1 


(5.33) 


£ 

and,  multiplying  both  sides  by  =  1*  this  becomes 


N  -  2 

PN-1,N  ^  Pi , N-l  +  1  =  ANPN-1 , N 


(5.3zO 


i.e. 


ANPN-1 , N  -  PN-1,NAN-1  =  1 


(5.35) 


Writing  (5.35)  in  terms  of  the  submatrices  labelled  rs,  pq  for  the  two 
cases  pq  =  rs  and  pq  =  sr,  we  obtain 


)P 


rs  ,rs 


-  P 


rs,rs 


Lrs,rs 


1 


(5.36) 
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and 


)P 


rs ,  sr 


P  ( c 
rs,srs  s 


c  ) 

sr' 


0 


(5.37) 


where  c  =  c_„.  Substituting  (5.21)  and  (5.22)  in  (5.36)  and  (5.37)  we 

X*S  Si 

get 


rs ,  rs 


c  -  2cr  +  °rs 


(5.38) 


and  A.„„  =  0  unless  c  -  c  -  ca  +  =  0.  (5*38)  together  with  (5.2^) 

are  the  equations  we  use  to  determine  ^  and  complete  our  calculation 

of  the  irreducible  representation  D.  An  arbitrariness  in  sign  that  will 
arise  in  the  solution  of  (5.24)  just  indicates  that  this  method  may 
yield  several  equivalent  representations. 

We  should  add  that  in  applying  this  technique  we  assume  a  knowl¬ 
edge  of  the  irreducible  representations  of  TT->  and  Tl-^.  They  are  quite 
trivial,  however,  and  can  be  determined  almost  by  inspection.  As  an  ex¬ 
ample  of  this  procedure,  we  shall  consider  the  six-dimensional  irreducible 
representation  (3»1^)  of  TT ~  corresponding  to  the  partition  (3»l»l)»  In 
(5.15)  then  we  have 

fi  =  3;  f2  =  1;  f3  =  1  •  (5.39) 

The  matrices  of  this  representation  belonging  to  the  subgroup  are  the 

direct  sum  of  the  irreducible  representations  (2, l2)  and (3.1) 

labelled  according  to  our  procedure  as  D-^  and  respectively  i.e.  of  the 


form  , 


5.18 


f  (2,12)  0 

\  0  (3,1) 


(5. '*«) 


where  both  submatrices  are  3-dimensional.  Similarly  the  matrices  of  the 
subgroup  TT^  are  of  the  form, 


11  13 


31  32 


11 

13 

31 

32 


0 

0 

(2,1) 

0 


(5.41) 


If  the  irreducible  representations  forming  the  submatrices  of  (5.14), 

D-q  =  (1^)  is  1-dimensional  and  of  the  type  a),  D-^  =  D^]_  s  (2,1)  is 
2-dimensional  and  of  type  b)  and  Dy^  =  (3)  is  1-dimensional  and  of  type 
c).  Also  in  (5.41)  we  have  subdivided  the  rows  and  columns  of  the  matrix 
and  labelled  them  according  to  our  procedure  outlined  above. 

To  determine  ,  we  use  (5.27)  and  (5. 30)  together  with  the  char¬ 
acter  tables  in  appendix  C  to  find  that 


c  =  0, 

c^  =  -2 ;  C-^  =  +2  ^ 

Ufi  =  -3;  =  o  —  c^][»  cyz  ~  * 


A 

y  (5.42) 

J 
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From  (5*38)  we  then  get 


hi, 11  -  +1:  h.3,13  ~  7 :  hi. 31  =  ’I 


(5.43) 


and  finally  using  (5.24)  that 


A  _  4-  V 15  _  a 

*13  .  3-1  "  2  4  ‘  *31.13  * 


(5.44) 


r\ 

A  possible  form  for  in  the  irreducible  representation  (3.u  )  of  TT- 
is  therefore  given  by, 


45 


13 


0 

1 

4 

0 

Vl‘ 

ZT 

0 


V15 

- 


0 


0 

0 

1 

4 

0 

_Vi^ 

4 

0 


31 

0  0 


ill 

4 

0 

1“ 

"  4 

0 

0~ 


_  L 


32 

0 


0 

V15 

4 

0 

1 

"  4 


0 
0 
0 

0  / 


(5.45) 


0  '  ' 


An  exactly  similar  procedure  was  followed  for  all  the  irreducible  rep 
resentations  y  for  which  n(y)  given  in  appendix  C  was  non-zero*  In 
fact,  the  whole  calculation  of  finding  the  explicit  form  of  the  irre¬ 
ducible  representations,  substituting  in  (5.14)  and  then  finding  the 
cumulants  was  performed  in  one  computer  programme. 
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Thus  far  our  discussion  has  been  aimed  at  deriving  the  log  Z 
series  expansions.  However,  in  determining  the  initial  susceptibility 
series,  no  further  traces  have  to  be  calculated.  Only  a  different  lin¬ 
ear  combination  from  (5. 14)  is  required  to  find  ji”  given  by  (4.12).  To 
make  this  point  clear,  it  is  convenient  to  visualise  the  transformation 
from  the  direct  product  representation  to  the  completely  reduced  repre¬ 
sentation  in  two  stages.  The  intermediate  stage  is  a  transformation  to 
a  representation  in  which  the  basic  states  are  simultaneous  eigenstates 
of  S2  =  (S^  +  S2  -l-  .  .  .  +  S^)2  and  Q.  Since 

[P,S2]  =  [P,Q]  =  [P,S+]  =  0  (5.^6) 

where  S+  =  (S*  +  s£  +  .  .  .+  S*)  +  i  (S^  +  s|  +  .  .  .+  s|),  P  will  have  no 

non-zero  matrix  elements  between  states  of  different  S2  and  Q  quantum 

numbers.  In  addition,  for  constant  S2,  the  non-zero  matrix  elements 

will  be  independent  of  the  Q  quantum  number.  In  this  representation,  P 

will  consist  of  a  direct  sum  of  matrices  each  associated  with  a  value  of 

the  quant 'urn  numbers  S  of  S2,  where  S  =  N,  N-l,  N-2,  ...  1,  or  0  with 

P(S)  being  repeated  2S+1  times.  The  notation  P(S)  is  introduced  here 

2 

to  denote  the  matrix  representing  P  in  the  representation  in  which  S  and 
Q  are  diagonal  with  S  fixed. 

The  dimensions  of  P(S)  are  given  in  the  'branching1  diagram  in 
figure  5.1  for  all  S  and  N  <  7,  It  is  constructed  according  to  the  well 
known  rules  for  addition  of  angular  momenta  in  quantum  mechanics.  We 


■ 


Fig.  5.1 


Branching  diagram  for  spin  one  particles.  N  is  the  number  of 
particles,  S  the  total  spin  and  the  numbers  in  brackets  indicate  the 
degeneracy  of  the  corresponding  state  for  fixed  S  . 


5.20a 


7  - 

4 

5 
4  - 

3  - 

2  - 


0 


a) 

(1) — (6) 

\/ 

(D  — (5)— (21) 


(jl) — (4)— (15)— (49) 


(1) —  (3)—  (10)-  (29)- (84) 


(1) - (2) — (6) —  (15)—  (40)- (105) 


(1)— (1) —  (3) —  (6)^ — (15) — (36) —  (91) 


(1)  (1)  (3)  (6)  (15)  (36) 


4 


7 


N 


5.21 


then  have  in  this  representation,  from  (4.12),  that 

=  -£4  E  traoeipr  (S)[s(S  +  l)(2S  +  X)  (5>7) 

r  3n+i  s  =  o  l  -3 

where  we  have  already  summed  (for  fixed  S)  over  all  Q  from  -S  to  -i-S. 

The  determination  of  trace  ^ Pr”2 3 ( S  is  further  simplified  by  transform¬ 

ing  to  the  completely  reduced  form.  For  a  given  value  of  N,  this  actu¬ 
ally  entails  distributing  the  known  numbers  of  the  irreducible  repre¬ 
sentations  among  the  various  P(S)  for  S  =  0,  1,  2,  ...»  and  N.  This  is 
achieved  by  a  process  of  elimination  satisfying  the  following  constraints. 

1)  The  irreducible  representation  y  must  occur  a  total  of  n(y) 
times  i.e.  if  it  is  associated  in  the  completely  reduced  form  with 

S  =  S-p  S £»  ...,  then  we  must  have 

n(y)  =  £  (2Sr  +  D  •  (5.W 

r  =  1 

2)  The  dimensions  of  the  irreducible  representations  associated 
with  a  particular  value  of  S  must  add  up  to  the  dimensions  of  the  P(S) 
as  given  in  the  branching  diagram. 

3)  For  fixed  N,  the  one  dimensional  irreducible  representation 
n(N)  will  occur  once  in  the  reduced  form  of  P(S)  for  S  =  N,  N  -  2, 

N  -  4....  etc.  We  must  have  therefore 


■  ■ 

1 


' 


* 
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N 

2 

n(N)  =  5Zl  (2N  -  4r  +  1)  (5.49) 

r  =  0 


where  the  summation  in  (5*49)  is  over  all  positive  integer  values  of 

r«f. 

This  last  condition  requires  some  explanation.  We  shall  illustrate 

for  the  case  of  W  =  4  but  the  argument  is  exactly  the  same  for  all  values 

2 

of  N.  To  find  eigenfunctions  of  S  and  Q  which  might  correspond  to  the 
irreducible  representation  (4),  we  require  combinations  of  products  of 
the  a(i),  $(i)»  and  y(i),  the  single  particle  eigenfunctions  of  S^,  that 
are  symmetric  under  all  the  operations  of  the  symmetric  group.  For  Q=4, 
there  is  only  one  possibility  namely  aacca  where  the  labelling  of  the 
spins  has  been  omitted.  This  must  also  be  an  eigenstate  of  and  hence 
(4)  must  be  associated  with  the  S  =  4,  Q  =  4  eigenstate.  Because  of  the 
independence  of  P(S)  on  Q,  we  therefore  associate  nine  of  the  occurrences 
of  (4)  with  the  S  =  4,  Q  =  +4~^  -4  eigenstates.  For  Q  =  3»  there  is 
likewise  only  one  possible  completely  symmetric  eigenstate  which  must 
also  be  then  a  simultaneous  eigenstate  of  S  with  S  =  4.  It  is  of  the 
form  ctcta3  +  aa3oc  +  aBaa  -1-  |3aaa.  For  Q  =  2,  there  iie  two  linearly  inde¬ 
pendent  eigenstates  namely  omy  +  aaya  -i-  ayaa  +  yaaa  and  aa33  +  33aa  + 
a3oc3  +  3cc|3ct  +  a33oc  +  3oca 3.  Either  one  of  these  or  a  linear  combination 
of  both  (it  does  not  matter  which  for  our  purpose)  must  be  a  simultane- 
ous  eigenstate  of  S  with  S  =  4  and  the  other  state  or  the  other  lin- 


bftwrv-is'TCO:?  3  : 
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early  independent  combination  must  be  an  eigenstate  with  S  =  2,  We  have 
therefore  an  additional  five  occurrences  of  the  irreducible  representa¬ 
tion  (4)  associated  with  S  =  2.  For  Q  =  1,  there  are  again  only  two 

possibilities  namely  ocatfy  +  aay3  +  .  .  .  a  total  of  —  =  12  different 

2  I 


orderings  of  the  four  symbols  and  a333  +  3a33  +  33a3  +  333a.  Hence  (4) 
can  not  be  associated  with  S  =  1.  For  Q  =  0,  however,  there  are  three 

linearly  independent  eigenstates  given  by  3333.  aayy  +  .  .  .  a  total  of 

4i/  4i 

=  6  different  orderings  and  a33y  -l-  ...  a  total  of  ~  =  12  differ- 

c  *4  I  us  \ 

ent  orderings.  One  of  these  states  or  one  linear  combination  must  also 

2 

be  an  eigenstate  of  S  with  S  =  0  and  so  we  have  one  more  occurence  of 
(4).  The  total  number  of  occurrences  is  fifteen  and  from  appendix  C  we 
see  that  this  concurs  with  (5.49)  for  N  =  4. 

Finally,  as  an  example  of  how  all  these  constraints  are  used,  we 
shall  consider  the  case  of  N  =  5.  From  appendix  C,  we  know  that 


n(5)  =  21  with  dimension 


n(4,l)  =  24  " 

n(3.2)  =  15  " 

n(22,l)  =  3  " 

n(3,l2)  =  6  " 


(5.50) 


and  all  other  n(y)  =  0.  This  distribution  has  to  be  fitted  to  the  block 
diagonalisation  given  in  figure  5.1  for  N  =  5  i.e. 


•r 

•/ 


' 
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for  S  =  5, 

"  S  =  4, 

"  s  =  3, 

"  S  =  2, 

"  S  =  1, 

»  S  =  0, 


the  dimension  of  the  block  is  1  , 


ii  n 


II  I! 


II  II 


II  II 


It  It  II  II 


»  "  4  , 

"  "  10  , 

"  "  15  , 

"  "  15  , 

6  . 
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From  constraint  3)  the  representation  (5)  occurs  eleven  times  in  the 
S  =  5  block,  seven  times  in  the  S  =  3  and.  three  times  in  the  S  =  1.  The 
total  Qf  21  agrees  with  (5.50).  The  only  possible  location  for 

r\ 

the  (2fl)  representation  is  the  S  =  1  block.  Then  from  (5.48)  we  must. 
have  (3,1  )  distributed  so  that  five  of  its  occurrences  are  in  the  S  =  2 
block  and  one  in  the  S  =  0  block.  Nine  of  the  (4,1)  must  go  to  the  S  =  4 
block  and  so  on.  By  this  procedure,  we  obtain  a  unique  distribution  of 
the  irreducible  representations  so  that  all  the  contrainsts  are  satis¬ 
fied  i.e. 


for  S  = 
"  S  = 
"  S  = 

"  s  = 
"  s  = 
"  s  = 


5,  the  irreducible  representations  are 


4,  " 


ti 


ti 


3»  " 


ii 


ti 


ii 


2,  " 
1,  ” 
0,  " 


ti 

ii 

ii 


ii  it 

ti  n 


ti 


it 


(5)  . 

(4.1)  , 

(3.2) +(4,l)+(5)  , 
(3>12)+(3,2)+(4,1)  , 

(3.2) +(22,l)+(4,l)-h(5) 
(3. I2)  . 


H5.52) 

9 
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where  the  suiomation  sign  in  (5.52)  means  direct  sum. 
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For  each  N  ^7,  it  was  possible  by  this  procedure  to  associate 
the  irreducible  representations  uniquely  with  particular  values  of  S. 

The  distributions  are  given  in  appendix  C.  Therefore,  essentially  the 
same  data  is  used  to  find  }i”  as  was  required  in  (5.1^)  to  determine  p.  , 
The  susceptibility  series  for  the  various  clusters  is  then  obtained  by 
making  use  of  the  identities  in  appendix  A  giving  the  in  terms  of 
the  ji'g  with  s  r.  Again  the  calculation  was  performed  on  the  computer. 
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CHAPTER  VI 

DISCUSSION  AND  ANALYSIS  OF  SERIES 


A.  Discussion  of  Series 

For  the  spin  one  model  of  ferromagnetism  defined  by  (3.2),  using 
the  procedures  outlined  in  (4.28-4.35)  together  with  the  group  theoretic 
simplifications  described  in  the  last  chapter,  it  is  possible  to  deter¬ 
mine  exactly  the  coefficients  of  Kr  for  r  ^  7  in  the  high  temperature 
series  expansions  of  the  zero  field  partition  function  and  the  initial 
reduced  susceptibility.  The  series  ^(K)  and  y^(K),  the  weights  associ¬ 
ated  with  the  cluster  Ch  in  the  partition  function  and  susceptibility 
series  respectively,  are  reproduced  in  appendix  D  together  with  all  the 
linear  combinations  embodied  in  (4.30)  and  illustrated  for  a  few  cases 
in  (4.32). 

A  number  of  features  of  the  series  is  worth  noting.  While  it  is 
necessary  to  consider  clusters  containing  a  maximum  of  seven  bonds  we  do 
not  need  to  perform  calculations  on  all  such  graphs.  For  certain  classes 
of  clusters  containing  b  bonds,  all  coefficients  of  Kr  in  ^(K)  or  y^(K) 
are  zero  for  r  <  b  -f  A.  where  the  positive  integer  k  is  known  as  the 
entry  parameter.  Thus,  for  example,  all  seven  bonded  clusters  for  which 
both  (f>  ^(K)  and  ^(K)  have  k  =  1  need  not  be  considered.  Similarly  all 
six-bonded  clusters  with  k  -  2  in  both  cases  make  no  contribution  up  to 
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K'  and  so  on.  The  entry  parameter  of  a  cluster  may  not  be  the  same  for 
both  the  ^(K)  and  the  yn(K)  series  and,  as  we  see  in  the  appendix  D, 
for  many  clusters,  only  yy(K)  was  calculated. 

To  see  how  we  can  determine  the  entry  parameter  of  any  cluster, 
we  recall  how  the  various  terms  in  the  series  ^(K)  and  Ij-H(K)  are  associ¬ 
ated  with  graphs  of  the  same  basic  type  as  CH.  In  the  case  of  the  ^(K), 
if  contains  b  bonds  the  coefficient  of  K  or  the  primary  weight  of 
^•(K),  as  it  is  called,  is  given  by  the  trace  of  spin  operators  corre¬ 
sponding  to  CH  itself.  In  our  discussion  of  the  direct  method  we  noted 
how  graphs  with  at  least  one  vertex  of  degree  one  had  zero  weight  i.e. 
the  corresponding  spin  variable  products  were  traceless.  Hence  we  can 
conclude  that  all  clusters  containing  one  vertex  of  degree  one  will  have 
no  primary  weight  and  an  entry  parameter  of  one  in  the  corresponding 
^(K).  For  example,  as  inspection  of  appendix  D  will  show,  ^(K)  cor¬ 


responding  to  the  cluster 


has  no  primary  weight. 


In  this 


case  the  first  non-zero  coefficient,  that  of  K  ,  will  be  associated  with 


the  multiline  graph 


degree  one  e.g. 


.  Similarly  graphs  with  two  vertices  of 


will  have  no  primary  weight  or  secondary 


weight  and  the  first  non-zero  coefficient,  which  occurs  for  K?,  is  associ¬ 


ated  with  the  multiline  graph 


Indeed  for  any  cluster  Cj_ 
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we  can  say  that  the  entry  parameter  of  its  yT(K)  is  at  least  equal  to 
the  number  of  vertices  of  degree  one  in  the  cluster. 

For  thel^(K),  however,  the  primary  weight  corresponds  to  the 
cluster  decorated  with  two  crosses,  these  crosses  indicating  the  ad- 
dition  of  S  variables  on  the  sites  they  occupy.  If  the  cross  is  placed 
on  a  vertex  of  degree  one,  then  the  product  of  spin  variables  on  that 
site  will  no  longer  be  traceless.  As  an  example  while  ^(K)  correspond¬ 


ing  to  the  cluster 


has  an  entry  parameter  of  two,y^(K) 


has  a  non-zero  primary  weight  associated  with  the  decorated  cluster 
* 

,  A  simple  extension  of  the  arguments  presented  leads  to 

* 

the  conclusion  that  the  entry  parameter  of  yy(K)  is  at  least  two  less 
than  the  number  of  vertices  of  degree  one  in  Cy. 

Inspection  of  the  series  in  appendix  D  will  show  that  these  entry 
parameter  rules  for  both  the  ^(K)  and  the  yy(K)  are  verified  by  calcu¬ 
lation  in  all  cases.  The  reason  for  the  ’at  least’  in  the  statements  of 
both  rules  is  that  we  have  no  guarantee  that  the  primary  weight  of  a 
graph,  which  we  expect  from  our  considerations  to  be  non-zero,  may  not  in 
fact  turn  out  to  be  zero  because  of  the  particular  topology  of  the  graph. 

For  example  the  yy(K)  corresponding  to  the  clusters 


etc.  which  we  expect  to  have  zero  entry 


' 


6,4 


parameters  actually  have  no  primary  weights.  However,  we  have  ennumerated 
our  clusters  on  the  basis  that  the  entry  parameters  are  either  equal  to 
or  two  less  than  the  number  of  vertices  of  degree  one  to  ensure  that  all 
those  contributing  up  to  are  included. 

An  important  feature  of  these  rules  is  that  they  provide  a  signifi¬ 
cant  check  on  our  calculations.  If  anycp^(K)  orCp/^CK),  whose  calculation 
involves  many  of  the(jb^(K)  and(^i(K)  previously  evaluated,  does  not  fol¬ 
low  the  rules  then  we  are  assured  a  mistake  has  occurred.  It  is  a  necessary 
though  not  sufficient  condition  for  the  correctness  of  our  calculation 
that  these  rules  be  followed  in  all  cases. 

Of  all  the  clusters  we  considered,  only  the  linear  chain  of  seven 

bonds  i.e,  ^/\^//%^/\^/*  contains  eight  vertices.  To  calculate 

the  corresponding  y^(K)  by  the  cluster  expansion  method  would  require 

o 

consideration  of  matrices  of  dimension  3  and  the  considerable  labour  in 
evaluating  the  irreducible  representations  etc.  of  TTg  just  for  the  sake 
of  one  graph.  To  avoid  all  these  difficulties  we  note  only  the  primary 

weight  of  the  ^(K)  given  graphically  by  has  to 

be  evaluated  and  this  can  be  done  directly  for  all  linear  chains.  For  a 
linear  chain  of  b  bonds,  the  coefficient  in  the  corresponding  ^(K)  of 
Kb/3bbl  is  given  from  (4.12)  and  (4.35)  hy 

§  2  traCe[(P12-i)(P23-i)...(Pbib  +  i-^)S*S^+1]  (6.1) 


3  MO 


6.5 


where  the  summation  in  (6.1)  is  over  all  the  b!  arrangements  of  the  b 
brackets.  Let  us  first  consider  the  product  P]_2P23P34  •**^b  b  +  1*  Pn 
the  direct  product  representation  the  only  non-zero  diagonal  elements  of 
this  permutation  operator  correspond  to  the  completely  symmetric  states 
given  symbolically  by  a(l)a(2)  ...  a(b  +  l),  |3(l)|d(2)  ...  3(b  +  l)  and 
y(l)y(2)  ...  y(b  +  l)  i.e.  there  are  three  diagonal  elements  equal  to 
unity.  The  eigenvalues  of  in  these  states  are  +1,  0+1  respec¬ 

tively  so  that 

trace(P12P23P34  •  •  •  pb,b+lSlSb+l)  =  2  •  (6.2) 

All  the  b!  arrangements  of  permutation  operators  in  (6.2)  will  have  the 
same  trace  so  that  the  contribution  to  (6.1)  of  this  product  is  4bt/3. 

By  similar  arguments  to  those  given  above,  we  can  show  that  all  other 
contributions  to  the  expansion  of  (6.1)  in  terms  of  products  of  permu¬ 
tation  operators  is  zero.  Each  product  contains  a  part  of  the  form, 


P12P23 


P  .Sf 
r,r+l  1 


(6.3) 


where  the  number  of  permutation  operators  varies  from  zero  to  b-1.  In 
the  direct  product  representation,  the  permutation  product  in  (6.3)  has, 
as  before,  three  non-zero  diagonal  elements  each  equal  to  unity.  With 
the  inclusion  of  the  single  Sz  matrix,  (6.3)  then  becomes  traceless. 

We  therefore  conclude  that  the  coefficient  of  Kb/31:)bt  in  3-4  (K) 


■ 

.  ' 


6.6 


corresponding  to  the  linear  chain  of  b  bonds  is  given  by  4b ! /3 .  This 
result  agrees  with  all  our  cluster  expansion  calculations  for  b  -<C  6  and, 
as  shown  in  appendix  D,  we  conclude  that  for  b  =  7  the  coefficient  is 
equal  to  6720. 

The  results  described  above  have  been  combined  according  to  (4.28) 
and  (4.35)  to  obtain  the  logarithm  of  the  partition  function  per  site 
and  the  reduced  initial  susceptibility  per  site  respectively  on  the 
face-centred  cubic  lattice.  The  reason  for  choosing  this  particular 
lattice  is  that,  because  of  its  high  coordination  number  (twelve),  so 
many  configurations  of  spins  occur  for  each  power  of  K  that  the  conver¬ 
gence  of  the  series  is  expected  to  be  most  rapid  for  this  lattice. 

If  we  write  each  thermodynamic  quantity  in  the  form, 


(6.4) 


the  coefficients  &n  for  log  Z/N  and  kT^/Nm2  are  reproduced  in  Table 
6.1.  Also  included  in  that  table  are  the  series  coefficients  a^  for 
the  specific  heat  per  site  c/hk  and  the  entropy  per  site  in  the  form 
(Soo-S)/Nk  obtained  from  log  Z/N  by  using  (3.6)  and  (3.5)  respectively. 
(S^/Nk  =  loge3  is  the  entropy  per  site  in  the  limit  of  infinite  tem¬ 
perature.  ) 

Many  of  the  lattice  constants  (C, ;L)  required  in  the  derivation 
of  the  above  series  are  listed  in  Domb  (i960).  For  the  last  term  in 
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6.7 


the  susceptibility  series,  however,  thirty  of  the  contributing  clusters 
do  not  have  their  lattice  constants  listed  in  the  literature.  All  of 
these  were  open  configurations  and  their  constants  can  be  obtained  by 
the  'symbolic  equation'  method  described  by  Domb  (i960).  As  an  example 


let  us  consider  the  cluster 


constant  on  the  face-centred  cubic  lattice  by 


and  represent  its  lattice 


can  be  formed  from  the  hexagon  by  adding  a  single  bond  to  any  one  of  its 
vertices  which  can  be  done  a  total  of  sixty  distinct  ways.  However, 
this  added  bond  may  be  a  chord  or  diagonal  of  the  hexagon  so  that  the 


clusters 


may  be  formed,  the 


number  of  times  per  site  that  these  two  configurations  occur  is  twice 
the  corresponding  lattice  constants  since  the  added  bond  could  have  been 
originally  attached  to  either  of  the  two  vertices  it  joins.  We  there¬ 
fore  set  up  the  following  symbolic  equation, 


(6.5) 


linking  the  unknown  lattice  constant  on  the  left-hand  side  to  a  linear 
combination  of  lattice  constants  for  the  face-centred  cubic  lattice 
listed  in  Domb  (i960).  Hence  we  obtain  in  this  case  that 


/  \ 


( 


)  =  51,z^8  . 


(6.6) 
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The  numerical  value  of  the  lattice  constants  of  all  thirty  required  open 
configurations  are  listed  in  appendix  E  together  with  the  symbolic  equa¬ 
tions  from  which  they  were  derived. 

B«  Analysis  of  Series 

To  locate  the  singularities  in  the  thermodynamic  functions  such 
as  magnetic  susceptibility  etc.  to  which  the  series  in  table  6.1  represent 
approximations  for  T  <  T  <  °o,  two  well  known  extrapolation  techniques  are 

V 

used.  Both  the  ratio  method  introduced  by  Domb  and  Sykes  (1957)  a^d  the 
Pade  approximant  method  resurrected  by  Baker  (1961)  are  attempts  to  ex¬ 
trapolate  to  the  singular  behavior  of  the  physical  functions  near  the 
critical  temperature. 

In  the  ratio  method ,  if  the  dominant  singularity  of  a  function 

f(x)  occurs  at  some  value  x  ,  it  is  assumed  that  near  x  f(x)  is  of  the 

c  ^ 

form, 


f(x)  «  a(i  -  #-rY  , 


(6.7) 


where  A  is  a  constant.  Then  the  coefficients  b^  in  the  series  expansion 
of  f(x)  in  powers  of  x  should  tend  in  the  limit  of  large  n  to  the  bi¬ 
nomial  coefficients  in  the  expansion  of  (6.7).  The  limiting  form  of  the 
ratio  of  coefficients  is  given  in  that  case  by, 


(6.8) 


- 
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Therefore,  if  the  ratio  bn/t>n  is  plotted  against  l/n,  the  graph  should 

become  linear  for  sufficiently  large  n.  The  projected  intercept  on  the 

b  /b  ->  axis  will  give  l/x_  and  the  limiting  slope  will  equal  (y-I/x  . 
n  n-i  c  “ 

i 

Similarly  a  plot  of  (b  /b  9)2  against  l/n  should  approach  linear  behav- 
ior  with  intercept  l/xc  and  slope  (y-l)/xc.  Also  from  (6.8)  if  we  plot 
nbn/(nH-y  -  l)bn  ^  against  l/n,  the  graph  should  approach  a  horizontal 
line  with  an  intercept  of  l/x  .  This  last  test  is  used  if  some  value  of 
Y-l  is  postulated. 

It  is  clear  that  this  method  can  only  be  applied  to  a  series  of 
all  positive  coefficients.  This  means  in  our  case  that  only  the  sus¬ 
ceptibility  series  can  be  analysed  by  the  ratio  method.  All  three  plots 
of  the  bn  given  by  (6.4)  and  Table  6.1  for  the  reduced  susceptibility 
are  shown  in  figure  6.1.  For  the  first  two,  the  linear  region  seems  to 
be  reached  and  we  estimate  that  l/Kc  =  3 .075  +  0.01  or  KQ = 0.325 +0.001 
and  the  critical  index  of  the  susceptibility  y  =  1.27  +  0.03.  There  is 
a  fashion  in  this  field  to  postulate  simple  fractions  for  critical  indices 
and  if  we  yield  to  this  and  postulate  y  =  5/4,  the  third  plot  of 
nbn/(n ■+• £)b  ^  against  l/n  does  seem  to  approach  linear  behavior  with 

an  intercept  agreeing  within  graphical  accuracy  with  the  values  quoted 
above . 

It  should  be  emphasized  that  in  the  first  instance  it  is  a  signifi¬ 
cant  assumption  to  postulate  the  form  (6.7)  for  the  susceptibility  in  the 
critical  region  for  K  <  Kc  (i.e.  above  the  critical  temperature).  How¬ 
ever,  the  remarkable  accuracy  of  the  predictions  of  such  a  simple  tech- 


.n 


Fig.  6.1 


Ratio  plots  for  the  spin  one  exchange  model  initial  suscepti¬ 
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nique  as  the  ratio  method  strongly  suggests  that  (6.7)  may  in  fact  rep¬ 
resent  the  asymptotic  behavior  of  the  susceptibility. 

The  Pade  approximant  method  can  be  applied  to  any  series  whether 
or  not  the  coefficients  are  all  positive.  The  Pade  approximant 

approximates  a  function  f(x)  by  the  ratio  of  two  polynomials 

P(x)  =  £  p_,xr  of  degree  N  and  Q(x)  =  £  q  xs  of  degree  M  with  qn  =  1, 

I*  *  s  s  u 

The  coefficients  of  the  polynomials  are  chosen  to  make  the  expansion  of 
the  [N,M]  approximant  agree  exactly  with  the  expansion  of  f(x)  to  the 
first  M+N  +  l  terms.  The  identity, 

f  (x)Q(x)  -  P(x)  =  ax^'^+^  +  bx^1*2  +  . . . ,  (6.9) 

defines  the  approximant  and  by  equating  coefficients  of  powers  of  x  in 

(6.9)  one  can  solve  a  set  of  M+N  +  l  simultaneous  equations  to  obtain 

the  coefficients  p  and  q  .  Such  an  operation  can  be  readily  performed 

I*  s 

on  a  computer. 

The  convergence  of  Pade  approximants  has  been  studied  by  Baker ,  Gam- 
mel  and  Wills  (1961)  and  Baker  (1965).  We  will  not  consider  the  prob¬ 
lem  here  except  to  say  that  there  is  considerable  evidence  not  all  mathe¬ 
matically  rigorous  to  suggest  that  in  most  problems  of  the  type  we  are 
considering  they  converge  very  well  when  sufficient  terms  of  the  series 
are  known  (Fisher  1965). 

The  Pade  approximant  can  be  applied  in  the  following  four  situa¬ 


tions. 
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a)  If  the  function  f(x)  is  of  the  form  (6.7) .  then  it  is  con¬ 
venient  to  examine  the  approximants  to  the  logarithmic  derivative  of  the 
series  for  f(x)  given  by, 


(6.10) 


which  now  has  a  simple  pole.  We  expect  the  smallest  positive  zero  in 
the  denominator  of  the  approximant  to  correspond  to  x  and  that  the  resi- 

V 

due  at  that  point  will  give  the  critical  index  y. 

b)  If  some  precise  independent  estimate  of  the  critical  point  xq 
is  known,  then  evaluating  the  approximants  of  the  transformed  series 
(6.10)  multiplied  by  (x  -x)  at  x=x  should  yield  further  estimates  of  y. 

c)  If  f(x)  is  the  form  (6.7)  and  some  independent  estimate  of  y 
is  known,  then 


(6.11) 


x 

c 


and  the  smallest  positive  zero  in  the  denominators  of  the  approximants 

to  (6.11)  should  correspond  to  x  .  Also  the  residues  at  that  point 

l/y 

should  reproduce  A 

d)  Baker  et  al.  (1967)  have  suggested  that  the  approximants  to 


f(x)  = 


the  expression  , 


(6.12) 


■  ' 


■ 

. 

i  J- 


-  ■  -  tr" 


6.12 


evaluated  at  x  =  x  ^  reproduce  values  of  y  which  are  very  insensitive  to 

c 

the  choice  of  x  . 

c 

All  four  procedures  were  applied  to  the  susceptibility  series. 

In  case  a)  the  results  were  inconclusive.  For  the  [M, N]  approximants 
for  which  M-»-N  =  6,  the  values  of  K  varied  from  0.320  to  0.329  and  y 
from  1.17  to  1.39.  While  they  are  consistent  with  the  much  more  accurate 
ratio  method  results,  the  scatter  is  much  too  large  to  draw  any  definite 
conclusions.  Considerable  improvement,  however,  occurs  in  method  b)  if 
we  use  the  value  of  K,  =  0.325  suggested  by  the  ratio  method.  The  values 

V 

of  y  given  by  all  the  approximants  for  which  M+N  =  5  and  6  are  repro¬ 
duced  in  Table  6.3.  Here  the  value  of  y  indicated  is  1.27  +  0.02,  pro¬ 
vided  we  regard  the  value  of  the  l_3»3j  approximant  as  spurious.  It  is  a 

matter  of  experience  that  occasionally  one  or  two  approximants  yield  re¬ 
sults  very  different  from  all  others.  This  can  often  occur  if  there 

exists  in  those  approximants  another  singularity  close  to  the  physically 
relevant  one.  This  value  of  y,  then,  is  impressively  consistent  with 
the  independent  estimate  of  the  ratio  method. 

In  method  c) ,  we  chose  the  value  of  y  to  be  5/^»  and  the  roots  of 
the  resulting  approximants  for  M  +  N  =  5*  6  and  7  are  reproduced  in 
Table  6.2.  There  is  very  little  scatter  in  these  results  which  yield  a 
mean  value  of  K  =  0.325  +  0.0005.  The  mean  value  of  the  residues  sug- 
gests  that  A  =  0.230  +  .01. 

Because  of  the  double  differentiation  involved  in  method  d),  it 
is  more  suitable  for  longer  series,  than  we  have  available  in  this  case. 


. 
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Table  6.2 

The  Pade  approximant  estimates  of  from  (Xq)^^ 


m/n 

1 

2 

34567 

0 

O.3233  O.3252  0.3244 

1 

0.3233  O.3235  0.3246 

2 

0.3259  0.3252  0.3245 

3 

0.3219 

0.3252  0.3273 

4 

0.3221 

0.3230 

0.3243 

5 

0.3249 

0.3247 

6 

0.3247 

Table  6.3 

-  d 

The  [M,N]  Pade  approximant  estimates  of  y  from  [(K^-K)  ^logX0Jj/_^ 

with  K  =  0.3252 
c 

M/N 

1 

2 

3  4  5  6 

0 

1.28  1.27 

1 

1.29  1.27 

2 

1.24  1.27 

3 

1.25 

1.41 

4 

1.28 

1.27 

6.13 


It  is  not  surprising  that  for  K  =  0.325,  the  values  of  y  obtained  ranged 
from  1.17  to  1.39  for  the  approximants  for  which  M  +  N  =  5.  It  is  inter¬ 
esting  to  note,  however,  that  the  mean  value  of  the  estimates  gives 
Y  =  1.27. 

The  very  close  agreement  between  the  results  of  these  two  inde¬ 
pendent  extrapolation  techniques  suggests  strongly  that  just  above  the 
critical  temperature  the  susceptibility  is  of  the  form  (6.7)  with  the  two 
main  critical  parameters  given  by, 


K  =  0.325  ±  0.001 
y  =  1.27  ±  0.03 


(6.13) 


We  emphasize  that  the  errors  bars  in  (6.13)  are  not  based  on  rigorous 
mathematics  but  rather  are  confidence  limits  we  have  set  on  the  results 
on  the  basis  of  the  consistency  we  have  obtained.  Perhaps  one  of  the 
disadvantages  of  the  Pade  approximant  is  that  it  is  not  possible  to  give 
rigorous  mathematical  bound  to  its  rate  of  convergence  except  in  very 
specialized  circumstances.  While  the  series  we  have  obtained  con¬ 
tain  a  relatively  small  number  of  terms,  they  incorporate  a  considerable 
amount  of  information  about  the  details  of  the  lattice  and  how  the  cor¬ 
relations  begin  to  increase  as  the  temperature  approaches  critical.  This 
is  reflected  in  the  contribution  of  larger  graphs  to  the  higher  order 
terms.  When  regarded  from  this  point  of  view,  the  accuracy  of  the  re¬ 
sults  (6.13)  is  perhaps  not  so  unreasonable  and  a  good  indication  of  how 
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the  ultimate  critical  behavior  is  already  incorporated  in  these  trun¬ 
cated  series. 

As  has  been  found  for  most  other  models,  similar  analysis  of  the 
specific  heat  series  yields  no  consistent  results.  This  is  perhaps  a 
reflection  of  the  fact  that  the  specific  heat  singularity,  if  it  exists, 
is  much  weaker  than  for  the  susceptibility  and  requires  much  longer  ser¬ 
ies  to  predict  it  precisely.  Only  very  recently,  for  the  three-dimen¬ 
sional  Ising  model  have  precise  predictions  of  the  specific  heat  behav¬ 
ior  been  made  on  the  basis  of  a  series  with  fourteen  terms  (Hunter  1967). 
Again  Baker  et  al.  (1967)  on  the  basis  of  considerations  we  shall  not 
discuss  here  have  indirectly  predicted  the  critical  behavior  of  the 
specific  heat  for  the  spin  one  half  Heisenberg  model.  They  find  that  the 
specific  heat  is  not  even  infinite  at  the  critical  temperature  but  rather 
has  an  infinite  slope.  If  such  behavior  is  characteristic  of  the  spe¬ 
cific  heat  in  quantum  mechanical  models,  it  is  not  surprising  that  with 
only  six  terms  we  can  not  obtain  any  consistency. 

It  is  possible  however  to  evaluate  the  Pade  approximants  to  the 
entropy  series  to  find  the  fraction  of  the  entropy  change  occurring 

above  T,.  This  is  an  experimentally  interesting  quantity  which  like  the 
o 

critical  index  is  independent  of  the  interaction  constant  J.  The  mean 
value  we  obtained  is  given  by, 

S  -S 


00  c 


=  0.33  +  0.02  . 


(6.14) 
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The  principal  results  of  our  calculation  on  the  spin  one  exchange 
model  are  displayed  in  comparison  with  other  relevant  models  in  Table  6.4. 
Many  of  these  results  can  be  found  in  the  review  article  of  Domb  and 
Miedema  (1964).  It  is  interesting  to  note  that  the  spin  one  exchange 
model  becomes  disordered  at  a  lower  temperature  than  the  spin  one  Heisen¬ 
berg  model,  as  is  to  be  expected  from  the  higher  degeneracy  of  the  ground 
state  of  a  pair  of  spins  interacting  with  the  pure  exchange  interaction. 
The  fraction  of  the  total  entropy  change  above  T  is  directly  comparable 
among  the  various  models.  We  see  that  the  Ising  model  differs  markedly 
from  all  the  quantum  mechanical  models  in  this  respect,  whereas  the 
(Soo-Sc)/S00  for  the  spin  one-half  and  spin  one  exchange  models  are  equal 
within  1  experimental' error , 

The  quantity  of  greatest  theoretical  interest  is  probably  the 
critical  index  y  which  is  expected  to  be  independent  of  lattice  struc¬ 
ture.  Here  there  is  a  clear  difference  between  the  spin  one  half  and  the 
spin  one  exchange  models,  with  the  value  for  the  spin  one  Heisenberg 
model  (the  least  well  known)  lying  in  between.  However,  the  similarity 
in  values  of  y  for  the  spin  one  exchange  model  and  the  spin  one  half 
Ising  Model,  two  totally  different  models,  suggests  that  we  cannot  always 
distinguish  between  quantum  and  classical  models  on  the  basis  of  suscepti¬ 
bility  behavior  alone. 

As  far  as  real  spin  one  f erromagnets  are  concerned ,  we  expect  the 
inclusion  of  biquadratic  terms  in  their  interaction  to  lower  their  criti¬ 
cal  temperatures  and  reduce  the  value  of  the  critical  index  y  in  compari- 
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Table  6.4 

Comparison  of  critical  values  for  three  lattice  models. 


Critical 

Quantity 

Spin  j 

Ising 

Model 

Spin  j 
Heisenberg 
or  Exchange 
Model 

Spin  1 
Heisenberg 
Model 

kTc/J 

4.90 

4.07 

12.0 

Y 

5/4 

1.43 

1.38 

(s  -s  )/s 

oo  c  00 

.147 

.317 

.27 

Spin  1 
Exchange 
Model 


3.07 


1.2? 


.33 
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son  with  the  pure  spin  one  Heisenberg  ferromagnets .  It  seems  then  that 
the  details  of  the  intermolecular  interaction  can  produce  detectable 
differences  in  critical  behavior.  However,  the  main  question  which 
still  remains  unanswered  is  precisely  what  features  of  the  model  are 
responsible  for  all  the  differences  and  similarities  particularly  when 
comparisons  between  quantum  mechanical  and  classical  models  are  made. 


PART  II 


THE  FREQUENCY-DEPENDENT  INITIAL  PERPENDICULAR 

i 

SUSCEPTIBILITY  OF  AND  THE  INELASTIC  NEUTRON 


SCATTERING  FROM  THE  TWO  DIMENSIONAL  ISING  MODEL 


7.1 


CHAPTER  VII 

THE  PERPENDICULAR  SUSCEPTIBILITY  OF  THE  IS IMG  MODEL 


A.  Introduction 

The  Ising  model  is  one  of  the  simplest  and  most  studied  models 
which  exhibits  a  critical  phase  transition.  In  the  presence  of  a  mag¬ 
netic  field  H,  its  Hamiltonian  can  be  written  in  the  form, 

^  =  -J  2  <J±a .  -  mH  2  a.  ,  (7.1) 

nn  J  i 

where  a  is  a  two-valued  variable  (+1  or  -l)  and  m  is  the  magnetic  mo¬ 
ment  associated  with  it.  We  shall  be  primarily  concerned  with  the  in¬ 
terpretation  of  a  as  a  spin  variable  but  as  we  mention  in  chapter  one 
there  are  several  other  interpretations.  Clearly  (7.1)  is  a  much  sim¬ 
pler  system  than  we  considered  in  Part  I  if  only  for  the  fact  that  there 
is  no  question  of  non-commuting  variables  occurring  in  the  system.  Un¬ 
til  recent  years  much  of  the  work  on  the  Ising  model  had  been  concerned 
with  evaluating  its  equilibrium  properties  and  a  number  of  excellent 
reviews  exist  in  the  literature  (e.g.  Kadanoff  et  al.  1967;  Fisher  1965 
and  1967 ;  Domb  i960 ) . 

A  notable  contribution  to  the  study  of  the  time  dependent  fea¬ 
tures  of  the  model  was  made  by  Glauber  (1963)  who  made  use  of  the  5 
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•master  equation*  approach  which  discusses  the  first  order  differential 
equation  governing  the  time  evolution  of  the  probability  function  associ¬ 
ated  with  any  given  state  of  the  system.  His  theory  involves  the  import¬ 
ant  assumption  that  the  transition  probability  that  any  spin  will  jump 
from  one  its  states  to  the  other  depends  only  on  the  interaction  of  the 
given  spin  with  its  nearest  neighbours.  This  probability  is  expressed 
in  terms  of  two  parameters  which  simultaneously  describe  the  interaction 
of  the  Ising  system  with  an  external  heat  bath  and  the  interaction  be¬ 
tween  the  individual  spins.  It  is  determined  in  part  by  the  detailed 
balancing  condition  obeyed  by  the  equilibrium  state  of  the  model. 

Glauber  considers  the  problem  in  detail  for  the  case  of  a  closed 
N-member  chain  and  with  further  approximations  and  assumptions,  the  valid¬ 
ity  of  which  we  shall  not  discuss  here,  he  obtains  for  instance  a  param¬ 
eterized  expression  for  the  complex  frequency  dependent  magnetic  sus¬ 
ceptibility  defined  by  the  relation, 

M(t)  =X(w)H0e'1U)t  ,  (7.2) 

where  M(t)  is  the  magnetization  of  the  system  at  time  t  and  HQe"iu)t  is  a 
sinusoidally  time  dependent  applied  magnetic  field. 

Our  aim  in  the  following  calculation  is  also  to  obtain  an  expres¬ 
sion  for  the  frequency  dependent  susceptibility  of  the  Ising  model  but 
the  details  of  our  problem  and  the  techniques  we  use  are  quite  different. 
In  the  first  place,  we  notice  that  the  magnetic  moment  M  given  by, 
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M  s  m  Z  , 


(7.3) 


commutes  with  the  Ising  Hamiltonian  (7.1)  and  so  an  isolated  Ising  sys¬ 
tem  has  no  dynamic  properties.  In  Glauber's  case  it  is  the  interaction 
with  the  external  heat  bath  which  produces  the  time  dependent  features. 

However,  if  we  introduce  the  quantum  mechanical  Ising  Model  de¬ 
fined  by  the  Hamiltonian, 


(7.4) 


where  the  cr's  are  now  the  usual  Pauli  spin  matrices,  we  see  that  the 
magnetic  moment  operator  given  by, 


M  =  I,  (m  of  ,  m.of  ,  m  <jf)  , 

~  ^  x  k  ’  y  k  ’  zk 


(7.5) 


does  not  commute  with  (7.4)  and  so  is  not  a  constant  of  the  motion.  In 


particular,  if  we  specialize  to  the  case  of  a  perpendicular  field  H 


x 


only  and  perpendicular  magnetic  moment  M  ,  (7.4)  and  (7.5)  become 

A 


y  =% 


£m\  •  yj  •  «*»  1J11  lrn\  \J  .  to 

i  j  x  r  i  ’ 
nn  l 


(7.6) 


, ;  , 


. 
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and 


M 


x 


(7.7) 


and  we  have  the  model  whose  static  properties  were  first  investigated 
by  Fisher  (1963).  He  found  exactly  for  the  honeycomb  and  plane  square 
lattices  the  initial  perpendicular  susceptibility  defined  by, 


(7.8) 


where  M  is  the  magnetization  of  the  system  in  the  x  direction. 

A 

In  our  case,  we  are  concerned  with  the  dynamic  or  frequency  de¬ 
pendent  initial  perpendicular  susceptibility  defined  by  (7.8)  where  H 

X 

is  now  a  sinusoidally  time-dependent  magnetic  field  of  the  form, 


H  (t)  =  H  e~iwt 
x'  x 


(7.9) 


For  an  isolated  Ising  system,  we  shall  apply  the  Kubo  linear  response 
theory  (1957)  to  obtain  an  exact  expression  for  the  frequency  dependent 
perpendicular  susceptibility  without  the  introduction  of  additional 
parameters.  Detailed  consideration  is  given  to  the  honeycomb  and  plane 
square  lattices  and  comparison  with  Fisher’s  static  results  leads  to 
interesting  conclusions  regarding  the  susceptibility  of  an  isolated  sys¬ 
tem  as  opposed  to  that  of  a  system  in  thermal  contact  with  a  heat  bath. 
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B.  The  Rubo  Formalism 

The  Kubo  formalism  is  chosen  in  this  problem  because  it  is  a  par¬ 
ticularly  convenient  way  to  study  the  response  of  a  system  to  an  external 
mechanical  disturbance  such  as  the  perpendicular  magnetic  field  in  (7.6). 
As  first  pointed  out  by  Kubo  (1957)  it  can  be  applied  to  a  variety  of 
different  physical  problems  but  for  our  purposes  it  will  be  sufficient 
to  develop  the  general  formalism  in  the  language  of  magnetism. 

We  make  the  initial  assumption  that  our  magnetic  system  in  the 
remote  past  is  in  thermal  equilibrium  with  its  surroundings  without  the 
external  magnetic  field.  If  then  a  weak  time  dependent  field  H(t)  is 
introduced,  the  magnetization  M(t)  will  be  linear  in  the  field  and  tak¬ 
ing  into  account  possible  time  lag  effects  we  write 


__  t 

M  (t,T)  =  2  /  dt'jl  (t-t',T)H,(t>) 

M  ■»  _oo  ^  ’ 


(7.10) 


where  =  x.y.z  and  <jt  .(t,T)  is  known  as  the  response  function.  T  is 

yi  v 

the  temperature  and  t  and  t*  the  time.  It  is  important  to  note  that 
(7.10)  is  reasonable  only  for  weak  magnetic  fields  and  Kubo ’ s  linear 
response  theory  based  on  it  cannot  be  expected  to  yield  accurate  results 
with  large  fields.  However,  in  our  case  we  are  interested  in  the  zero 
field  limit  for  which  we  expect  (7.10)  to  be  asymptotically  exact.  The 
response  function  characterizes  then  the  temporal  behavior  of  our  system 
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and  plays  in  this  time  dependent  problem  much  the  same  role  as  the  par¬ 
tition  function  does  in  static  problems  such  as  in  Part  I.  As  with  the 


partition  function,  our  main  object  is  to  calculate  ^^(t,T)  since  it 


can  be  directly  related  to  the  experimentally  significant  macroscopic 
variables  of  the  system. 

In  particular,  to  obtain  an  expression  for  the  frequency  dependent 
susceptibility,  we  consider  the  response  of  the  system  to  a  sinusoidal 
field  of  the  form, 


H(t)  =  He-iut+6t 

^  rsj 


(7.11) 


Substitution  of  (7.11)  in  (7.10)  gives 


(7.12) 


=  L  He 


-i(to+i6)t  7ds^(s>T)eik+i6)s 


0 


(7.13) 


(7.1^) 


X,^(W,T)  =  7 dt^(t,T)ei(u+i6)t  . 


(7.15) 


,  .  n  q 
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The  infinitesimal  positive  constant  6  has  been  introduced  to  satisfy  our 
initial  condition  of  H(t  =  -  °°)  =  0.  We  always  finally  take  the  limit 
6—^  +  0.  (7.15)  indicates  that  the  frequency  dependent  susceptibility 

is  simply  the  one-sided  fourier  transform  of  the  response  function. 

To  obtain  ^  (t,T)  we  first  identify  the  magnetization  M(t,T)  of 

our  system  at  time  t  with  the  statistical  average  of  M,  the  magnetic 
moment  operator,  i.e. 


M(t,T)  =  trace[p(t,T)M]  ,  (7.16) 

with  M  given  by  (7*5) •  p(t,T),  the  time  dependent  density  matrix  of 

the  system,  is  obtained  by  solving  the  von  Neumann  equation  of  motion, 
namely , 

dpCt.T)  =  _i  (ff(t),p(t,T)]  ,  (7.17) 

dt  -ft  1 

where ^(t)  is  the  total  Hamiltonian  of  the  system  and  the  following 
initial  conditions  are  satisfied ; 


H(t  =  -  00)  -  0 


».T)  =  p 


eCi  trace(e~^°) 


(7.18) 


)^q  is  the  Hamiltonian  of  the  system  in  the  absence  of  the  external  field 


-j  .  - 
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H(t)  i.e.  we  write 


V(t)  =  \  +  \{t)  (7.19) 

with  Wj(t)  =  -  M.H(t)  .  (7.20) 

In  that  ease,  (7.17)  becomes 

=  -  i  [%,p(t,T)]  -  i  L^1(t).p(t.T)]  .  (7.21) 

and  since  we  are  interested  only  in  the  linear  response  of  the  magnetiza¬ 
tion  to  the  field,  we  replace  p(t,T)  in  the  second  term  on  the  right 
hand  side  of  (7.21)  by  p  given  in  ( 7 .18)  so  that 

*  ©C[ 

Mhli,.iL>0.p(t.T)]-i[Vl(t).peq]  .  (7.22) 

Hence  we  obtain  the  solution, 

t  )/*  ife,(t-t')/fi 

p(t,T)  =  pe  -i  /  dt'e  £^j_(t' )  »peqj  ©  .  (7.23) 

^  fi  — OO  ^ 


which  can  be  verified  directly  by  substitution  in  (7.22). 

Finally  to  obtain  an  expression  for  the  response  function,  we 
substitute  (7.23)  in  (7.16)  and  comparing  the  result  with  (7.10)  we  can 
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make  the  identification, 

^(t-f,I)  =  itraoe|e  QVpeqJ  e  '  (7.24) 

provided 

trace(peqM^)  =  0  .  (7.25) 

In  our  case  where  }i  =  x  and  is  given  by  (7.6),  (7.25)  is  certainly 
satisfied. 

(7.24)  is  the  important  relation  we  require  to  link  the  macroscop¬ 
ic  variables  of  our  system  to  the  microscopic  dynamics  as  provided  by 
the  Hamiltonian.  Using  the  invariance  property  of  the  trace  under  cyclic 
transformation,  namely, 

trace(ABC  ...  Z)  =  trace(ZAB  ...  Y)  etc.,  (7.26) 

and  setting  t1  =0,  (7.24)  can  be  more  compactly  written  in  the  form, 

,(t.T)  =  i  trace  £|>rpeq]  M^t)]  (7.27) 

=  i  trace  [tyeqy  t)  -  peqMvy  t)j  (7.28) 

=  |  trace{  Peq(1t(t)M^-MA(t))i  (7.29) 


' 
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_i 

-ft 


(7.30) 


In  (7.27  -  7.30)  we  have 


yt)  =  e  M^e 


(7.31) 


In  conclusion,  we  should  note  from  (7.30)  one  of  the  important 


features  of  Kubo's  linear  response  theory.  It  expresses  a  time-depend¬ 
ent  property  of  the  system,  namely  the  response  function,  as  a  statisti¬ 
cal  average  of  operators  for  the  equilibrium  system  i.e.  in  the  absence 
of  the  external  disturbance. 

C.  Application  of  the  Formalism 

To  obtain  the  frequency  dependent  initial  perpendicular  suscepti¬ 
bility  of  the  Ising  model  with  the  Hamiltonian  given  by  (7.6),  we  first 
must  evaluate  (7.30)  for  the  particular  case  of  p =  x  and  then  sub¬ 
stitute  in  (7.15). 


Firstly  we  note  that 


e 


(7.32) 


e 


(%  =  2J/-K) 
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-I  i 

=  He 
nn 


azaz 
i  j 


=  IT  (cos  -j-  + 
nn 


i0i°j 


sln  — ' 1 


where  we  have  used  the  identity, 


f 


for  n  any  positive  integer.  (7.3*0  is  completely  analogous  to 
miliar  van  der  Waerden  identity. 


/  z  z  N 

TF  (cosh  K  +  a. a.  sinh  K) 

nn  J 


t 


(K  =  j/kT) 


used  in  the  expansion  of  the  Ising  model  partition  function. 
Substituting  (7.3*0  in  (7.30)  for  p  =  =  x  and  using 


obtain 


im 

* 


(cos 


4  TT(1  -  iaV  tan  f  t) 
2  i',j'  \  nn  J  c 


(7.33) 

(7.34) 


(7.35) 

the  fa- 


(7.36) 


.7).  we 


nri 


-  •.  r.  L- 
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.  Zt  z 

]T  (i  +  ia  "tan 
nn 


(7.37) 


where  q  is  the  coordination  number  of  the  lattice  and  N  the  number  of 
lattice  sites.  Using  the  familiar  anticommutation  properties  of  the 
Pauli  spin  operators,  (7.3?)  becomes 


%-(cosf  t)1^  < 

*  2  i»,j» 


(i\) 


[  IT  (1-iCTCT  tan~t) 
s  nn  ^ 

i.  j  4  v 


(7.38) 


(1  +  iaVtanf  t)  TT  (1  -  io^.tanf  t)o£,  ,o*  T) 

k  —  1 


/Nc 

=  ^  Coos  if  t^cdd+t^t)  ■ 2 


(^-  q) 


IT  (1  -  tan2  ~t  -  2iaZaZ  tan^  t)  aX  oX  .  (7.39) 

t  2  k  i*  2  l 1  .1 1  / 


k  =  1 


Making  use  of  the  following  simple  trigonometric  identities, 


x  +  2COo  , 

a;  1  +  tan  — -  t  = 


2  0)o  + 

sec  —  t 

2 


5 


(7.40) 
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b)  u  =  tan  LOq t 


2  tan 


WO 

2 


t 


1  -  tan2  ^ 
2 


(7.41) 


c)  1  -  tan^  —  t 
2 


cos  OJ^t 


sec 


(7.42) 


we  obtain 


where  in  obtaining  (7.44)  we  have  made  use  of  the  commutator  identity, 

[AB.CJ  =  A[B,C]  +  [A,C]B  .  (7.45) 

The  commutator  in  (7.44)  can  obviously  be  non-zero  in  only  two 
cases;  a)  i *  =  J 1  —  j  and  b)  i*  and  j*  are  nearest  neighbour  sites. 
Considering  the  latter  case,  we  get 
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JT  t ,  d  -  ia^a^) [(1  -  ia* , a* ,u) , a*  ]  a *, 


(7.46) 


(7.4?) 


We  notice  in  (7.47)  that  oh,  is  the  only  variable  at  the  j'th  site. 

J 

Thus  when  the  thermal  average  is  taken  to  obtain  the  response  function 

i.e.  when  p  is  included,  only  terms  of  the  form  (o^,)na^t  (n  a  positive 

eq  J  J 

integer)  will  occur  at  the  j’th  site.  For  the  Pauli  spin  operators,  all 
such  terms  are  traceless.  The  thermal  averaging  involves  taking  the 
trace  over  the  direct  product  space  of  all  N  spins  which  will  be  zero  if 
the  trace  over  any  one  spin  is  zero.  Case  b)  then  gives  no  contribution 
and  (7.44)  can  now  be  written  as, 


^(t.T)  =  i~~(cosUJ0t)q  E  <^[  TT  (1-  ia*a  ji),a*]a* 

.1  k  —  1 


(7.48) 


Using  a  generalisation  of  (7.45)  we  then  obtain 


^(t.n 


2  q  /r  Q.  k  - 1 

=  iJL.(oo«Ul.t)q2<V  S  IT 

-n  v  j  =  l  i  - 1 


(1  -  ia.  cr.u)  * 


[(l-iaVu).ah  F  (1-iaVu)]  a*)>  (7.49) 

J i-k+1  1  J  J 


lo  eurc#i  ijlftb  ,Jb$bi/Iort±  el 

♦  cries  si:  rtxqa  erto  ynn  ievo  moB’iS  exi i 
.  /j  nr  ,r.  ••  «:  i  . 
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2im 

•ft 


(cosU^)t)q"1sinUJ0t  r,<^f  y~!  If  (l-iaZ<JZu)cy^ 


TT  (1  -  iaZaZu)|  aX^>  . 
i  J  J  CK 


(7.50) 


i  =  k  +  l 


We  now  expand  the  expression  within  the  j  ;  brackets  in  (7.50) 
in  powers  of  u.  Associated  with  u  is  a  sum  of  distinct  products  of  r  +  1 
pairs  of  spin  operators  of  which  r  are  z,z  pairings  and  one  is  of  the 
form  a*0^.  To  enumerate  all  the  products  we  first  of  all  consider  the 
q^r  +  1  c^s’^nc*'  selections  of  r  +  1  nearest  neighbour  sites  to  the  jth 
site.  Corresponding  to  each  selection  there  will  be  r+1  distinct  prod¬ 
ucts  of  the  form 


(a,!  orZ)(^Z  CTj)...(aZo^)...((j,  a\)  , 


kl  3 


k2  ^ 


(7.51) 


r+1 


z_y 


each  differing  only  in  the  location  of  the  pair.  Here  k]_k2,,,kr+l 

represents  a  particular  selection  of  r+1  nearest  neighbour  sites. 
(7.51)  can  be  rearranged  to  become. 


z  z  z  z  ,  z.a  y.  z. 

a  a,  ...a  ...a  (a.)  a. (a.) 

1  2  K  Kr+1  3  J  J 


3 


(7.52) 


where 


a+3+l=r+l 


(7.53) 
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If  we  examine  the  product  of  spin  variables  on  the  jth  site  in  (7.52) 

and  recall  the  anticommutation  properties  of  the  a's,  we  see  that  the 

r  +  1  distinct  products  will  cancel  in  pairs.  In  particular  if  r  +  1  is 

an  even  integer  there  will  be  no  resultant  contribution  to  the  coef- 

ficient  of  u  *  For  instance  for  the  simple  quadratic  lattice  for  which 

q  =  4  if  we  take  r  =  l,  then  the  two  products  associated  with  one  of  the 

six  possible  choices  of  two  sites  from  four  would  be  given  by  <j  ct  a.cr. 

2  J  J 

Z  Z  V  z 

and  a]_a2  jaj  w^ere  we  have  arbitrarily  labelled  the  two  nearest  neigh¬ 
bour  vertices  1  and  2.  These  two  terms  anticommute  and  their  sum  is 
zero.  Similarly  the  other  five  pairs  of  terms  cancel  so  that  the  coef¬ 
ficient  of  u  for  the  simple  quadratic  lattice  is  zero. 

If,  however,  r  +  1  is  an  odd  integer  then  there  will  always  be 
one  spin  product  uncancelled  which  we  can  chose  to  be  of  the  form, 


z  z  z 

ak  ak  •••ak 
1  2  r+1 


°^  =  °k  ak  •”ak  °j 
J  J  ^  k2  Kr+1  j 


(7.54) 


using  (7.35).  For  example  for  the  simple  quadratic  lattice,  there  are 
four  spin  products  contributing  to  the  coefficient  of  u^  in  (7.50) » 


namely, 


zzzy  zzzy  zzzy  zzzy 
°i°2a3aj  ’  ala2a4aj  *  al°3a4aj  ’  °2a3a4a,j 


(7.55) 


one  each  for  the  four  possible  selections  of  three  nearest  neighbour 
sites  to  the  jth  (labelled  1,2, 3, 4). 


v;  2  ‘o  ©rid*  10  aoiJ  otq  a  ten  .  •  -  ©,••••'•  ilaoo'x  b..n 
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From  all  these  considerations  we  can  now  write  (7.50)  in  the  form, 


^(t.T)  =  SiSU(cOSUJ0t)q-1sinOJQt(  £ 

2r +  l£q 


(7.56) 


2 


(cosL0nt)q"1sin(jOnt(  y~!  (-l)ru2r> 

*  u  u  r  =  0 

2r  +  1  <  q 


q  ?-^+  V  z  z  z  z 

2_.  <ak  ak  •••°k  a0 

k  x  K1  k2  K2r  + 1  u 


» 


(7.57) 


where  in  (7.56)  and  (7.57)  we  have  invoked  the  translational  symmetry  of 
the  lattice  to  replace  the  summation  over  j  by  a  factor  of  N  and  re¬ 
labelled  the  jth  site  as  the  zeroth.  The  first  summation  in  (7.56) 
and  (7.57)  is  over  all  positive  integer  values  of  r  such  that  2r  +  l^q. 
The  second,  summation  is  over  all  possible  selections  of  2r  + 1  sites 
from  the  q  nearest  neighbours  to  the  zeroth. 

The  equation  (7.57)  is  the  desired  exact  relation  between  the 
response  function  and  the  equilibrium  spin  correlation  functions,  valid 
for  all  regular  lattices.  The  form  of  the  relationship  depends  only  on 
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the  coordination  number,  not  on  the  particular  lattice  nor  even  on  the 
dimensionality  of  the  lattice.  For  two  dimensional  lattices  all  the 
required  correlation  functions  are  either  known  or  can  be  evaluated  by 
known  methods  (Montroll,  Potts  and  Ward  1963).  In  the  next  chapter,  the 
particular  form  of  the  response  function  will  be  determined  for  the 
honeycomb  and  square  lattices  and  used  in  conjunction  with  (7. 15)  to 
obtain  the  frequency  dependent  initial  perpendicular  susceptibility  of 
these  two  lattices. 


..t  „  b  tc  •  r.n~.  b 
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CHAPTER  VIII 

APPLICATION  TO  HONEYCOMB  AND  PLANE  SQUARE  LATTICES 


A.  Honeycomb  Lattice 

Putting  q  =  3  in  (7.57)  the  response  function  for  the  honeycomb 
lattice  is  given  by, 

(8.1) 

2 

^(t.T)  =  cos2(f>0t  sintLQt(3  ^aoal^~  » 


where  spins  1,  2,  and  3  are  the  nearest  neighbours  to  spin  0.  Tak¬ 
ing  the  one-sided  fourier  transform  of  ^JL(t,T)  according  to  (7.15)*  we 
obtain  the  frequency  dependent  initial  perpendicular  susceptibility 
X^(W,T)  in  the  form, 


=  5  c  |  <a0°x>  +  i  <a0ala2°3>^  — 


(w/3 


+  l  t<°0a3>  -  <®0<W*3/3  LTV 


(Uj/(jl30) 


2  * 


(8.2) 


where  a  =  aZ.  For  the  honeycomb  lattice,  (8.1)  and  (8.2)  provide  exact 
solutions  when  the  known  values  of  the  correlations  are  inserted. 

Any  correlation  of  a  finite  number  of  spins  in  a  two  dimensional 
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Ising  model  can  be  obtained  by  the  dimer-Pfaffian  method  (Mantroll, 

Potts  and  Ward  1963)  which  we  will  not  discuss  here.  However,  the  near¬ 
est  neighbour  correlation  <(jJo(Jl/>  course  directly  proportional  to 

the  energy  and  was  first  calculated  for  the  honeycomb  lattice  by 
Houtappel  (1950)  to  give 


<?0°j>=  |  l>th2K  +  aX(kx)]  . 


(8.3) 


in  which  }£(k^)  is  the  complete  elliptic  integral  of  the  first  kind  de¬ 
fined,  for  example,  by  Jahnke  and  Emde  (1945)  as, 


TT/2 


f 


sin  x 


(8.4) 


0 


Also 


K  >  K 


K  <  K 


c 


c 


(8.5) 
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and 


cosh  2K(cosh  2K  -  2) 
4  TT  sinh4  K 


K  >  K 


a  = 


(8.6) 


sinh  4K(cosh  2K-2)(cosh  K) 


4  TT  sinh2  K(cosh  3K)1^2 


1/2) 


K  <  K 


where,  again,  K  =  j/kT. 

The  four  particle  correlation  is  most  conveniently 

calculated  by  using  a  transformation  theorem  developed  by  Fisher  (1959). 
This  theorem  shows  that  a  correlation  function  involving  a  spin  can 
be  expressed  as  a  linear  combination  of  correlation  functions  involving 
the  q  spins  a^,  ...  which  are  nearest  neighbours  to  <Jq.  In 

our  proof,  we  shall  assume  initially  that  the  interaction  constant  J  may 
not  necessarily  be  the  same  for  all  bonds  in  the  lattice.  In  that  case, 
the  correlation  product  between  the  spins  C7q,  <j^,  0^,  ...»  o^,  is  given 


by. 


<£oag°k  *•*  Qm)>= 


trace [a0aak 


•  •  QL® 
m 


(f-j 


] 


{  £  Kijaiaj} 

trace  e 


(8.7) 


where  the  trace  is  over  the  direct  product  space  of  all  the  spins  in  the 
lattice.  Let  us  consider  in  particular  the  trace  over  the  zeroth  spin 
space  in  the  numerator  of  (8.7)  i.e. 


■ 
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To(al,a2’***,aq)  = 


(Kla0al  +  K2°0c,2  + 


+  VoV 


]  (8.8) 


where  we  now  regard  a^,  (Jg ,  o^*  •  •  • » aq  as  scalar  parameters  which  can 
equal  to  -hi  or  -1,  We  may  try  to  represent  (8.8)  in  the  form. 


T  0  ( *  a2  ’  •  •  *  *  °q )  —  j 


C 

q  q  3 

E  “ J 7  +  E  Ma  °b  °c 

r  =  l  t  =  1  tt  t 


(8.9) 


qS  ->  r  (  K1°0°1  +  Vo°2  +  •  •  •  Kq°Oaq)  q 

+  ErA  °b  +...°  Uraoe^  J  , 

U=1  u  u  u) 


where  the  q  coefficients  ar,  the  coefficients  0^,  the  coef¬ 
ficients  yu*  etc.,  are  functions  only  of  the  q  interaction  parameters 
Kq,  Kg,  ...»  K^,  between  the  spin  aQ  and  its  nearest  neighbours 

a>,  ...»  a  respectively.  The  summations  in  (8,9)  are  successively 
d.  q 

over  all  distinct  selections  of  1,  3>  3*  etc.  spins  from  the  q  nearest 
neighbours  to  aQ.  Because  (8.8)  changes  sign  when  the  signs  of  the 
spins  Uq,  Og,  ...»  0^  are  simultaneously  altered,  we  do  not  need  to 
consider  combinations  of  even  numbers  of  spins  in  (8.9).  The  number  of 
unknown  parameters  in  (8.9)  is 


o  ® 


=  2 


qci  +  qc3  +  qc5  + 


(8.10) 
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There  are  2q  distinct  joint  spin  states  of  the  nearest  neighbour  spins 
and.  the  identity  of  (8,8)  and  (8,9)  must  hold  for  all  of  them.  However, 
since  both  expressions  change  sign  on  simultaneous  reversal  of  the  q 
nearest  neighbour  spins,  we  obtain  only  (l/2)2^  =  2^“^  independent  equa¬ 
tions  between  (8,8)  and  (8.9)  which  therefore  will  exactly  determine  the 
set  of  coefficients  a  f  |3^,  yu  •  •  •  ensuring  that  (8.9)  is  a  faithful 
representation  of  (8.8).  On  inserting  this  transformation  into  (8,7)* 
we  get 

(8.11) 

^<\Vab+0otagah-"am>+ 

U  U  Tj 


'aAa  a 

\  0  g  h 


q  q  3 

. . .  am>  =  £  ar<ara  ah . . . +  77 
r=l 


which  expresses  the  correlation  involving  <Jq  as  a  sum  of  correlation 
functions  including  the  nearest  neighbour  spins  Ug*  •••  Oq  in  place 


We  shall  use  this  theorem  to  determine  in  terms  of 

the  second  neighbour  correlation  ^or  h°neycortlk  lattice 

with  a  unique  interaction  constant.  In  that  case  (8.8)  becomes 


K(<JqCFi  +  cJqC?2  + 

T0^al,a2,a3^  s  trace0[a0e  J  * 


and  from  the  symmetry  of  the  lattice,  (8.9)  is 

C  1  K  (a0a1+a0a2+a0a3) 

Tq(  al »  a2 * a3 )  =  f  a(  ai+  a2+  Q3 ^  +  ^  pla2a3 M  traced-e 


(8.12) 


(8.13) 

]. 


Wow  there  are  only  two  distinct  joint  states  of  the  three  nearest  neigh¬ 
bour  spins  corresponding  to  a)  o1  =  a,  =  =  +1  and  b)  typically 

a]_  =  -1»  a2  =  ^3  =  whi°h  will  be,  however,  sufficient  to  determine  the 
two  unknown  parameters  in  (8,13).  The  identity  of  (8.12)  and  (8.13)  for 
case  a)  leads  to  the  equation, 

3ot  +  3  =  tanh  JK  ,  (8.14) 


and  for  case  b)  to 


0C  -  |d  =  tanh  K  .  (8.15) 

Solving  (8.14)  and  (8.15),  substituting  (8,13)  in  (8.7)  and  making  use  of 
(7*35)  together  with  the  symmetry  of  the  lattice,  we  finally  obtain 

(8.16) 

\a0ala2a3/  =  \  ( tanh  3K  + tanh  ^aiCT2^  + 1  ( tanh  3K  -  3  tanh  K)  • 

Using  an  exactly  similar  procedure  for  <^aoal/  »  we  also  get 

^1^  =  ^cothK  +  3  tanh  K)  <^aQa^)>  -  i  (1  +  tanh2  K)  .  (8.17) 

Substitution  of  (8.17)  in  (8.16)  gives  after  some  trigonometric  simpli¬ 
fication, 

<^0ala2a3^  =  tanh2  K)  "  2  tanh  K^tarih2  K  +  3).  (8.18) 
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(8.1d)  represents  the  exact  known  value  of  the  four  particle  correlation 
since  we  have  already  stated  the  exact  solution  of  <\aQal/>  equations 
(8.3)  to  (8.6). 

The  most  obvious  check  on  X^(U),T)  given  in  (8.2)  is  to  compare 
the  zero  frequency  limit  with  Fisher's  (1963)  static  susceptibility  re¬ 
sult  for  the  honeycomb  lattice.  Substituting  (8.18)  in  (8.2)  we  find 
for  LO  =  0  that 

jxV.t) 

- - -  =  —  tanh  K(tanhi  K  +  3)  -  tanh^K  Fana-,\  .  (8.19) 

Nirr  3  x  u  -1-/ 

Specializing  Fisher's  equation  (4.17)  to  the  honeycomb  lattice  and  mak¬ 
ing  use  of  (8.17),  we  find  perfect  agreement  between  ours  results  and 
Fisher's. 

n 

The  most  striking  feature  ofXj_(tO,T)  consists  of  the  two  reso¬ 
nances  occurring  at  Go  =  2J/h  and  6j/h.  It  is  easy  to  see  the  physical 
origin  of  these  resonances.  In  the  honeycomb  lattice,  a  given  spin  will 
be  found  in  one  of  four  environments.  In  two  environments,  all  three  of 
its  nearest  neighbours  are  oreintated  either  parallel  or  antiparallel  to 
the  central  spin.  In  these  configurations,  the  energy  required  to  flip 
the  central  spin  will  be  respectively  -»-6j  and  -6J  and  the  00  =  6j/K  reso¬ 
nance  corresponds  to  an  absorption  of  energy  from  the  magnetic  field  by 
those  spins  in  the  all  parallel  state.  In  the  other  two  environments, 
the  central  spin  has  either  two  of  its  neighbours  parallel  and  one  anti- 
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parallel  or  vice  versa.  Clearly  the  QJ  =  2j/<tt  resonance  occurs  because 
of  spins  in  the  two  parallel,  one  antiparallel  configuration.  It  is 
interesting  to  note  that  for  critical  temperatures  of  the  order  of  a 
few  degrees  Kelvin  these  resonances  are  in  the  microwave  region. 


It  should  be  pointed  out  that  theX^_(tO,T)  which  we  have  been  cal¬ 


ling  the  perpendicular  susceptibility  is  in  fact  the  real  part  of  a  com¬ 
plex  frequency  dependant  susceptibility  of  which  the  imaginary  part  is 
given  by, 


r<P>)(&(3GL^-u>)  -  6(3u^4U))) 


and  corresponds  to  the  energy  absorption  by  the  spin  system.  The  reso¬ 
nances  in  (8,20)  are  of  zero  linewidth  because  the  only  interactions  in 
this  model  are  the  nearest  neighbour  Ising  interactions.  Introducing 
additional  near  neighbour  interactions  would  simply  increase  the  number 
of  zero  linewidth  resonances.  In  a  real  crystal,  of  course,  additional 
long  range  forces  such  as  magnetic  dipole-dipole  interactions  and  spin- 
lattice  forces  will  result  in  a  finite  linewidth  to  the  resonances. 
However,  if  a  crystal  approximates  an  Ising  lattice  at  least  the  loca¬ 


tion  of  the  resonances  should  be  predicted  by  this  analysis. 

Returning  to  consideration  of  (8.2),  the  real  part  of  the  fre- 


. 
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quency-dependent  susceptibility,  it  is  interesting  to  examine  the  temper¬ 
ature  dependence  of  the  amplitude  of  the  resonances.  At  T  =  0,  where 
all  even  correlations  are  unity,  corresponding  to  the  completely  mag¬ 
netized  state,  the  resonance  atUG  =  U)q  in  (8.2)  disappears.  This  is  be¬ 
cause  with  all  spins  aligned  the  energy  necessary  to  overturn  a  single 


spin  can  only  be  +6J.  At  the  other 


»  1,  and  the  amplitudes  of  the  U)  =(j0q  and  Lo  =  3^  resonances  are  pre¬ 
dicted  by  (8.2)  to  be  in  the  ratio  of  3  to  1.  This  simply  corresponds  to 
the  relative  probabilities  in  a  random  array  of  spins  for  the  arrangement 
of  two  first  neighbour  spins  parallel  and  one  antiparallel  to  the  arrange- 
nebt  of  all  spins  parallel. 

We  shall  also  consider  the  critical  region  given  by  sinh^K^  =  l/2 
for  the  honeycomb  lattice  (Houtappel  1950)*  From  (8.5)  K^K  corresponds 
to  1  and  ar'-'0  and,  from  Jahnke  and  Emde  (19^5)*  to 


(8.21) 


1 


Expanding  these  terms  in  a  Taylor  series  about  K=KC  and  retaining  only 
lowest  order  terms,  we  find 


(8.22) 


c 


c 


and  substituting  (8.22)  and  (8.18)  in  (8.2)  that 


j  .  ixbnoc  ••  ■  j  jj  • 

•  . 

£ 

&icqijn&  ©no  '  le  I  q  .  or  •  :er  :•  •  1  o  i  :o 
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-4d£(w.t) 

Nkt 


(8.23) 


±2  k“  L _ i _ 

TT  c  l-(U)/u0)2 


-111 _ ]U- 

i-(io/3<^)2 


1- 


We  see  that  for  all  frequences  X  f(w»T)  is  finite  in  magnitude  but  has 

an  infinite  slope  at  T  .  In  particular  at  low  frequencies  the  slope 

c 

is  positive  infinite.  This  behavior  contrasts  sharply  with  the  much 
stronger  singularity  exhibited  at  T  by  the  susceptibility  in  the  spin 
one  exchange  model.  This  is  very  much  related  to  the  fact  that  the 
latter  susceptibility  can  be  shown  to  be  equivalent  to  an  infinite 
sum  of  spin  correlation  functions  whereas  from  (8.2)  and  (8.18)  the 
perpendicular  susceptibility  of  the  Ising  model  is  related  to  only  one 
correlation. 

To  illustrate  the  frequency  and  temperature  dependence  of  the 

perpendicular  susceptibility  of  the  honeycomb  lattice,  the  behavior  of 

X[(OJ,T)  as  a  function  of  temperature  for  various  fixed  values  of  UJ 

approaching  00  =OOq  is  given  in  figure  8.1.  In  figure  8.2/X/^(w*T)  as 

a  function  of  UJ  for  various  values  of  T  is  shown.  The  numerous  features 

of’X^(ljO,T)  mentioned  above  are  manifest  in  these  illustrations. 

*■1— . 

B.  The  Plane  Square  Lattice 

An  expression  for  the  frequency  dependent  initial  perpendicular 
susceptibility  of  the  plane  square  lattice  can  be  found  in  just  the  same 


. 

\  xiaqieq 


Fig  8.1. 


The  temperature  dependence  of  the  initial  perpendicular 
susceptibility  of  the  quantum  mechanical  Ising  model  on 
the  honeycomb  lattice  for 

a)  LO  =  0  b)  03  =  0.6cOq  and  c)  10  =  0.8  03^. 
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Fig  8.2 

The  frequency  dependence  of  the  initial  perpendicular 
susceptibility  of  the  quantum  mechanical  Ising  model  on 


the  honeycomb 


lattice  for 

=  T  and 
c 


T  ->  oo  . 


T 


0,  T 


i(«,T)/xU<O.T) 
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way  as  for  the  honeycomb  lattice.  For  any  regular  q  =  4  lattice  (plane 
square,  Kagome,  diagonal),  the  response  function  expression  (7.57)  be¬ 
comes 


^(t.T)  -  C0S3U)J  sinco0t[ 


-  tan  coQt 


aOalaZa% 


(8.24) 


where  now  spins  1,  2,  3  are  three  of  the  four  nearest  neighbours  to  spin 
zero.  Again  from  (7.15)  we  obtain 


(w/2co0)- 


(8.25) 


For  the  plane  square  lattice  the  nearest  neighbour  pair  correla¬ 
tion  was  first  evaluated  by  Onsager  (1944)  and  is  given  in  convenient 
form  by  Fisher  (1963)  as> 

<a0a^>=  |  coth  2K[l+~  (2  tanh2  2K-l)]^(k1)]  ,  (8.26) 

where  now  ?  k-^  =  2  tanh  2K/cosh  2K  ,  (8,27) 

for  all  temperatures.  Again  the  four  particle  correlation  may  be  calcu- 
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lated  by  the  transformation  theorem  described  in  the  last  section.  By 
an  exactly  similar  procedure  we  obtain  first 


a0ala2a3>=  £  ( tanh  4K  -  2  tanh  2K)  +~  tanh 


(8.28) 


+ i  tanh  4K<(a1<x^>  +  ^  ( tanh  4K  +  2  tanh  ZK)<(o^<J2a^0l\^> , 


and  also, 


<^la2a3a4/>  =  1  +  ^  coth2  2K  -  4  coth  2K(1  +  coth2  2K)  * 

<a0a-j^>  +  2  coth2  2K[2<(a1a2)>  +  <^ia3/0  * 


(8.29) 


where  and  <\^Jla3/>  rePresent  the  second  and  third  neighbour  Ising 

correlations  respectively.  Then  substituting  (8.29)  in  (8.28)  we  get 


4oW3>  “  COth  2KL2<^ia2  ^  +  <^ala3^)+1J  -  (i  +  2  coth2  2K)<^r0a:i^>  . 

(8.30) 

The  second  and  third  neighbour  correlations,  first  worked  out  by  Kaufman 
and  Onsager  (1949),  are  also  given  in  a  convenient  form  by  Fisher  (1963) 
which  is  valid  for  all  temperatures  i.e. 


<0i023=  i  o°th2  2K[E(ki)  +  k^(ki)]  , 


(8.31) 
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and 


?103>  =iooth2  2K-<fl|-)2LE2(kl)  -  2k£)£(k1)E(k1)  +  (k^)3X2(k1)H  ■ 

(8.32) 


where  k-^  is  given  in  (8.2?)  and 


k"  =  2  tanh2  2K  -  1  . 


(8.33) 


E(k-^)  is  the  complete  elliptic  integral  of  the  second  kind  defined  by 
Jahnke  and  Emde  (1945)  as, 


E(kx)  = 


sin^x  dx 


0 


(8.34) 


Equations  (8.25)  to  (8,34)  constitute  an  exact  known  solution  to 
the  perpendicular  susceptibility  of  the  plane  square  lattice.  In  many 
respects  it  is  qualitatively  very  similar  to  that  of  the  honeycomb  lat¬ 
tice.  There  are  resonances  at  U]  =  4J/ft  and  8j/h,  just  as  we  might  ex¬ 
pect,  corresponding  respectively  to  configurations  in  which  three  neigh¬ 
bours  to  a  given  spin  are  parallel  and  one  antiparallel,  and  all  four 
spins  parallel.  (We  shall  turn  later  to  the  case  of  two  spins  pointing 
in  each  direction.)  Again  the  resonances  have  zero  linewidth. 

At  T  =  0,  the  resonance  at  2U^  disappears  as  expected,  and  only 
the  resonance  at  4(jl^  corresponding  to  all  spins  parallel  persists.  As 


1 
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T— >  oo  the  amplitudes  of  the  2u)q  and  4o)q  resonances  approach  the  expected 
ratio  of  4  to  1.  Each  of  the  above  remarks  applies,  of  course,  to  any 
q  =  4  lattice. 

Wear  the  critical  temperature,  given  for  the  plane  square  lattice 
2  s 

by  sinh  2Kc  =  1,  we  again  perform  Taylor  series  expansions  about  Tc  and 
obtain  from  (8.26),  (8. 31)  and  (8.32) 


<°Ocl>=^-^K-Kc)los|K-Kc 


(8.35) 


and 


K*)log|K-K= 


(8.36) 


(K-Kplog 


k-  k: 


(8.37) 


Substituting  in  (8,30)  gives 


<a0ala2°3>  =  V2 


4(TT-1)  ...  1 

-  JL  (8  -TT)(K  -  KS)log 

K-  Ks 

TT  2  2_ 

IT  2  0 

c 

(8.38) 


f 


and  hence  for  the  analogue  of  (8.23)  we  get 


-i-i;x>,T)  -■*>,  t®)]  = 


(8.39) 


2K5 
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4  -  IT  1 
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We  see  that  again  the  perpendicular  susceptibility  is  finite  at  T  ^  and  has 
infinite  slope.  However,  there  is  now  the  important  difference  that  for 
low  frequencies,  the  slope  is  negative  infinite. 

This  provides  only  a  premonition  of  the  discrepancies  to  follow. 
If  we  now  put  CO  =  0  in  (8.25)  and  substitute  (8. 30),  we  obtain 


8JX^(0,T) 


(8.40) 


-  =  2(4  +  3  coth2  2K)<a0a1^>  -  3  coth  2Kj_2<^y1a^>+<^J1a^>+  lj  , 


whereas  Fisher's  static  perpendicular  susceptibility  result  is 


8jXx(T)  f  2  2  /  \ 

- - —  =  K<  6(1  +  cotn  2K)  -  12  coth  2K(l  +  coth*  2K)/<Jna,  > 

Nm2  V  \  0  (8.41) 


+2(3  coth2  2K 


8jX3_(o,T) 

Nm2 


where  in  (6.41) X1(0,T)  is  given  by  (8.40).  In  other  words  Fisher's  static 
result  for  the  square  lattice  consists  of  our  result  for  zero  frequency 
plus  an  extra  term  containing  a  factor  K.  This  is  in  striking  contrast 
to  the  case  for  the  honeycomb  lattice  where  the  zero  frequency  result  and 
Fisher's  static  result  agree  exactly.  The  resolution  of  this  discrepancy 
will  be  expounded  both  physically  and  mathematically  in  the  next  chapter. 
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CHAPTER  IX 

THE  ZERO  FREQUENCY  LIMIT 


We  must  now  tackle  the  question  of  why  our  zero  frequency  re¬ 
sults  agree  with  Fisher’s  static  results  for  the  honeycomb  lattice  and 
not  for  the  plane  square  lattice.  More  generally  it  seems  that  a  static 
calculation  of  the  type  carried  out  by  Fisher  (1963)  would  agree  with 
the  U)~K)  dynamic  result  for  any  regular  lattice  of  odd  q  and  disagree 
for  any  regular  lattice  of  even  q.  Indeed  an  examination  of  his  equa¬ 
tion  (4.10)  reveals  how  terms  involving  a  factor  K  can  arise  for  even 
q  and  not  odd,  whereas  the  Kubo  formalism  as  applied  in  previous  chap¬ 
ters  seems  unable  to  yield  such  factors  for  any  lattice. 

First  we  must  realize,  as  Kubo  (1957)  himself  emphasizes,  that 
we  have  just  calculated  the  susceptibility  of  an  isolated  system  on 
which  the  magnetic  field  H  has  been  imposed  adiabatically  (because  of 
the  e6i:  factor).  On  the  other  hand,  Fisher's  calculation  is  of  the 
isothermal  susceptibility  of  a  system  in  thermal  contact  with  a  heat 
reservoir.  Thus  they  need  not  be  the  same. 

To  see  how  this  distinction  might  be  important  for  even  q  lat¬ 
tices  and  not  for  odd  q,  we  note  that  the  former  will  contain  at  all 
non-zero  temperatures  configurations  in  which  a  central  spin  has  half 
of  its  neighbours  aligned  parallel  to  it  and  the  other  half  aligned 
antiparallel.  Such  a  ’sublattice’  of  spins  are  effectively  free  spins. 
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It  is  well  known  that  for  a  system  of  free  spins  there  exists  a  great 
difference  between  the  adiabatic  and  isothermal  susceptibilities.  In 
the  former  case,  a  free  spin  merely  executes  a  Larmor  precession  about 
the  direction  of  the  applied  field  with  no  change  in  the  component  of 
the  magnetic  moment  of  the  spin  along  the  direction  of  the  field.  The 
adiabatic  susceptibility  is  simply  zero  and  it  is  no  surprise  then  that 
our  result  (8.25)  Tor  the  square  lattice  contains  no  contribution  associ¬ 
ated  with  such  free  spins.  On  the  other  hand,  in  a  calculation  of  the 
isothermal  susceptibility,  it  is  always  implicitly  assumed  that  there 
exist  weak  interactions  with  the  heat  bath  that  are  important  in  provid¬ 
ing  a  mechanism  by  means  of  which  the  system  can  approach  its  equilib¬ 
rium  state  in  the  presence  of  the  magnetic  field,  but  that  can  be  neg¬ 
lected  in  determining  the  nature  of  that  equilibrium  state.  The  isother¬ 
mal  susceptibility  of  a  system  of  free  spins  is  non-zero  because  of  this 
assumption. 


Our  conclusion,  then,  is  that  the  termXQ(T)  given  by, 


=  K  )  6(1  +  coth2  2K)  -  12  ooth  2K(1  +  coth2  2K)<a0oJ> 


(9.1) 


which  represents  the  difference  between  Fisher* s  result  and  ours  can  be 
considered  as  the  contribution  of  the  sublattice  of  free  spins  to  the 
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isothermal  initial  perpendicular  susceptibility  of  the  Ising  model  on 
the  plane  square  lattice.  The  contribution  (9.1)  appears  because  of  the 
weak  interactions  discussed  above  and,  as  expected,  does  not  appear  in  a 
calculation  of  the  adiabatic  perpendicular  susceptibility,  where  such 
interactions  are  explicitly  omitted. 

To  illustrate  further  the  validity  of  this  assertion,  we  shall 
introduce  a  mechanism  by  means  of  which  it  is  possible  to  calculate 
Xq(T)  by  the  Kabo  formalism.  We  now  consider  the  problem  where ^  is 
given  not  by  (7.6)  but  by 


$5  =  J  2  aV  -  mH z  2  o?  ,  (9.2) 

nn  J  i 

i.e,  an  external  static  magnetic  field  is  introduced  in  the  z  direction. 
The  perpendicular  susceptibility  can  also  be  calculated  with  (9*2)  as 
the  unperturbed  Hamiltonian.  We  can  see  that  the  formerly  free  spins 
can  now  be  expected  to  give  a  contribution  to  the  susceptibility  as  de¬ 
termined  by  the  Kubo  formalism  and  in  particular,  if  theio->0  limit  is 
taken  first  and  then  the  H  — >  0  limit,  (essentially  the  reverse  of  what 
we  did.  in  chapter  eight)  we  should  be  able  to  make  a  direct  comparison 
with  Fisher's  result.  By  so  ordering  our  limits  we  are  doing  precisely 
what  is  required,  namely,  considering  an  interaction  which  is  important 
in  guiding  the  system  to  its  equilibrium  state  but  unimportant  in  de¬ 
termining  its  details  (such  details  depending  primarily  on  the  Ising 
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interaction. ) 

The  calculation  is  very  similar  to  that  outlined  in  previous 
chapters  and  we  shall  only  emphasize  the  important  details  here.  Cor¬ 
responding  to  the  identity  (7. 34),  we  make  use  of  the  relationship, 

+  it%  (9.3) 

e’  _  TFLcos(-^)  +  iazaz  sin(^~)]  H[cos(-^)  +  ia^  sin(~^)]  , 

nn  2  J  j  J 

where  7^  is  given  by  (9.2),  6U)=  2mH ^/fi  and  the  second  product  is  over 
all  lattice  sites.  Corresponding  to  (7.38) »  we  get 


^xx(t,T,Hz)  =  ig-(cos  OJ0t)q  cos 


6wt 


(9.4) 


(1’iakai»u),<Jv^(1”iai'ui)ai' 


q 

+  TT 

k  =  1 


/  z  z 
(1  -  ia,  a. 
k  i 


u)[(l 


iaz  u  )  cr 
i»  ly  i> 


where  u^  =  tan  6tot.  For  the  plane  square  lattice,  using  the  same  pro 
cedures  as  before,  we  then  obtain 


8Nm 


cos 


(6out) 


cos^aJQt  sin  uo0t  <(^oaly>H 


9.5 


-  cos  CuQt  sin^u^t  <^oal°2a3/H 


y  +  Mb...  sin  6cut  * 
h 


(oos  V  ^0>H  -  sinVKW^H  + 


+  ,inV45i?3w 


(9.5) 


where  the  subscript  H  indicates  the  presence  of  and,  as  before,  a=<J. 
The  term  in  (9.5)  proportional  to  cos(6ujt)  corresponds  to  (8,24).  When 
the  fourier  transform  is  taken  to  obtain  its  contribution  to  the  per¬ 
pendicular  susceptibility,  the  zero  frequency  and  zero  field  result 
(8.40)  is  obtained  no  matter  in  which  order  uu-tO  and  Hz~>  0  limits  are 
taken. 

Simple  trigonometry  shows  that  the  remainder  of  (9.5)  separates 
into  parts  proportional  to  sin  6oat,  sin  6 out  cos  2(jjQt  and  sin  6 cot  cos  4u^t . 
The  first  is  given  by 


si„  tatl)^ .  +  2(g0ol(i^)  +  3^^ ^Vk  ]  •  (9-6) 

It  can  easily  be  demonstrated  that  the  other  two  terms  give  no  contribu¬ 
tion  to  the  zero  field,  zero  frequency  susceptibility.  However  the  term 
(9.6)  requires  careful  attention.  Using  (7.15)*  this  gives  the  following 
contribution  to  the  perpendicular  susceptibility, 
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%®(Ou,T,Hz) 


Nm  6  (jo 
445  (6U/)2-U ? 


(9.7) 


(440°1o2>H  +  2^0°la3^)  +  34ala233a4^]  • 


We  see  that  this  term  can  be  associated  with  the  formerly  free  spins 
where  now  an  energy  -ft6uj is  required  to  overturn  one  of  them. 

We  now  apply  Fisher's  transformation  theorem  with,  however,  a 
number  of  important  differences  from  the  procedure  previously  described. 
Corresponding  to  (8.8)  we  now  have 


T0(a. 


1»°2,w3 


V 


=  traceQ^a^e 


«°0al +  °0a2  +  aoa3  +  °oV  +aao-]  (9>8) 


where  a  =  mH  /kT.  However  because  of  the  aaA  term  in  (9.8)*  it  is  212, 

Z  U 

longer  antisymmetric  when  all  four  nearest  neighbour  spins  are  simul¬ 
taneously  reversed.  We  must  therefore  write  (8.9)  in  the  form, 


V’i'WV  = 


4  6  4 

£  +a  Z  °r  +  *  Z  aa  °b  +Y  2  0 

r=l  t=l  t  t 


a  a 
_  a  b  c 
u  =  1  u  u  u 


(9.9) 


6<Jl°2a3CT4 


traceQe 


K(  aQOl  +  a0a2  +  +  aQa4)  +  aaQ 


To  determine  oc,  |3,  Y*  6  and  £.  *  we  obtain  five  independent  equations  by 
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identifying  (9.8)  and  (9.9)  when  a)  all  four  nearest  neighbour  spins 
are  +1,  b)  all  four  -1,  c)  three  -i-l  and  one  -1,  d)  three  -1  and  one 
+1,  and  e)  two  +1  and  two  -1.  Solving  these  equations  and  then  sub¬ 
stituting  in  the  obvious  extension  of  (8.1l)  to  obtain  ^o^H* 

<\J0alCJ2^>H  »  and  <^aoala2a3a/As  we  finally  obtain 

3<5>>H  “  4<^0ala2^H  "  2<^a0ala3^i  +  3<^0aia2a3a^ 

(9.10) 

=  tanh  a [3  -  “  2<^1°3/H  +  ^  ^1  a2 a3 

Inserting  (9.10)  in  (9.7)  and  taking  limits,  we  get 


lira  limX-oKT  ,H  )  =  ftr  [3  -  2(2<V2>  + 

H  — >  0  uj->0 
z 


+  3<a1a203<:,v>]  ^ 


mH 

tanh  — -5- 
kT 

2mHr 


(9.H) 


3  H  ->  o 

z  *c 


i.e.  —  =  K[3  -  2(2<V2>  +  <W>  +  3<W3^.  (9.12) 

Nrn 


and  the  correlations  in  (9.H)  and  (9.12)  are  evaluated  in  zero  field. 
Then  substituting  (8.29)  in  (9.12) »  we  immediately  obtain. 
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Xg(O.T.O)  =%®(T),  (9.13) 

where^^T)  is  given  in  (9.1). 

This  result  confirms  our  conclusions  regarding  the  distinction  be¬ 
tween  the  adiabatic  susceptibility  calculated  by  the  Kubo  formalism  and 
the  isothermal  susceptibility  evaluated  by  Fisher  for  the  plane  square 
lattice.  It  is  important  to  note  that,  for  lattices  of  odd  coordination 
number,  one  does  not  obtain  terms  proportional  to  sin  60jt  alone  as  in 
(9.5)  and  hence  no  additional  contributions  will  be  found  by  reversing 
the  limits. 

In  Figure  9*1  we  plot  for  the  plane  square  lattice  the  isothermal 
susceptibility  given  by  Fisher's  (1963)  equation  (5.^2),  the  zero  fre¬ 
quency  adiabatic  susceptibility  in  our  equation  (8.25)  together  with 
their  differenceX^T) .  We  see  that  the  difference  is  appreciable  only 
in  the  critical  region.  Certainly  that  the  two  susceptibilities  are 
equal  at  zero  temperature  is  consistent  with  the  non-existence  in  the 
completely  magnetized  state  of  any  'free5  spins. 


The  temperature  dependence  of  the  zero  frequency  initial  per¬ 
pendicular  susceptibility  of  the  quantum  mechanical  Ising  model  on  the 
square  lattice 

(a)  the  isothermal  susceptibility fX^_(T) 

(b)  the  adiabatic  susceptibility%^(0,T)  and 

(c)  the  difference  between  the  two  susceptibilities. 
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CHAPTER  X 

INELASTIC  NEUTRON  SCATTERING  FROM  AN  ISING  MAGNET 


As  a  final  example  of  the  application  of  the  techniques  we  have 
just  described  we  shall  consider  the  problem  of  the  inelastic  scattering 
of  neutrons  from  an  Ising  system. 

Since  Van  Hove  (195*0  first  point  out  and  explicitly  stated  the 
connection  between  neutron  scattering  from  ferromagnetic  crystals  and 
time-dependent  spin  correlations,  there  has  been  considerable  theoreti¬ 
cal  interest  in  establishing  microscopic  theories  of  neutron  magnetic 
scattering.  One  of  the  primary  aims  has  been  to  describe  the  phenome¬ 
non  of  critical  scattering,  the  large  increase  in  the  scattering  ampli¬ 
tude  that  occurs  near  the  critical  temperature  of  a  ferromagnet  (and 
antiferromagnet ) .  Notable  contributions  have  been  made  by  Elliott  and 
Marshall  (1958)  and  more  recently  by  Marshall  (1965)  and  Fisher  and 
Burford  (1967)  all  of  which  contain  many  other  references  to  relevant 
work. 

The  problem  of  inelastic  scattering  has  been  specifically  con¬ 
sidered  for  the  Heisenberg  model  by  De  Gennes  (1958)  and  for  the  Ising 
model  by  Collins  and  Windsor  (1967)  and  Sears  (1967a).  Both  papers 
deal  only  with  the  limiting  case  of  infinite  temperature.  Very  recently, 
by  considering  the  effect  of  the  short  range  order  on  the  scattering 
from  a  Heisenberg  magnet,  Sears  (1967b)  has  extended  De  Gennes' s  original 


work.  However,  it  has  not  been  possible  thus  far  to  make  accurate  pre¬ 
dictions  regarding  the  critical  scattering  from  a  Heisenberg  system 
where  fluctuations  in  the  long  range  order  are  important.  In  contrast, 
we  shall  demonstrate  that  in  the  much  simpler  Ising  model,  it  is  possible 
to  determine  the  exact  energy  distribution  of  the  inelastic  scattering 
for  the  particular  case  of  two  dimensional  spin  one  half  systems  at  all 
temperatures. 

For  a  system  of  N  spins  rigidly  fixed  to  lattice  sites,  the  cross 
section  for  the  scattering  of  neutrons  per  unit  solid  angle  O. per  unit 
energy  E  is  given  by  (Van  Hove  195*0 » 


,2  v_2  2..  k. 
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where  e  and  m  are  the  electron  charge  and  mass,  y  =  1.91  is  the  neutron 
magnetic  moment  in  nuclear  magnetons,  c  is  the  velocity  of  light,  k^ 
and  k^  =  JkQ  -  K/  are  the  initial  and  final  wave  vectors  of  the  neutron, 
-ftoo  = -h(kQ  -  k^)/2iriQ,  mQ  being  the  mass  of  the  neutron,  is  the  energy 

A 

transferred  from  neutron  to  spin  system  and  K  is  the  octh  direction 
cosine  of  K' .  f(K' )  is  the  usual  atomic  form  factor,  the  fourier  trans- 
form  of  the  electron  spin  density,  and  S^K*  ,w)  is  given  by, 


+  00 

g°c3( ,^)  =  J=-  2  /exp[i(Kj.R-a)t)]<c^(0)^(t^>dt  .  (10.2) 
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(t),  the  |d  component  of  the  spin  on  lattice  site  R  at  time  t  is  of  the 


form, 


a^(t)  =  exp(i}/0t/h)a^exp(-i7^t/ft)  ,  (10.3) 


with^  given  by  (7.6). 

As  pointed  out  by  Collins  and  Windsor  (1967 )*  for  the  Ising  model 
gCXb^1  ^  for  a  ^  (3  and  so  the  "off  diagonal"  elements  make 

no  contribution  to  the  scattering  cross  section  (10.1).  Simple  inspec¬ 
tion  also  shows  that  SZZ(J£*  ,Uj)  is  responsible  for  the  critical  elastic 
scattering  and  has  been  very  completely  investigated  by  Fisher  and  Bur- 
ford  (1967).  In  contrast  to  the  Heisenberg  model,  the  inelastic  scatter¬ 
ing  is  solely  determined  by  SXX(K/I  ,uS)  =  S^(KJ  ,ib)  which  is  essentially 
what  Collins  and  Windsor  calculated  in  the  infinite  temperature  limit. 

We  shall  proceed  to  determine  S  (K3  ,uu)  exactly  for  all  temperatures  by 
techniques  very  similar  to  those  used  in  previous  chapters. 

Substituting  (7.34)  in  (10. 3)  we  obtain,  after  some  manipulation, 


<a0aXR(t)>=  cosV<(a0aR  TV1 


Z  Z 

+  ia  a  tan  tuAt 
k  R  0 


(10.4) 


where,  as  before,  U)q  =  2J/K.  (10.4)  and  all  future  remarks  unless 

otherwise  stated  refer  to  the  particular  case  of  the  plane  square  lat¬ 
tice  (or  any  lattice  for  which  q  =  4).  However,  an  exactly  similar 
analysis  could  be  undertaken  for  any  known  lattice.  It  is  clear  that 
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for  the  Ising  model  as  given  in  (7.6)  non  zero  correlations  can  exist 
only  between  the  z  components  of  spins  on  different  lattice  sites.  Hence 
(10.4)  is  non  zero  only  for  R,  =  0  and  on  expanding  the  right  hand  side, 
we  obtain 

<aoao(t)>  =  cosV^t  +  4i  cos-\oQt  sinouQt<aQCJ^> 

-  2  cos2a^t  sin2(j^t [Z (0^2^  +  (10.5) 

-  4i  cosCu0t  sin^t<&0a1o2a^)>+  s in^  t ^ a2 a ^ 

where  as  before  we  have  dropped  the  superscript  z  from  the  spins  on  the 
right  hand  side  and  1,  2,  3.  4  label  the  nearest  neighbour  sites  to  the 
zeroth.  Substitution  of  (10.5)  in  (10.2)  then  gives 

SXX(JL'  .u).T)  =  6(co)[|-i  (2^2)  +  <^a^>)  +|  <^0,030,53 

+  \  a(co-2<^)[i+  <oQo^-  i  <0^2°?^)- <^aoal°2a3/^ 

+  \  ^ -  V4  +  Ko°3>  +  l  (2<?la2>  +  <ala3> 

+  |  <a0ala2a3>  +  I^al02a3‘74^  +  2  6  (t0+2^)t2  -  <O0aV> 


-  |<^la2°3a4>+  <a0al°2a3^  +  l6  <iUJ+4lii)t8_Ka0al> 
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+  l  (2<?la2>  +  <ala3>  +l<$la20304>  -  2  <a0ala2a3>^  •  (10<6) 

The  temperature  dependence  of  all  the  spin  correlations  in  (10.6)  for  the 
plane  square  lattice  have  been  discussed  in  chapter  eight  and  (10.6)  then 
constitutes  an  exact  solution  of  both  the  temperature  dependence  and 
energy  distribution  of  the  inelastic  scattering  of  neutrons  from  the 
plane  square  spin  one  half  Ising  lattice. 

In  contrast  to  the  Heisenberg  model,  the  energy  distribution  con¬ 
sists  of  a  number  of  discrete  levels,  each  level  corresponding  to  a  par¬ 
ticular  orientation  of  the  four  nearest  neighbours  to  any  given  spin. 

For  instance,  as  before,  the  6(Cu-46i^)  contribution  is  associated  with 
the  configuration  in  which  all  four  spins  are  parallel  with  the  central 
one  only  now  the  four  antiparallel  neighbouring  spin  environment  makes  a 
physically  meaningful  contribution  to  the  6(60+40^)  peak.  The  former 
represents  absorption  of  energy  by  the  spins,  the  latter  by  the  neutrons. 

The  temperature  dependent  amplitudes  of  these  levels  given  in 
(10.6)  in  terms  of  spin  correlation  functions  are  proportional  to  the 
probability  that  the  associated  spin  orientation  may  occur.  For  instance, 
as  we  might  expect  by  now,  at  T  =  0  only  the  Qj  =  +  4(J^  level  has  non 
zero  amplitude  since  in  the  completely  magnetized  state  only  the  all 
parallel  configuration  will  arise.  Also  the  fact  that  the  energy  distri¬ 
bution  is  symmetric  in  only  the  infinite  temperature  limit  is  consistent 
with  this  physical  interpretation.  Indeed  since  all  spin  correlations 
tend  to  zero  in  this  limit,  we  see  by  inspection  of  (10.6)  that  theUj=0, 
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60=  +  2U^  and  00=  +  4lJ^  levels  are  in  the  ratio  of  6  to  4  to  1.  This 
is  precisely  the  ratio  of  the  probabilities  of  occurrence  of  the  corre¬ 
sponding  spin  configurations  in  a  random  array  of  spins.  From  equations 
(8.35)  to  (8.38),  we  recall  that  all  the  spin  correlations  have  infinite 
slope  at  the  critical  temperature  and  vary  very  little  outside  a  rela¬ 
tively  narrow  region  around  T  .  This  means  that  the  energy  distribution 
varies  appreciably  only  in  the  critical  region.  Many  of  the  above  fea¬ 
tures  are  illustrated  in  Figure  10.1  where  Sxx'( K*  ,UJ,T)  is  plotted  versus 

CO  for  T  =  0,  T  =  0.8  T  ,  T  =  T  ,  T  =  2  T  and  T4». 

c  c  c 

Much  that  we  have  said  about  the  scattering  from  two  dimensional 
Ising  magnets  and  the  plane  square  lattice  in  particular  will  also  be 
characteristic  of  scattering  by  three  dimensional  Ising  magnets.  Firstly 
we  should  emphasize  once  again  that  an  expression  of  the  form  (10.6)  for 
SXX(KI  ,uj,T)  can  be  found,  ■without  difficulty  for  a  regular  lattice  of 
any  coordination  number  and  irrespective  of  dimension.  For  three  dimen¬ 
sional  lattices,  however,  the  correlations  cannot  be  found  exactly  but 
the  high  and  how  temperature  series  expansions  for  them  could  be  calcu*. 
lated.  Then  the  entire  temperature  dependence  of  each  would  be  well 
represented  by  Pade  approximant  interpolation  and  extrapolation.  We 
know  further  from  the  work  of  Sykes,  Martin  and  Hunter  (1967)  that  close 

to  T  for  T  >  T  , 
c  c 

CH  a(T  -  T0)-l/8  .  (10.7) 

H40 

This  means  that  for  the  nearest  neighbour  correlation which  is 


Fig.  10.1 


Inelastic  neutron  scattering  amplitudes  from  a  plane 
square  Ising  lattice  for  T  =  0,  T  =  0.8T  ,  T  =  T  ,  and 
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proportional  to  the  energy, 


(10.8) 


for  the  same  temperature  range  (C-^  and  C ^  are  constants).  The  same  type 
of  rapid  temperature  variation  is  expected  for  all  nearest  neighbour  cor¬ 


relations  in  three  dimensions  and  hence  for  S^CK*  ,uj,T)  . 


As  in  two  dimensions,  energy  can  be  lost  or  gained  by  the  scat¬ 
tered  neutrons  only  in  amounts 


(10.9) 


E  =  +diUJ=  2qJ,(2qJ  -  4)  J, 


0  9  0  9 


We  therefore  expect  that  qualitatively  the  temperature  and  energy  depend¬ 
ence  of  the  scattered  neutrons  will  be  the  same  in  three  dimensions  as 
in  two. 

Finally  the  effect  of  interactions  with  neighbours  farther  than 

the  first  is  simply  to  introduce  further  sharp  resonances.  On  the  other 

hand  the  effect  of  small  crfa^  and  type  interactions,  to  be  expected 

1  J  i  J 

in  any  real  crystal,  would  be  principally  to  round  out  the  sharp  reso¬ 
nances  in  energy  distribution  and  also  introduce  slight  KJ  dependence  in 
SXX(K’ ,oj,T)  because  in  that  case  correlations  between  the  x  or  y  com¬ 
ponents  of  spins  on  different  lattice  sites  would  no  longer  be  zero. 
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CHAPTER  XI 

FINAL  REVIEW  AND  DISCUSSION 


We  have  considered  in  this  thesis  two  problems  whose  many  dif¬ 
ferences  illustrate  both  the  variety  of  lattice  models  of  critical 
phenomena  and  the  many  distinct  techniques  that  have  been  developed  to 
investigate  them. 

A.  Spin  One  Exchange  Model 

In  the  spin  one  exchange  model,  we  have  examined  the  equilibrium 
critical  behavior  of  a  model  with  a  completely  isotropic  Hamiltonian. 

We  have  already  noted  in  chapter  one  that  short  range  interaction  mod¬ 
els  (usually  just  nearest  neighbours)  give  results  which  approximate 
very  closely  to  experiment.  Also  these  results  are  quite  different  from 
those  obtained  from  the  mean  field  calculations  which  are  equivalent  to 
systems  with  infinitely  long  range  interactions.  From  the  theoretical 
point  of  view,  it  seems  clear  that  the  range  of  the  interaction  is  an 
important  parameter  controlling  the  macroscopic  critical  behavior  of  the 
system.  The  question  still  remains,  however,  of  how  relevant  are  the 
details  of  the  interaction,  i.e.  whether  it  is  pure  Heisenberg  or  spin 
one  exchange  for  instance,  when  the  range  is  kept  constant.  This  is 
largely  the  interest  in  our  calculation  on  this  model  by  means  of  which 
we  hope  to  shed  some  light  on  this  question. 


. 
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Because  of  the  inherent  complexity  of  the  calculation,  it  has  not 
been  possible  to  obtain  an  exact  solution.  The  technique  of  high  tem¬ 
perature  series  expansions  is,  however,  a  very  precise  means  of  obtain¬ 
ing  successive  approximations.  The  series  in  powers  of  inverse  tempera¬ 
ture  represent  an  expansion  about  the  infinite  temperature  limit  where 
succeeding  terms  take  into  account  the  increasing  order  in  the  system 
as  the  temperature  is  lowered.  The  aim  is  to  determine  exactly  as  many 
of  the  coefficients  of  the  various  series  as  is  reasonably  possible. 
These  truncated  series  in  themselves  fail  to  represent  precisely  the 
corresponding  thermodynamic  variable  for  temperatures  very  close  to 
critical  but  extrapolation  techniques  such  as  the  Domb  ratio  method  and 
the  Pade  approximant  method  have  succeeded  in  estimating  the  critical 
behavior  with  considerable  precision. 

There  exist  two  methods  of  obtaining  the  series  coefficients, 
the  direct  method  and.  the  linked  cluster  method,  and.  we  have  shown  how 
the  latter  is  essentially  a  rearrangement  of  the  former.  Nevertheless, 
as  we  have  emphasized,  it  is  a  much  less  complicated  and  more  system¬ 
atic  technique  to  apply  particularly  with  the  further  simplifications 
introduced  by  using  the  representation  theory  of  the  symmetric  group. 

By  this  method  we  have  obtained  longer  series  than  have  previously  been 
calculated  for  a  quantum  mechanical  spin  one  system. 

The  numerical  values  we  obtained  for  the  critical  temperature, 
the  critical  entropy  and.  y,  the  high  temperature  critical  index  of  the 
magnetic  susceptibility  clearly  indicate  that  there  are  small  but  sig- 


. 
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nificant  differences  in  critical  behavior  which  can  be  directly  attributed 
to  the  details  of  the  nearest  neighbour  interaction.  It  would  seem  then 
that  any  theory  based  on  the  assumption  that  such  details  are  not  im¬ 
portant  (e.g.  Kadanoff  1966)  can  only  achieve  qualitative  and  not  quanti¬ 
tative  consistency  with  the  model  calculations. 

Using  Domb  and  Miedema's  (1964)  and  Heller's  (196?)  experimental 
reviews  as  our  main  references,  we  know  of  no  spin  one  insulating  ferro- 
magnet  whose  experimental  behavior  might  be  compared  to  our  results. 

This  is  not  to  say  that  such  substances  do  not  exist  or  may  not  be  dis¬ 
covered.  If  such  experiments  were  performed,  we  certainly  would  expect 
that  some  of  the  trends  our  results  show  in  comparison  with  the  spin  one 
Heisenberg  model  would  be  manifest  in  some  of  these  substances. 

A  recent  controversy  has  surrounded  the  possibility  of  a  non-zero 
critical  temperature  in  the  two-dimensional  Heisenberg  model.  On  the 
one  hand,  Mermin  and  Wagner  (1966)  have  proved  rigorously  that  a  two- 
dimensional  Heisenberg  model  cannot  have  a  spontaneous  magnetization  or 
sublattice  magnetization.  On  the  other,  Stanley  and  Kaplan  (1966)  have 
pointed  out  that  the  evidence  from  high- temperature  series  expansions 
for^Q  for  the  Heisenberg  model  in  two  dimensions  favours  a  non-zero 
critical  temperature,  especially  for  S  >  ■§•,  We  have  examined  the  high 
temperature  series  of  the  spin  one  exchange  model  for  the  triangular 
lattice,  the  two  dimensional  lattice  with  the  highest  coordination  num¬ 
ber  and,  therefore,  presumably  the  most  regular  behavior.  The  results 
of  the  analysis,  based  on  longer  series  than  Stanley  and  Kaplan  have 


. 
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available  for  the  Heisenberg  model,  are  presented  in  appendix  F.  Our 
interpretation  is  that  with  the  series  available,  there  is  no  evidence 
for  the  existence  of  a  critical  temperature  in  two  dimensional  spin  one 
exchange  models. 

There  remains  at  least  one  interesting  theoretical  question  about 

spin  one  systems.  Just  as  there  are  two  linearly  independent  isotropic 

2 

interactions  for  spin  one  particles  e.g.  S..S.  and  (S..S.)  so  there  are 

**  i  ^ j  **  1 

two  independent  long  range  order  parameters  suggested  by  these  inter¬ 
actions.  The  Heisenberg  interaction  S . .S.  between  nearest  neighbour 

l  j 

pairs  suggests  the  usual  long-range  order  parameter. 


h  =  <§0-^>=  3<sos»> 


(11.1) 


=  3 


-2 

M_ 

m2 


where  the  subscripts  indicate  the  limit  of  infinite  separation  of  the 
two  spins  and  M  is  the  spontaneous  magnetization  per  spin.  The  second 
long  range  order  parameter  may  be  taken  to  be  <(^Sq.  but  the  prefera¬ 

ble  definition  for  both  computational  and  theoretical  reasons  is 

R2=<\P0°°/>  (11.2) 

where  P0oo is  the  exchange  operator  for  two  particles  in  the  limit  of 
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infinite  separation.  For  spin  one  half  systems  the  two  definitions  are 
of  course  linearly  related. 

In  the  high  temperature  region  (T  >  T  )  all  long  range  order  is 

c 

absent  but  the  fluctuation  in  the  long  range  order  parameters  is  non¬ 
zero  and  becomes  very  large  just  above  Tc*  In  fact  the  fluctuation  Ar ^ 
is  directly  proportional  to  the  initial  susceptibility  which  we  have 
just  determined.  It  would  be  interesting  then  to  investigate  AR^,.  Pre- 
sumablyARp  also  has  a  singularity  at  T  ^  with  some  critical  index  y  but 
it  is  not  clear  that  it  would  be  equal  to  y  for  the  spin  one  exchange 
model . 

B,  The  Perpendicular  Susceptibility  of  the  Ising  Model 

In  complete  contrast,  in  this  problem  we  have  calculated  a  time- 
dependent  macroscopic  variable  for  a  system  whose  nearest  neighbour 
interaction  is  highly  anisotropic.  Indeed  it  is  this  anisotropy  that 
is  largely  responsible  for  the  time  dependent  behavior.  Again,  in  con¬ 
trast  to  the  spin  one  exchange  model,  the  Ising  interaction  is  suffi¬ 
ciently  elementary  that  we  have  been  able  to  calculate  the  frequency 
dependent  initial  perpendicular  susceptibility  exactly  for  two  dimen¬ 
sional  systems. 

The  formalism  we  have  used,  the  Kubo  linear  response  theory,  is 
particularly  useful  for  studying  the  response  of  a  system  to  a  weak 
external,  mechanical  disturbance.  The  linear  approximations  embodied 
in  the  theory  become  exact  in  our  case  of  the  external  field  tending  to 
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zero. 

In  comparing  our  results  with  Fisher's  static  calculations,  we 
have  emphasized  the  important  distinction  between  the  susceptibility  of 
an  isolated  system  and  of  one  in  contact  with  a  heat  bath  or  some  agency 
that  allows  the  system  to  reach  its  equilibrium  state  in  the  presence  of 
the  external  magnetic  field.  In  particular  we  found  that  for  the  simple 
quadratic  lattice,  and  most  likely  all  even  coordination  number  lattices, 
there  exists  a  quantitative  difference  which  can  be  directly  attributed 
to  whether  the  system  is  isolated  or  not. 

The  point  can  also  be  made  that  in  comparison  with  Glauber' s  (1963) 
master  equation  approach,  our  results  contain  no  additional  unknown  pa¬ 
rameters.  Any  time  dependent  variable  of  a  system  in  contact  with  a 
heat  bath  will  contain  properties  that  depend  not  only  on  the  basic 
interaction  between  the  constituents  of  the  system  (e.g.  the  Ising  inter¬ 
action)  but  also  on  the  nature  of  the  external  contact.  This  is  reflected 
in  the  presence  of  unknown  parameters  in  Glauber's  results.  However,  our 
results  both  for  the  perpendicular  susceptibility  and  the  neutron  scat¬ 
tering  cross-section  are  wholly  dependent  on  the  nature  of  the  Ising 
interaction.  As  a  consequence,  similar  results  say  for  the  Heisenberg 
model  will  be  completely  different  in  nature  and  will  provide  a  very 
good  basis  on  which  to  decide  whether  a  given  system  has  or  has  not  any 
anisotropy  in  its  basic  interaction.  Indeed,  the  frequency  dependent 
perpendicular  susceptibility  differentiates  between  such  models  to  a  much 
greater  extent  than  do  the  static  parallel  susceptibilities. 
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Again  the  experimental  data  on  the  perpendicular  susceptibility 
of  Ising  systems  is  very  limited.,  Fisher  (1965)  points  out  some  work 
performed  by  Lasheen  et  al.  (1958)  on  the  isothermal  perpendicular 
susceptibility  of  the  antiferromagnet  MnCl^H^O.  At  zero  frequency  there 
is  qualitative  agreement  only  between  his  results  (and  therefore  ours) 
and.  experiment.  However,  much  of  the  discrepancy  may  be  due  to  the  fact 
that  MnClgH^O  inevitably  only  approximates  to  an  Ising  system. 

It  seems  unlikely  that  exact  results  will  ever  be  obtained  both 
for  the  frequency  dependent  susceptibility  and.  the  neutron  scattering 
cross  section  for  three  dimensional  Ising  or  Heisenberg  systems.  How¬ 
ever,  it  would,  seem  a  useful  approach  to  apply  the  high  temperature 
series  expansion  techniques  we  have  described  in  Part  I  together  with 
the  Kubo  formalism  we  introduced  in  Part  II  to  make  much  more  precise 
approximate  calculations  than  have  thus  far  been  attempted. 

The  work  we  have  just  described  in  this  thesis  is  also  reported 
in  a  series  of  four  papers,  one  of  which  deals  with  the  spin  one  exchange 
model  (Allan  and  Betts  1967),  one  with  the  inelastic  neutron  scattering 
(Allan  and  Betts  1968c)  and  two  with  the  frequency  dependent  initial 
perpendicular  susceptibility  (Allan  and  Betts  1968a),  b)). 

We  should  also  like  to  mention  four  recent  papers  by  other 
authors  which  either  make  reference  to  the  above  publications  or  have 
arisen  partly  as  a  result  of  private  communication  between  the  authors 
and  ourselves.  Seagraves  (1966)  gives  an  alternative  derivation  of  the 
spin  representation  of  the  permutation  operator  for  spin  l/2,  1,  3/2 
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and  2.  Joseph  (1967)  has  considered  the  problem  of  the  exchange  inter¬ 
action  for  general  spin  in  comparison  with  Heisenberg  interaction.  Falk 
(1968)  demonstrates  that  the  adiabatic  susceptibility  is  a  lower  bound 
to  the  isothermal  susceptibility  and  uses  our  calculation  in  part  II  of 
this  thesis  as  an  example.  Finally  Kim  (1968)  constructs  a  frequency 
dependent  susceptibility  which  equals  the  static  isothermal  suscepti¬ 
bility  in  the  sero  frequency  limit  and  agrees  with  our  adiabatic  sus¬ 
ceptibility  for  non-zero  frequencies. 
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The  curnulant-moment  identities  for  p.  ~  0, 
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^  -  36)i?p2  -  84^3 

-  7560>i3>i| 


2  3 

126^14  +  756)1^2  +  2520^)1^2  +  5  6  Op., 


A.2 


*5-)* 2 


,  II  II 

^3  ~  ^3 


^4  =  H 


^  ^  -  10(^2  +  Vz?p 

A6  =  ^6  "  15  +  ^4^  “  20>13)13  +  90>^)i2 


= 


= 


8 


A9  ~ 


y'y  -  21(^2  +  y-2^)  -  35(^3  +  313)14)  +  210(2)13)12)12  +  YZ^ 

>4  “  2 8 ^6^2  +  >12>16)  “  36(313} 13  ■!*  >13)13)  -  70)iZj}i4  -i-  420(^1 

11  2  11  11  o  11 

+  2>i2p/l>l2 )  +  560  +  2^1^2)13)  -  252031^2 

F9  -  36(}1^2  +  ~  +  ^3)  -  126 (>15)14  + 

+  756(^30^  +  2313)12)12)  +  252  +  v-ipjYz') 

+  1680)13)13  -  7560(31^)12  +  3313)12 P2) 


B.l 


where 


APPENDIX  B 


Decomposition  of  the  Permutation  Operator 


n„n 
a.  a . 

i  J 


4  =  (^)2  - 1 


a2  =  (Sy)2  -  2 
i  v  i'  3 


4 


z\2 


(sp 


2 

3 


a.  = 


sK 


a,  = 


SySX 

ii 


8 

a.  = 


z  y 
s.s. 

i  l 


a.  = 


Qx  z 
SiSi 


a10 

ai 


= 


sysz 

1  1 


6 

ai 


Z  X 

sisi 


S 


z 

i 


SI 

C.l 


APPENDIX  C 


Character  Tables  of  the  Irreducible  Representations  of  TT  ^ — *  TT  ^ 


Notation:  j~ ^  represents  the  3^  dimensional  reducible  direct  product  rep¬ 
resentation  of  TF^.  Only  non  zero  n(y)  are  listed. 


ti2 


CLASS 

- 

2 

ORDER 

1 

1 

(2) 

1 

1 

d2) 

1 

-1 

r 

1  2 

9 

3 

n(2)  =  6 

n(l2)  =  3 


CLASS 

i3 

- —  —  . — ' 

1,2 

3 

ORDER 

i 

3 

2 

(3) 

i 

1 

1 

(2,1) 

2 

0 

-1 

(l3) 

1 

-1 

1 

r3 

27 

9 

3 

n(3)  =  10 

n(l3)  =  1 
n(2,l)  =  8 


CLASS 

l4 

12»2 

1,3 

4 

22 

ORDER 

1 

6 

8 

6 

3 

(4) 

1 

1 

1 

1 

1 

(3,1) 

3 

1 

0 

-1 

-1 

(22) 

2 

0 

-1 

0 

2 

(2,12) 

3 

-1 

0 

1 

-1 

(l4) 

1 

-1 

1 

-1 

1 

[7 

81 

27 

9 

3 

9 

nW  =  15 
n(22)  =  6 


n(3,l)  =  15 

n(2,l2)  =  3 


CLASS 

l5 

13,2 

l2  *3 

1,4 

1.22 

2,3 

5 

ORDER 

1 

10 

20 

30 

15 

20 

24 

(5) 

1 

1 

1 

1 

1 

1 

1 

(4,1) 

4 

2 

1 

0 

0 

-1 

-1 

(3,2) 

5 

1 

-1 

-1 

1 

1 

0 

(3.12) 

6 

0 

0 

0 

-2 

0 

1 

(22,1) 

5 

-1 

-1 

1 

1 

-1 

0 

(2.13) 

4 

-2 

1 

0 

0 

1 

-1 

(l5) 

1 

-1 

1 

-1 

1 

-1 

1 

r, 

243 

81 

27 

9 

27 

9 

3 

n(5)  =  21 

n(4,l)  =  24 
n(22,l)  =  3 
n(3,2)  =  15 
n(3.12)  =  6 


C.4 


Tf6 


CLASS 

I6 

14,2 

13,3 

12,4 

1^22 

1.2,3 

1,5 

6 

2,4 

23 

32 

ORDER 

1 

15 

40 

90 

45 

120 

144 

120 

90 

15 

4 

(6) 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

(5.D 

5 

3 

2 

1 

1 

0 

0 

-1 

-1 

-1 

-1 

(4.2) 

9 

3 

0 

-1 

1 

0 

-1 

0 

1 

3 

0 

(4.12) 

10 

2 

1 

0 

-2 

-1 

0 

1 

0 

-2 

1 

(32) 

5 

1 

-1 

-1 

1 

1 

0 

0 

-1 

-3 

2 

(3.2.1) 

16 

0 

-2 

0 

0 

0 

1 

0 

0 

0 

-2 

(23) 

5 

-1 

-1 

1 

1 

-1 

0 

0 

-1 

3 

2 

(3.13) 

10 

-2 

1 

0 

-2 

1 

0 

-1 

0 

2 

1 

(22,12) 

9 

-3 

0 

1 

1 

0 

-1 

0 

1 

-3 

0 

(2.14) 

5 

-3 

2 

-1 

1 

0 

0 

1 

-1 

1 

-1 

(l6) 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

"*  1 

-1 

1 

f6 

729 

243 

81 

27 

81 

27 

9 

3 

9 

27 

9 

n(6)  =  28 

n(5.1)  =  35 

n(4,2)  =  27 

n(4,l2)  =  10 

n(32)  =  10 

n(23)  =  1 


n(3,2.l)  =  8 


.  o 

• 

C.5 


CLASS 

17 

1?2 

4 

173 

1?4 

1?22 

1?5 

1,6 

124 

123 

132 

2,5 

2 

2,3 

3,4 

7 

ORDER 

1 

21 

70 

210 

105 

420 

504 

840 

530 

105 

280 

504 

210 

£0 

720 

(7) 

l 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

(6,1) 

6 

4 

3 

2 

2 

1 

1 

0 

0 

0 

0 

-1 

-1 

-1 

-1 

(5,2) 

14 

6 

2 

0 

2 

0 

-1 

-1 

0 

2 

-1 

1 

2 

0 

0 

(5.12) 

15 

5 

3 

1 

-1 

-1 

0 

0 

-1 

-3 

0 

0 

-1 

1 

1 

(4,3) 

14 

4 

-1 

-2 

2 

1 

-1 

0 

0 

0 

2 

-1 

-1 

1 

0 

(4,2.1) 

35 

5 

-1 

-1 

-1 

-1 

0 

1 

1 

1 

-1 

0 

-1 

-1 

0 

(32.1) 

21 

1 

-3 

-1 

1 

1 

1 

0 

-1 

-3 

0 

1 

1 

-1 

0 

(4,l3) 

20 

0 

2 

0 

-4 

0 

0 

0 

0 

0 

2 

0 

2 

0 

-1 

(3,22) 

21 

-1 

3 

1 

1 

-1 

1 

0 

-1 

3 

0 

-1 

1 

1 

0 

(3,2, l2) 

35 

-5 

1 

1 

-1 

1 

0 

-1 

1 

-1 

1 

0 

-1 

1 

0 

(23,1) 

14 

-4 

-1 

2 

2 

-1 

-1 

0 

0 

0 

2 

1 

-1 

-1 

0 

(3. I4) 

15 

-5 

3 

-1 

-1 

1 

0 

0 

-1 

3 

0 

0 

-1 

-1 

1 

(22.13) 

14 

-6 

2 

0 

2 

0 

-1 

1 

0 

-2 

-1 

-1 

2 

0 

0 

(2,15) 

6 

-4 

3 

-2 

2 

-1 

1 

0 

0 

0 

0 

1 

-1 

1 

-1 

(l7) 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

-1 

1 

m? 

2187 

729 

243 

81 

243 

81 

27 

9 

27 

81 

27 

9 

27 

9 

3 

n(7)  =  36 

n(6,l)  =  48 

n(5,2)  =  42 
n(5.12)  =  15 


n(4,3)  =  24 

n(4,2,l)  =  15 
n(32,l)  =  6 
n(3,22)  =  3 


C.6 


Distribution  of  Irreducible  Representations 

2 

Among  the  Eigenstates  of  S 


N  a 

2 

For 

S  =  2, 

the  irreducible  representations 

are 

(2). 

ii 

S  =  1, 

it 

it 

ti 

ii 

(I2). 

n 

S  =  0, 

it 

it 

ti 

n 

(2). 

N  = 

2 

For 

s  =  3, 

the  irreducible  representations 

are 

(3), 

ti 

S  =  2, 

ii 

ti 

n 

ti 

(2.1), 

it 

S  =  1, 

n 

ii 

ti 

ii 

(2.1)  + (3), 

ii 

S  =  0, 

ti 

it 

ii 

ti 

(l3). 

N  a 

4 

For 

S  =  4, 

the  irreducible  representations 

are 

(4), 

ii 

S  =  3, 

n 

it 

ii 

ti 

(3.D. 

it 

S  a  2, 

ii 

n 

ii 

ti 

(3,l)  +  (22)  +  (4), 

ti 

s  =  1, 

it 

it 

it 

ti 

(3,1) + (2. I2). 

n 

S  =  0, 

ii 

it 

ii 

it 

(2,l2)  +  (4). 

N  = 

A 

For 

s  =  5, 

the  irreducible  representations 

are 

(5), 

ii 

S  =  4, 

ii 

it 

ii 

ti 

(4,1). 

it 

s  *  3. 

ii 

ti 

ti 

ii 

(3,2)  + (4,1)  + (5) 

n 

S  =  2, 

n 

it 

ii 

ti 

(3.12) + (3.2) + (4 

ii 

S  a  1, 

ti 

ii 

n 

it 

(3,2)  +  (22,1)  +  (4 

n 

S  a  0, 

ii 

n 

ti 

ii 

(3.12). 

C.7 


N  =  6 

For  S  =  6,  the  irreducible  representations  are  (6) , 


It 

s  =  5, 

tl 

it 

ii 

it 

(5.1). 

It 

S  =  4, 

ll 

ti 

it 

it 

(4,2)  +  (5.1)  +  (6), 

II 

C3 

II 

C/3 

li 

ii 

n 

ii 

(4,l2)  +  (4,2)  +  (5.1)  +  (32). 

II 

S  =  2, 

ll 

ii 

tt 

ii 

(3,2,l)+2(4,2)  +  (5*l)  +  (6), 

II 

Sal, 

ll 

ii 

ii 

it 

(3.2,1) +  (4,12)  + (5,1) +  (32) 

It 

S  =  0, 

ll 

ti 

ii 

ii 

(4,2)  +  (2-))  +  (6) , 

N  = 

2. 

For 

s  =  7, 

the  irreducible  representations 

are 

(7), 

it 

s  =  6, 

it 

ti 

ti 

it 

(6,1), 

tt 

S  a  5. 

it 

it 

n 

it 

(5,2) + (6,1) + (7) , 

ti 

S  =  4, 

ti 

ii 

it 

it 

(5,12) + (5.2) + (4,3) + (6,1), 

it 

S  a  3, 

n 

it 

ii 

it 

(4,2,1)  +2(5,2)  +  (4,3) 

+  (6,1) + (7) , 

ti 

S  =  2, 

ti 

n 

ii 

it 

(4,2,l)  +  (32,l)  +  (5,l2) 

+  (5,2) + (4,3) + (6,1), 

ti 

Sal, 

it 

it 

it 

ii 

(4,2,1) +  (3,22)  + (5,2) 

+  (4,3) + (6,1) + (7) » 
"  (32,D  +  (5.l2). 


"  S  =  0, 


(£,!')+  (S,e>s+  (.: .  5,  V) 


It 


D.l 


APPENDIX  D 


Cluster  Expansion  Data 


7 


K' 


n 


_  n  n 

n  =  0  3  n! 


Note 


1)  Where  two  series  are  given  for  any  cluster  C^,  the  upper 

series  is  and  the  lower  series  is  If  only  one  ser¬ 

ies  is  given,  then  it  isl^. 

2)  Only  non-zero  a^  are  given  in  the  table. 

3)  For  the  single  point  cluster  Cq,  we  have  IjJq  =  2/3  and 

C =  log  3»  For  the  single  bond  cluster  C^,  a^  =  4/ 3  in^. 


LABEL 

GRAPH 

a2 

3 

a4 

a5 

a6 

3.y 

ci 

»  - • 

-8/3 

+8 

-16 

-16 

+160 

-96 

+448 

+960 

-18,816 

+2,688 

+40,704 

-112,896 

C2 

A 

+8/3 

-  16 

-24 

-144 

+3,440 

+1,440 

-38,904 

+19,440 

-702,576 

-514,080 

C3 

A 

-  -  - 

+24 

+48 

-1488 

+288 

+1,080 

-13,680 

+513.360 

+10,080 

-2,005,416 

+3,922,128 

c4 

A/ 

+  8 

-48 

-720 

+17,040 

+7,776 

-266,112 

-181,440 

D.2 


D.3 


LABEL 

GRAPH 

a6  !  a7 

C2  7 

j 

N 

— — • 

A 

-6,144 

-234,528 

-169,344 

C28 

> 

-1-6,720 

-397,824 

-60,480 

C29 

<> 

+27,360 

+5.760 

-491,904 

+120,960 

1 

C30 

- 

\ - - 

► 

42,400 

-141,120 

C31 

f 

i _ , 

>  ■■ * 

I 

+2,400 

-184,800 

1 

C32 

► - 

i - 

t 

-2 , 880 

1 

-126,000 

j 

_ 

c33 

i 

—  i 

- 1 

A 

+4,320 

-209,832 

I 

- .  ,  .  .  .  .  t 

c34 

M 

+2,400 

-141,960  ; 

( 

1 

- [ 

°35 

-4,032 

C36 

N-J 

-2,880 

-41,328  | 

C37 

ZvV 

+1,920 

. 

-202,272  j 

D.5 


LABEL 

GRAPH 

a6 

a7 

oo 

r'v 

o 

X 

I 

1 

......  -  ... 

+272,160 

C39 

A> 

-45,360 

cw 

-84,672 

C4l 

AV\ 

+960 

0 

C42 

0 

C43 

AA/1 

-20,160 

C44 

/yv 

-10,080 

c45 

/ 

/ 

-  —  -  —  — 

-357,504 

-  66,528 

c46 

\ 

X 

—  —  -  —  - 

-2,135,952 

-665,280 

C47 

X 

> 

—  —  —  -  - 

-955,416 

+73.584 

c48 

X 

-272,160 

• 

•  ' 

D.6 


-.0'  t+  ■ 

.  v. 

- 

- 

.. 

D.7 


LABEL 

GRAPH 

a7 

LABEL 

GRAPH 

a7 

C71 

> - 

'\A 

+30,240 

C79 

A 

-20,160 

c72 

i - 

\ 

0 

C80 

X 

-20,160 

C73 

A 

V\ 

+13,' 440 

C8l 

+16,800 

C74 

> — 

( > - 

7 

-20,160 

C82 

:  x 

+16,800  : 

j 

—  1 

C75 

i — 

ii — 

X 

— 

1 

-• 

+16,800 

C83 

L 

>: 

1 

I 

+16,800  | 

C76 

0 

c84 

A 

-20,160 

C77 

■t- 

•  — - o- 

"l 

• 

+16,800 

C85 

/A 

+16,800 

C78 

A " 

- » 

0 

C86 

/WV 

+6,720 

Linear  Combinations  of  ^  » s 

Xq(C1(K)  =  2V0  +Vl 

X0(c2>k)  =  3Vo  +  2Vi  +v2 

OC/qCCj.k)  =  3Vo  +  3Vi  +  %  +  V3 

Xq(c4>k)  =  Wo  +  %  +  %  +(pk 


\i.  - 

-  -  ■  •• 

■  • 

-  ‘ 


D.  8 


'X^Cj.K)  : 

■  % 

+  3^ 

+  3fe 

+  V5 

X0(c6.iO  ; 

-  % 

•1-  Wj. 

4  4^2 

+  H 

+  V6 

X0(c7-k)  : 

‘  % 

4  4^ 

+  3V2 

+  V3  + 

2V4  "h 

V5  +y7 

'X0(Cg,K)  ! 

=  % 

+  % 

+  3^2 

+  *V4 

4  y8 

TCqCCj.k) 

=  5Vo 

+  ,% 

+  4(^2 

+  2^4 

4  V5  + 

V9 

X0(C10'K) 

=  5(fb 

4  4^ 

-1-  61^2 

+  % 

+  Yio 

k) 

=  l% 

+  % 

+ 

+  2yj 

+  6V4  ' 

f  21^  +  <V6 

4 

+  vn 

'Xq(c12,k) 

~  5Vo 

+  % 

+  6'V2 

+ v3  + 

■  + 

V5  +V7  + 

2^8  + 

V9  -^12 

=  % 

+  % 

+  7Vz 

4-  4 

•  5(^4  + 

2l|4j  4  2^r; 

4^8  ‘ 

f  2^9  4  V13 

=  5^0 

+  % 

+  a'h 

+  V3  -* 

■  /jV/j,  + 

4V5  4  2^ 

4  2^9 

+  Y10  +  Y14 

9(0  (ci5  *  k) 

=  5Vo 

-  % 

+  5V2 

+  5V4 

4  5Vg  ■ 

+  V15 

~  5Vo 

+  % 

+  % 

+  61^ 

4  V3  + 

^6  +  2^8  ' 

4  2^  • 

+  yi6 

^0(C17’K) 

=  6V0 

+  5V! 

4  41^ 

+  3V4 

4  2ljJg 

4  ljJ1? 

^d(ci8»k) 

=  6Yo 

+  514 

+  5y2 

+  3V4 

4  4 

2^8  4  ' 

+ Vis 

=  6V0 

+  5V! 

+  5V2 

+ 

+  ^5  + 

■  y8  4  2(jJ9 

+  ’+'19 

^C20’K^ 

=  % 

+  % 

+  7V2 

+  3iy4 

4  4l|}^ 

+  3^9  +  V10  +^20 

■ 

. 

■I  + 

J2 

D.9 


^(c21,k) 

=  6V0  +  %  +  6(f'2  +  L%  +  2V5  +  %  +  V21 

^0^22’^ 

=  4V0  +  6Vi  +  12  V 2  +  %  +  ^4  +  %  +  3%  +  12  V7  +  V8 

+  6l|J11  -^22 

X0(C23.k) 

=  3^0  +  6lti  +  9^2  +  12^4  +  2^5  +  3V‘6  +  6V8  +  6 ^9  +  6V'i 6  *h X  23 

VC24’K) 

=  5^)q  +  6^  +  9^  +  ^3  +  1°V4  +  2^5  +  ^6  +  2^7  +  ^Vs  +  ^9 

+  2^12  +  V13  +  V15  +  2Vi6  +  ^24 

Xq(C25»k) 

=  5(^o  +  +  lO^  +  2^3  +  8(j/^  -1-  4^  -i-  4^7  +  4yg  +  4y^  +  V^q 

+  4l4i2  +  2(44  +  ^25 

Xd(C26»k) 

=  5^o  +  ^Vi  +  -*-^2  +  2^3  +  IOY4  +  5y«j  +  V5  +  ^V7  +  2^3  +  ^9 

+  Vio  +  Vll  +  2Vl3  +  2^L4  +  Hi6  +  ^20 

^(C27’K) 

=  5^)q  +  6^_  +  10y2  +  2^3  +  10^  +  3^3  +  ^5  +  5(|^  -i-  4y8  +  4y^ 

%)(c28.k) 

+  ^11  +  2lVi2  +  2  ^3  +  Vl6  +  ^27 

=  6y0  •+■  +  7^2  +  7^4  +  ^3  +  7^3  +  2i|^  +  ^3  +  2iy^r,  +  2y^g 

+  V19  +V28 

<Xq(^'29  *  ^ 

=  6Vo  +  +  +  ^  +  ^8  +  ^t*17  +^29 

X0(c3o»k) 

=  6y0  +  +  8(^  +  8^ip  +  2y^  +  +  6y8  +  ^9  +  2^i6  +  Z|,^8 
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+  ^30 

^0^31 ’*0  =  +  8^2  +  9^4  +  2V5  +  V6  +  L%  +  ^9  +  2Vl6  +  v17 

+  2^9  +  V21  +  V31 

X0(C32,K)  =  6C|>0  +  6(|^  +  9(^2  +  8(A  +  ^5  +  V6  +  4V8  +  ^9  +  V10  +  2t^l6 

+  2vl8  +  2y2Q  +  (^2 

Xq(C^,K)  =  6(^0  +  ^  +  7y2  +  8^  +  +  z<Vg  +  3^  +  V]_6  +  2(^1? 

+  2^  +  ^33 

XqCc^.K)  =  6%  +  6^  +  Q(f>2  +  (^3  +  7^4  +  2  ^5  +  2(^7  +  z%  +  ^9  +  ^12 

+  ^13  +  +  V18  +  V19  '*'^34 

%0(C^,K)  =  6(^q  -1-  61^  +  8(^2  +  ^  +  6^4  +  2 +  (/?  +  +  L%  +  2^12 

2Vl8  +  ^21  +  V35 

fX0(c36«K)  =  6%  +  6l^l  +  +  ^3  +  ?^4  +  +  2(^7  +  2^8  +  %  +  V10 

+  Vi2  +  V14  +  2W.9  +  V20  +  ^36* 

X0(C3?,K)  =  6(^0  +  %  +  7(^2  +  ^3  +  31^4  +  l/5  +  V7  +  z%  +  V9  +  ^12 

+  2y17  +  yl8  +  y  37 

%0(C3g,K)  =  6ty0  +  6^  -I-  12t^2  -I-  1^3  +  6y4  -I-  10y5  +  3(^7  +  6^9  +  5^10 
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+  %4  +  2feo  +  V38 

X0^c39,k^  =  fyo  +  ^1  +  10<^2  *  ^3  +  8^4  +  3V5  +  ^  +  2V8  h'  7^9  +^10 
+  2V13  +  y14  +  y18  +  y20  +  y21  +  ^39 

Xo(c4o,k^  =  +  9V2  +  ^3  +  9(V4  +  ^5  +  ^7  +  3te  ■|'  6<t9  *H  3i//13 

+  31^  +  Vzk) 

0^(c4rK)  =  7%  +  %  +  5y2  +  4y4  +  3y8  +  4*17  +  V41 

XoCCie.K)  =  7g0  +  6yx  +  6^  +  6y4  +  y5  +  3y)8  +  3y9  +  3y19  +  y42 

'Xq(C4j.k)  =  ?y0  +  6y1  +  6y2  +  4y4  +  14  +  3y8  +  y^  +  2l+i7  +  Vis  +  V43 

%,(cft4.«  =  7Vo  +  6Vi  +  4*2  +  4*4  +  't'5  +  2’t'g  +  2'rl9  +  Vj.7  +  4  8  +  49  +l+Vt 

%0(C45-K)  =  4b  +  7Vl  +  13te  +  2V3  +  184  +  %  +  3V6  +  %  +  li4V8 
+V11  +  10V9  +  6Vl2  +  4Vl3  +  24.5  +  846  +  V23  +  '^24 
+  2V27  +  4-5 

'XqCc^.k)  =  5y0  +  ?V1  +  i5y2  +  3yij  +  isyi,.  +  8y5  +  3y^+  12V7  +  6Vg 

+  I2y9  +  2(^0  +  3(4!  +  6y13  +  6(44  +  6<46  +  y23  +  6y26 

+  Vi>6 

%0(C4?,K)  =  %)  +  ?Vi  +  li(V2  +  3V3  +  1?V4  +  %  +  +  1(¥?  +  1C% 


■ 
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+  lO^  +  ^10  +  2^n  -i-  6 ^2  +  ^Vi3  +  3(^4  -i-  V15  +  /jVi6 

+  2l|)24  +  +  2(^6  +  2^2?  +  y4? 

7.0(C48,K)  =  5(jJ0  +  +  15(^2  +  ^3  +  I8y4  +  ?y5  +  3y6  +  15y7  + 

+  +  ^LO  +  +  3yi2  +  6^L3  +  3^14  +  3^5  +  ^22 

+  ^'26  +  ^27  + 

OCqCc^.k)  =  6<^  +  7V-L  +  10^2  +  1^4  +  2^5  +  2y6  +  ioy8  +  8^  +  4yl6 
+  %7  +  &f>19  +  ^21  +  V29  +  2V31  +  4V33  +  V49 

X0(c30*k)  =  6y0  +  +  10^2  +  131fJ4  +  2V5  +  ^6  +  l4(^  +  6l^9  +  2(rt-5 

+  2^6  +  6^7  +  6(f^8  -1-  4ifij_9  +  4y28  +  y30  +  2(^3  +  y50 

ICqCc^.K)  =  6y0  +  7y>1  +  1^2  +  V3  +  12V4  +  %  +  4^6  +  2(^7  +  84fe 

+  8^  +  yi0  +  21^2  +  yi4  +  2l)46  +  4(|iL7  +  2yi8  +  4yi9 
+  2^0  +  V32  +  2^33  +  ^36  +  2^37  +  ^51 

7v(c52,k)  =  6y0  +  ?(+)- L  +  ioy2  +  y3  +  1^4  +  2^5  +  ^7  +  -^Ys  +  % 

+  21^22  +  y13  +  y15  +  8yL?  +  4yl8  +  2t^9  +  z^ZQ  +  2y34 

+  2V37  +  V32 

fX0(C33>K)  =  6(^0  +  7^  +  loy2  +  2V3  +  8^4  +  24^3  +  21^  +  8y8  +  4^ 
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+  4^  +  4V1?  +  4t^8  +^22  +  2(fj5  +  4^37  +  V53 
'"X0(C54,K)  =  6C|)q  +  7VX  +  12y2  +  2V3  +  ^4  +  5%  +  +  8(^  +  9ty 

+  V10  +  +  2(43  +  2(^L4  +  2  ^17  +  ^Vi8  +  Z^L9  +V'2  0 

+  ip21  +  iy25  +  2V34  +  +  ^36  +  V3?  +  y39  +  ^4 

""XqCC^.K)  =  6^j  +  7yi  +  1^2  +  2^3  +  ^4  +  1CV5  +  ^H7  +  ’l*  12(^9 

+  5y10  +  ^y12  +  6y14  +  4yL9  +  4y2Q  +y25  +  4^36  +  2^3g 

+  ^55 

rX0(C56 ,  K)  =  6  Y0  +  ?V1  +  12^2  +  Vj  +  14l+4  +  5^5  +  V6  +  ^7  +  10l^8 

+  11V9  +  V10  +  2^12  +  2Vl3  +  Vi4  +  4.5  +  3^16  +  2% 

+  3V18  +  3Vi9  +  2(fco  +  V21  +  V24  +  V28  +  V31  +  V32  +  V34 

+  ^36  +  *+39  +  V56 

X0(CJ7,K)  =  6%  +  7VX  +  1^2  +  Y3  +  WV4  +  3V5  +  %  +  2V7  +  ^ 

+  8y9  +  3^  +  V13  +  V15  +  3^6  +  3y17  +  ‘t^g  +  24'19 

+  ^21  +  ^24  +  V'28  +  V30  +  V31  +  V34  +  V35  +  ^37  +  ^57 

'X0(c58,k)  =  6^0  +  7y1  +  uy2  +  y3  +  14|)i  +  3Y5  +  Ve  +  3V7  +  1]V8 

+  8V9  +  2Y12  +  3y13  +  Vu  +  2^.6  +  4Hl7  +  2YL8  +  %9 
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‘^24  +  '^28  +  2^33  +  2^34  +  Vzj-0  +  V58 

^0^C59’K^  =  6Vo  +  +  -^-3^2  +  2Y3  +  +  +  Yg  +  6^7  +  8(^g  +  12  (^ 

+  V10  +  ^ll  +  ZV12  +  4Yl3  +  Viz*  +  2Vi6  +  4  ^8  +  2(^20  +  2^21 

+  2(^  +  +  2(^.  +  2(^2  +  ^9 

^0^c60’*^  =  %  +  *JVi  +  13^2  +  3(^g  +  15^4  +  6(j^  +  (^g  +  7^  +  7^g  +  12<^ 

+  V10  +  Vn  +  2Vi2  +  ^3  +  %4  +  2Yl6  +  Vl7  +  ^18  +  ^19 

+  ^20  +  V21  +  V26  +  ^27  +  %1  +  ^34  +  ^36  +  ^39  +  ^40  +  %0 

^(CgiiK)  =  6y0  +  7yx  +  15y2  +  2  ^3  +  1^4  +  11(^3  +  Yg  +  ^7  +  4{^  +  -*-4y^ 

+  5Vio  +  Yll  +  4W_3  +  6Vi4  +  2Yl6  +  29^8  +  4(^0  +  ^21 

+  21^5  +  Y32  +  2y38  +  2(^39  +  y6l 

0(o(C62,K)  =  6y0  +  7^  +  1^2  +  2^  +  15^  +  8(£  +  (^6  +  8^  +  6^  +  l4ty 

+  2YlO  +  t4l  +  %3  +  Z)Vl4  +  2Lti6  +  4Vi8  +  2(^20  +  2fel  +  2  (^6 
+  2  ^30  +  %9  +  ^62 

^(CgylC)  =  6  (^0  +  7^  +  12(^2  +  2^3  +  I4y4  +  4(^  +  <^g  +  ^  +  10^8 

+  +  4^12  +  4^3  +  2^  +  21^  +  4t^l8  +  2^9 

+  24*27  +  V30  +  4V3Zf  +  ^63 


%j(c64'k) 

XqCc^.K) 

'X0(c66,k) 

X0(c67.K) 

X0(c68,k) 

Xo^c69’*^ 

XqCc^q.K) 

X0^c7i 1  ^ 
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=  6y0  +  7y^  -i-  11^2  +  2V’3  "i'  12i44  +  3y^  +  yy  +  5y7  +  8y'g  + 

+  ^11  +  3^  +  ^13  +  Vi6  +  /fVL7  +  2^18  +  2l49  'i‘^27  +  V33 
+  ^34  +  2V37  ‘h  V64 

=  6y>Q  ■+  7^  -I-  121^2  +  2^3  ■+■  i3yy  +  5^  +  ^6  +  ^y7  **■  ^8  +  ®^9 

+  ^10  +  Yn  +  2^L2  +  2  ^3  +  2Vi4  +  Vl6  +  %7  +  Vl8  +  % 

-ny20  +V26  *h  V 33  ’h  2HJ34  +  2 ^36  +V65 

=  +  ^V]_  "**  ^^2  ■*"  3y^  +  31^  +  10yg  "*"  12^9 

+  8yl6  +  4yL?  +  6y19  +  2^  +  y23  +  z^31  +  2ify3  +  V66 

=  6(^q  +  7V-L  +  12^2  +  ^-^4  +  5y3  *1-  31^  +  l2tyg  +  I2y9  +  y^ 

+  9^5  +  9^3  +  3^0  +^23  +  3V3 0  +  ^32  +  V67 

=  7y0  +  7y-L  +  7(32  +  7(^4  +  7(jJq  +  7!y1?  +  7y4l  +y,68 

=  ?y0  ’*’  7Vi  ***  ^2  *"  ’l"  v3  ■*■  ®Vq  ■*"  2yy  ■*■  ■*■  2'Vi8  +  V 19 

+  ^29  *h  2^31  'h  2H 33  *  24*34  +l^69 
=  7\y0  +  71^  +  ay2  +  9^4  +  y5  +  91^  +  3y9  +  yy5  +  ^Vl?  +  2Vi8 

+  31f)19  +  V28  +  2V31  +  V32  'h  2^33  +  ^70 
=  7y0  +  7yx  +  8y2  +  9y4  +ly5  +1^6  +  6V8  +  ^ 9  +tVi6  +  4(^7 


+  ^18  -l-  2V19  +  2V’44  +V71 


^0^72*  ^  =  ^Vq  +  ^V]_  +  9^?  +  +  2X?  +  V5  +  8Vg  +  6V9  +  V-j_g  +  ^^7 

+  2Vi8  +  %  '*■  Y  21  +  2Hj33  +  ^43  +  V72 
^0(^73»^)  =  7iy0  +  ?H1  +  8Y?  +  HJ3  +  8  ^4  +  ^5  +  ^7  'H  ^8  +  V9  +  V12 

+  /^i7  "*■  V18  *  Y37  '**  2'^4i  "*"  Y  43  +  V  73 
^q(C74*k)  =  7Vo  +  7Vl  +  10V2  +  llvV4  +  H5  +  ^6  +  ^8  +  ^9  +  ^10 

+  2Hl6  +  2Vl7  +  2¥i8  +  Z¥l9  +  3^21  +  y32  +  Y33  +  2V44  ^  V  74 

XqCc^.k)  =  74)0  +  7^  +  9v2  +  11H)/|.  +  2y3  +  ^6  +  7^8  +  ^9  +  2Yi6 

+  3^]_2  +  ^18  +  34j_9  '*■  V21  +  X  31  +^33  +^4l  +^42  +  V44 

75 

X0(C76»k^  =  7Vo  +  7Yi  +  9^2  +  V3  +  ^4  +  %  +  V?  +  6Y8  +  2Y^  +^12 

+  4^7  +  4'^g  -1-  2^  +  2^3  +  Y76 

~X0(c77,k)  =  7y0  +  7^-l  +  9i|)2  +  10V4  +  2^3  +ty6  +  8y8  +  5iy9  +  2yl6  +  ^17 
+  Z^18  +  2Vi9  ■'■430  'b  X33  h  2^43  +  2Wj-  +X'77 

tX0(c78-k)  =  7Vo  +  7Vl  +  9Y2  +  V3  +  +  2Y5  +  Y7  -h  %  +  5Y9  +  2Yl2 

+  2Yl7  +  3^8  +  2Yl9  +  ^21  "h  V35  +  Y37  'H  2Y44  'b  Y78 
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^o(c79,k) 

=  7y0  +  7ty-]_  +  10^2  +  y3  +  8V4  +  ^  +  2V-^  +  5y8  -i-  5y^  +  y10 

+Vi2  'h(+'i4  +  2Vi 7  +  3yl8  +  2i^9  +y20  +  V36  +  2  ^44 

+V79 

X)(C80'K) 

=  7y0  +  71^-l  +  l0y2  +  9y4  +  ^  +  9y8  +  6y9  +  y1Q  +  y 13  -1-  Z*V17 

+  6lVi8  +  2  ^19  +  2^20  +  2^28  +  2^43  +  ^80 

Xo(Cgi,K) 

=  7y0  -1-  7y1  +  9y2  +  <y 3  +  8V4  +  2Y3  +  2Vy  +  6V8  +  3y^  +  V12 

+^13  +  4^7  +  2^3  +y19  +y34  +y3?  +y41  +  y43  +V44 

+V  81 

^•0^C82’K^ 

=  7y0  +  7y-]_  +  9y2  +  10yZf  +  2y^  +  9y8  +  6<y9  +  y^3  +  ^V17 

+  4ylg  +  4y19  +V21  *1-  2y28  +  y4l  +  2^44  +Y82 

X0(c83,k) 

=  7y0  +  7y-L  +  9y2  +  9y4  +  2y5  +  ioy8  +  4y9  +  y15  +  5^ 

+  6^18  +  2Hi9  +  2  ^2  8  +  2  ^43  +  2^Zj4  +V83 

X0(C84’K^ 

=  7y0  +  7y ^  +  ioy^  +  y3  -i-  ioy4  +  Z*V3  +  2y7  +  3y8  +  8y9  +  y10 

+  2^  +iyl4  +  6yL9  +  2y2Q  +  2yg6  +  2y42  +y84 

X0<VK) 

=  7y0  +  7^  +  9y2  +  y3  +  9y4  +  +  2y7  +  ?V8  +  z<y9  -i-  2V12 

+  V13  +  3yj_7  +  2l48  +  2^L9  +  2V34  +  V41  +  2H4  +  ^85 
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The  series  for  log  ZCC^K)  in  terras  of  <j>  •  j  <  i  are,  of  course, 
cal  to  the  above  series. 
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APPENDIX  E 


Symbolic  Equations  for  the  Face-Centred  Cubic  Lattice 
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Numerical  Valdes  of  Lattice  Constants 
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APPENDIX  F 


The  Triangular  Lattice 

In  the  notation  of  appendix  D,  the  zero  field  susceptibility 
series  is  given  by, 

a0  =  +2/3 

a1  s  +4 

q.^  =  +32 

9.j  =  +312 


The  Domb  ratio  plot  gives, 


a^  =  +3240 
a5  =  +31.584 
a6  at  +724,032 
a?  =  +34,031,232 


F.2 


and  the  Pade  approximant  table  for  Kc 


from  A.  log^(K) 


is  of  the  form. 


where  N.R.P,  =  no  real  root  and  N.P.R.  =  no  positive  real  root. 

The  manner  in  which  the  above  ratio  plot  contrasts  with  Figure 
6,1  and  the  complete  lack  of  consistency  of  the  Pade  approximant  table 
does  not  permit  us  to  draw  any  conclusions  regarding  the  existence  of 
a  critical  temperature  in  the  triangular  lattice. 
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